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The  processes  that  control  the  useful  life  time  and  performance  rates  in  non- 
flowing .zinc  electrodes  for  electrical  storage  batteries  have  been  investigated 
In  particular,  the  chemical  reaction  rates,  concentration,  and  electrical  po- 
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and  distributions  have  been  measured  and  correlated  for  RAI  P2291  ion  exchancex. 
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This  Final  Report  was  prepared  by  Won  G.  Sunu,  Manojit  Sinha, 
Siddhartha  Ghosh,  Bruce  M.  Eliash,  Douglas  N.  Bennion  (Principal  In- 
vestigator) and  John  Newman  (Co-Principal  Investigator) . The  report 
covers  work  done  April  15,  1976,  through  September  30,  1978,  on  con- 
tract number  F44620-76-C-0098.  Advice  and  direction  of  Dr.  David  Pickett 
and  other  personnel  at  AFAPL,  Wright  Patterson  Air  Force  Base,  is  grate- 
fully acknowledged. 

The  principal,  completed  results  of  this  work  are  reported  in  the 
Ph.D.  dissertation  of  Dr.  Won  G.  Sunu  and  the  M.S.  thesis  of  Mr.  Manojit 
Sinha.  Dr.  Sunu's  dissertation  is  included  in  its  entirety  as  a part  of 
this  report.  Three  papers  are  being  prepared  based  on  Dr.  Sunu's  re- 
search. They  will  be  submitted  for  publication  in  1979.  Mr.  Sinha's 
thesis  was  condensed,  modestly  revised,  and  published  in  the  J.  ElzcJtAa- 
diem.  Soc.  A reprint  of  that  paper  is  included  as  a part  of  this  report. 

Mr.  Siddhartha  Ghosh  was  supported  primarily  through  work  study 
funds  from  the  University  of  California.  Between  10  and  20%  of  his  sup- 
port was  from  this  contract.  His  thesis  is  included  as  part  of  this  re- 
port. A paper  based  on  that  work  is  to  be  submitted  for  publication  in 
1979. 

Two  projects  are  not  yet  complete.  Mr.  Bruce  Eliash  is  continuing 
the  work  of  Dr.  Sunu.  Mr.  Eliash' s work  should  be  finished  in  late  1979 
or  early  1980.  Publication  of  that  work  is  expected  in  1980.  Mr.  Mano- 
jit Sinha  has  begun  his  Ph.D.  research  work  on  the  performance  of  the 
nickel  oxide  electrode  in  alkaline  media.  A preliminary  theoretical 
model  is  proposed  and  a computer  code  has  been  written  which  allows  cal- 
culations based  on  the  model.  The  results  of  that  work  are  included 
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with  this  report.  Experimental  work  to  test  and  improve  the  model  is  be 
ing  started.  Publication  of  these  results  is  expected  in  1980. 
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ABSTRACT 


The  processes  that  control  the  useful  life  time  and  performance 
rates  in  non- flowing  zinc  electrodes  for  electrical  storage  batteries 
have  been  investigated.  In  particular,  the  chemical  reaction  rates,  con- 
centration, and  electrical  potential  distributions  have  been  predicted 
theoretically  and  measured  experimentally  as  functions  of  position  in 
the  electrode  and  time  for  both  charge,  discharge,  and  multiple  cycles. 
Changes  in  crystal  structure  and  solid  distributions  have  also  been  mea- 
sured and  examined  theoretically.  Similar  rates  and  distributions  have 
been  measured  and  correlated  for  RAI  P2291  ion  exchange  membrane  separa- 
tors used  in  zinc  secondary  batteries  such  as  the  nickel-zinc  battery 
(Zn/K0H,H20/Ni00H) . Distributions  in  flow  through  electrodes  were  also 
observed  experimentally  and  calculated  theoretically,  using  the  copper 
deposition  reaction  for  experimental  convenience.  Work  was  begun  on  the 
theoretical  description  of  the  nickel  oxide  electrode.  Experimental 
work  is  being  planned  to  validate  and  revise  our  understanding  of  how 
the  nickel  electrode  operates.  Work  is  continuing  on  measurement  of  con- 
centration profiles  in  zinc  electrodes. 

A number  of  factors  have  been  found  to  be  important  in  the  operation 
and  control  of  zinc  electrodes.  During  discharge,  zinc  metal  (Zn)  reacts 
with  hydroxyl  ions  (0H~)  to  yield  water  and  electrons  which  are  "pumped" 
into  the  external  electrical  circuit.  Initially  0H~  ions  within  the 
porous,  zinc  electrode  are  consumed.  When  those  OH"  ions  are  consumed, 
further  reaction  must  be  supported  by  transport  of  OH"  across  the  separa- 
tor from  the  counter  electrode,  nickel  oxide,  for  example.  For  RAI  P2291 

2 

membrane  separators,  it  was  found  that  reaction  rates  of  50  mA/cm  could 


v 


— - jj) 

2 

not  be  s imported  by  the  ion  transport,  but  20  mA/cm  could  be  supported. 

The  supply  of  OH  ions  across  the  separator  is  one  factor,  along  with 
the  electric  potential  distribution,  which  causes  the  discharge  reac- 
tion  to  concentrate  in  a narrow  zone  close  to  the  membrane.  At  20  mA/cn.' 
the  zone  is  about  0.3  mm  across.  At  50  mA/cm  the  zone  shrinks  to  about 

0.1  mm.  The  reaction  zone  thickness  has  important  implications  for  elec- 

trode design. 

During  charging  of  the  zinc  electrode,  OH"  ions  are  produced  within 
the  porous  electrode.  This  situation  causes  the  OH"  ion  concentration 
to  increase  in  the  depth  of  the  porous  electrode  favoring  a broad  reac- 
tion zone  and  lower  potential  losses  at  higher  current  densities  as  com- 
pared to  the  discharge  reaction.  The  implication  is  that  the  zinc  elec- 
trode can  be  safely  charged  at  higher  rates  than  discharged.  There  are 
other  considerations  which  enter,  however.  The  nickel  electrode  becomes 
OH  ion  "starved"  during  charging  and  wild,  dendritic  type  zinc  deposits 

get  worse  at  high  charge  rates.  Thus  high  charge  rates  are  not  really 

possible. 

Two  types  of  zinc  electrode  passivation  have  been  identified,  sur- 
face coverage  of  otherwise  reactive  zinc  area  and  pore  blockage.  Pore 
blockage  will  be  discussed  in  later  paragraphs.  The  exposed,  reactive 
zinc  is  gradually  covered  up  by  ZnO  precipitation  as  the  discharge  reac- 
tion proceeds.  At  modest  current  densities  and  reasonably  high  OH"  con- 
centrations, above  about  1.0  M,  the  precipitate  is  a porous,  fluffy  type 
deposit  which  slowly  inhibits  access  to  the  active  zone.  However,  this 
form  of  surface  coverage  allows  large,  useful  fractions  of  active  zinc 
to  continue  to  react.  It  appears  that  when  the  surface  concentration  of 
OH  ions  get  too  low,  the  form  of  the  precipitate  becomes  monolithic. 


vi 


I 

impervious,  and  totally  blanks  off  any  further  reaction.  The  electrode, 
or  at  least  the  affected  area  of  the  porous  electrode, is  passivated.  More 
detailed  observations  of  this  passivation  mechanism  are  needed;  but  the 
evidence  is  that  OH  ion  concentration  is  a significant  factor. 

Porosity  of  the  porous  zinc  electrode  is  shown  to  be  a critical 
parameter.  Zinc  oxide  (ZnO)  has  a much  larger  molar  volume  than  zinc. 

As  zinc  is  converted  to  zinc  oxide,  the  local  volume  fraction  of  solid 
increases  and  the  volume  fraction  of  electrolytic  solution  (KOH,K2ZnOH4, 

H20)  decreases.  As  has  been  shown,  the  reaction  rate  is  largest  at  the 
face  of  the  electrode  leading  to  pore  plugging  at  the  face  if  there  is 
not  sufficient  space  for  the  ZnO  to  precipitate.  When  the  face  of  the 
electrode  plugs,  further  discharge  is  impossible.  For  a porosity  (frac- 
tion of  volume  filled  with  electrolytic  solution)  of  0.33  pore  plugging 
at  the  face  was  the  failure  mode.  For  a porosity  of  0.6,  pore  plugging 
at  the  face  did  not  occur.  However,  tests  on  cycling  performance  showed 
that  with  increasing  number  of  cycles  more  and  more  zinc  and  zinc  oxide 
accumulated  in  the  high  reaction  rate  region,  that  is,  in  the  face  of 
the  electrode.  Long  term  cycle  tests  have  not  yet  been  done,  but  the 
clear  inplication  of  our  results  is  that  as  cycling  proceeds  zinc  trans- 
ports from  the  central  region  of  the  electrode  to  the  surface.  This 
transport  of  zinc  causes  increasing  porosity  in  the  central  region  and 
decreasing  porosity  at  the  face.  When  the  face  porosity  goes  below  a 
critical  value  of  about  0.5,  failure  by  pore  plugging  will  occur  and 
electrode  life  can  be  expected  to  end. 

The  0H~  concentration  distributions  have  been  shown  to  control 
many  factors  which  govern  electrode  performance.  Experiments  to  measure 
experimentally  the  actual  KOH  and  K^ZnOH^  concentrations  are  now  under 
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results  will  allow  validation  and  possibly  revisions  to  our  present  model 
of  transport  in  a porous  zinc  electrode. 

It  can  be  shown  that  once  zinc  electrodes  are  produced  which  approxi- 
mate the  practical  limit  implied  by  our  model  calculations,  the  nickel 
oxide  electrode  will  become  performance  limiting.  A theoretical  model 
of  the  nickel  electrode  has  been  completed  which  shows  the  interrelation- 
ships of  water  and  OH"  ion  transport  along  with  potential  and  reaction 
rate  distributions  and  their  relation  to  properties  of  the  nickel  oxide 
film  and  associated  charge  transport  steps.  Future  experiments  are 
designed  to  help  investigate  the  nickel  valence  states  which  can  be 
achieved  and  what  the  limitations  on  charge  transfer  capacity  are  for 
practical  electrodes.  In  addition  the  role  of  water  and  affects  of  water 
balance  on  the  nickel  oxide  electrode  will  be  considered. 

Flow  through  porous  electrodes  are  an  alternative  design  to  the  more 
customary  electrodes  with  limited  or  no  convective  flow.  In  convection 
designs,  the  working  electrode  is  inert,  and  the  ZnO  precipitates  in  a 
storage  tank  external  to  the  electrode  stack.  Experimental  and  theoreti- 
cal work  has  been  completed  demonstrating  the  interrelationships  between 
flow  rates,  applied  potential,  input  and  output  concentration,  electrode 
length,  and  electrode  thickness  with  potential,  current  and  reaction  dis- 
tribution and  cost  effectiveness.  The  experimental  work  was  done  with 
copper  deposition  at  600  ppm  input  concentration.  The  optimum  conditions 
were  found  to  be:  electrode  thickness  of  0.6  cm,  applied  potential  of 
1.4  V,  and  superficial  velocity  of  about  11  cm/s.  The  minimum  in  velocity 
was  very  broad. 


■ 


viii 


UNIVERSITY  OF  CALIFORNIA 


Los  Angeles 

Transient  and  Failure  Analyses  of  Porous 
Zinc  Electrodes 


A dissertation  submitted  in  partial  satisfaction  of  the 
requirements  for  the  degree  Doctor  of  Philosophy 
in  Engineering 


by 


Won^Guen  Sunu 
and 

Douglas  N.  Bennion 
Thesis  Dire' ior 


1978 


TABLE  OF  CONTENTS 


LIST  OF  SYMBOLS  

LIST  OF  TABLES  

LIST  OF  FIGURES  

ACKNOWLEDGMENTS  

VITA  

ABSTRACT  

I.  INTRODUCTION  

II.  THE  MATHEMATICAL  MODEL  

II-l.  Zinc  Electrode  Reaction  

II-2.  Ternary  Electrolyte  Theory  

II-3.  Description  of  the  Model  

II- 4.  Computation  Procedures  

III.  ZINC  ELECTRODE  BEHAVIOR  PREDICTED  FROM  THE  MATHEMATICAL 

MODEL  

III- l.  Discharge  Behavior  

III-2.  Behavior  on  Cycling  

III- 3.  Zinc  Electrode  Failure  Modes  

IV.  EXPERIMENTS 

IV-  1.  Preparation  of  the  Pressed  Zinc  Electrode  . 

IV-2.  Potential  Measurements  

IV-  3.  Determination  of  Reaction  Profiles  .... 

V.  RESULTS 

V- l.  Reaction  Profiles  

V-2.  Overpotentials  of  the  Zinc  Electrode  .... 

VI.  DISCUSSION 


76 


4 


VI-1.  Failure  During  Discharge  91 

VI-2.  Depletion  of  Electrolytes  92 

VI-3.  Pore  Plugging  95 

VI-4.  Passivation  99 

VI-5.  Failure  During  Charge  103 

VI-6.  Failure  During  Cycling . 105 

VII.  CONCLUSION  HO 

REFERENCES  U2 

APPENDIX  A 

A-l.  Ternary  Transport  Parameters  as  Functions  of 

Friction  Coefficients  120 

A-2.  Friction  Coefficients  as  Functions  of  Six 

Transport  Parameters  121 

A-3.  Transformation  of  the  Transport  Equations  into 

the  Other  Four  Component  Systems 122 

A-4.  Ternary  Transport  Parameters  Calculated  from 

Two  Sets  of  Binary  Data 12) 

* 

APPENDIX  B Mass  Transfer  Coefficients  (k^,  kg,  and  kg)  . 128 

APPENDIX  C 

C-l.  Solubility  of  ZnO  in  KOH  Solution  ....  133 

C-2.  Conductivity  of  KOH  Electrolytic  Solution  . . 134 

C-3.  Activity  Coefficients  of  Potassium  Hydroxide 

and  Potassium  Zincate 135 

APPENDIX  D 

D-l.  Input  Data  136 

D-2.  Programming  Parameters 139 

D-3.  Computer  Program 142 

iv 


LIST  OF  SYMBOLS 


2 

cross  sectional  area  of  zinc  electrode,  cm 

2 3 

specific  surface  area  per  unit  volume  of  the  electrode,  cm  /cm 

active  surface  area  for  charge  transfer  reaction  per  unit  volume 

2 3 

of  the  electrode,  cm  /cm 

active  surface  area  for  precipitation  or  dissolution  of  ZnO  per 

2 3 

unit  volume  of  the  electrode,  cm  /cm 

3 

concentration  of  species  1,  mole/cm 

3 

Initial  concentration  of  species  1,  mole/cm 

average  concentration  of  species  1 within  the  boundary  layer, 

3 

mole/cm, 

3 

equilibrium  or  saturation  concentration  of  zincate  ion,  mole/cm 

3 

concentration  of  species  1 at  the  active  metal  surface,  mole/cm 

thickness  of  the  solution  or  membrane  boundary  layer,  cm 

2 

binary  diffusion  coefficients,  cm  /sec 

2 

diffusion  coefficients  with  respect  to  solvent  velocity,  cm  /sec 

diffusion  coefficients  across  the  boundary  layer  (solution  or 
membrane) , cm^/sec 

Faraday's  constant,  96,487  coul/equlv. 

2 

superficial  applied  current  density,  amp/cm 

2 

superficial  current  density  in  matrix  phase,  amp/cm 

2 

superficial  current  density  in  solution  phase,  amp/cm 

jxchange  current  density  based  on  initial  concentrations  of  c° 
and  c°,  amp/cm^ 

2 

exchange  current  density,  amp/cm 

local  transfer  current  density  based  on  unit  area  of  active 
zinc  surface,  amp/cm^ 


v 


friction  coefficients,  interaction  of  species  1 and  j,  Joule- 
sec/cm^ 

mass  transfer  coefficients  of  potassium  zincate  from  active 
sites  to  bulk  in  the  pores  or  vice  versa,  cm/sec 

mass  transfer  coefficients  of  potassium  hydroxide  from  active 
sites  to  bulk  in  the  pores  or  vice  versa,  cm/sec 

rate  constant  for  precipitation  or  dissolution  of  ZnO,  cm/sec 

rate  constants  combining  mass  transfer  of  zincate,  k^,  with 
chemical  rate  constant,  k^^,  cm/sec 

electrode  thickness,  cm 

2 

multicomponent  transport  parameters,  mole  /sec-cm- Joule 

mass  transfer  coefficients  of  species  1 across  the  boundary 
layer  used  in  Table  I and  equation  (55),  cm/sec 

chemical  symbol  for  species  i 

molality  of  species  i,  mole/ Kg  of  water 

2 

flux  of  species  1 in  the  x-dlrection,  mole/cm  -sec 

2 

flux  of  species  1 across  the  boundary  layer,  mole/cm  -sec 

2 

flux  of  species  1 across  the  electrode  surface  (x**L) , mole/cm  - 


number  of  electrons  transferred  in  electrode  reaction 
number  of  moles  of  species  1 in  k th  sectioned  sample 
shape  factor,  dimensionless 

gas  constant  multiplied  by  absolute  temperature.  Joule/mole 

3 

source  term  of  species  1 in  equation  (35),  mole/cm  -sec 

stoichiometric  number  of  species  1 for  charge  transfer  reaction, 
defined  in  aquation  (10) 

stoichiometric  number  of  species  i for  precipitation  or  dissolu- 
tion of  ZnO,  defined  as  in  equation  (10) 


t,  t tortuosity  factor 

t°  transference  number  of  species  i referred  to  solvent  velocity 

t^  transference  number  of  species  i across  the  boundary  layer 

3 

partial  molar  volume  of  species  1,  cm  /mole 

V volume  of  electrolytic  solution  reservoir  between  the  membrane 

3 

and  the  zinc  electrode  surface,  cm 
v^  velocity  of  species  i in  solution,  cm/sec 

v^  superficial  volume  average  velocity,  SjC^V^v^  , cm/sec 
x distance  from  the  backing  plate,  cm 

z charge  number  of  species  i 


Greek  Symbols 

a kinetic  parameter 
ot^  a kinetic  parameter 

y exponent  in  zincate  concentration  dependence  of  the  exchange 

current  density 

y N mean  molal  activity  coefficient  of  species  i 

6 characteristic  diffusion  length  between  the  active  metal 

surface  and  the  bulk  in  the  pores,  cm 

e electrode  porosity 

£ ^ volume  fractions  of  solid  species  1 

C exponent  in  0H~  concentration  dependence  of  the  exchange 

current  density 

k conductivity  of  solution,  mho/cm 

k™  conductivity  of  solution  in  the  membrane,  mho/cm 


number  of  ionic  species  1 per  molecule  of  electrolyte  "k" 


vii 


v.  number  of  Ions  produced  by  the  dissociation  of  one  molecule  of 

* A A B B 

electrolyte  k.  v^“vi+  ^or  potassium  zlncate  and  VB*V1+V3 

■2  for  potassium  hydroxide, 
chemical  potential  of  species  1,  Joule/mole 
conductivity  of  solid  species  i,  mho/cm 
potential  in  metal  phase,  volts 
$2  potential  in  solution  phase  , volts 

g 

<f>2  potential  of  solution  at  the  active  metal  surface,  volts 
Subscripts 

1 any  arbitrary  species  1;  1 for  potassium  ion,  2 for  zlncate  ion, 

3 for  hydroxide  ion,  A for  potassium  zlncate,  and  B for  potass- 
ium hydroxide. 

o solvent 


LIST  OF  FIGURES 


Title  Page 

Geometric  features  of  the  porous  electrode 

(for  anodic  reaction)  19 


Overpotential,  (o)-  <J>_(L),  of  zinc  electrode  during 
constant  current ^discharge  for  type  I membrane  bound- 
ary (1-50  mA/cm^,  L-0.1  cm,  e°-0.4,  e?  -0.5,  e°  -0.1)  . 39 

zn  znu 

Concentration  profiles  during  constant  current  dis- 
charge for  type  I membrane  boundary  (1*50  mA/cm^, 

A for  one  minute  polarization  and  B,C,  and  D for  10Z, 

20%,  and  30%  depth  of  discharge,  respectively)  ....  40 

Concentration  profiles  of  KOH  after  10%  depth  of  dis- 
charge (1-50  mA/cnr  and  A,  B,  and  C for  solution 
boundary,  type  I membrane  boundary,  and  type  II  mem- 
brane boundary,  respectively)  42 

2 

Distribution  of  local  charge  transfer  current,  j(mA/cm  ), 
for  type  I membrane  boundary  (1-50  mA/cm^,  L-0.1  cm, 
and  A,B,C,and  D represent  0%,  10%,  20%,  and  40%  depth 
of  discharge,  respectively)  43 

Distribution  of  Zn  and  ZnO  plotted  as  volume  fractions 
(cm^/cm^)  for  type  I membrane  boundary  (1-50  mA/cm^, 
and  A,  B,  C,  and  D represent  10%,  20%,  30%,  and  40% 
depth  of  discharge,  respectively)  44 

Distribution  of  Zn  and  ZnO  plotted  as  volume  fractions 
(cm^/cmS)  for  solution  boundary  (1-50  mA/cm^,  and  A, 

B,  C,  and  D represent  10%,  20%,  30%,  and  40%  depth  of 
discharge,  respectively)  45 

Concentration  profile  during  constant  current  charge 

for  type  I membrane  boundary  (1-50  mA/cm^,  L-0.1  cm, 

and  A for  one  minute  polarization  and  B and  C for 

10%  and  20%  depth  of  charge,  respectively) 48 

Distribution  of  local  charge  transfer  current  density, 
jftaA/cnr),  during  the  first  cycle  for  type  I membrane 
boundary  (1-50  mA/cm^,  L-0.1  cm,  and  A,  B,  and  C re- 
present 0%,  10%,  20%  depth  of  discharge  or  charge, 
respectively)  49 

Overpotential,  #.(o)-  $_(L),  of  zinc  electrode  with 
time  during  three  cycles  for  type  I membrane  boundary. 

A,  1st  cycle;  B,  2nd  cycle;  C,  3rd  cycle  (1-50  mA/cmS 


x 


50 


11. 

12. 

13. 

14. 

15. 

16. 

17. 


18. 


19. 


20. 


21. 


L-0.1  cm,  e°“0 . 4 , s^-0.5,  e^-0.1) 


Distribution  of  Zn  and  ZnO  plotted  as  volume  fractions 
( cm3 / cm^)  during  two  cycles  for  type  I membrane  bound- 


lg  two  C’ 

ary  (1-50  mA/cm2,  L-0.1  cm,  e°  -0.5,  e°  -0.1)  . 

Lti  LtiU 

Distribution  of  Zn  and  ZnO  plotted  as  volume  fractions 
(cm3/cm3)  during  two  cycles  for  solution  boundary  (I- 
50  mA/cm2,  L-0.  1 co»  eZn“°*5’  EZnO"° * X) 

Distribution  of  Zn  and  ZnO  plotted  as  volume  fractions 
at  two  current  densities  for  type  I membrane  boundary. 
, after  one  cycle;  , after  two  cycles. 


Test  cell  assembly  

Electrical  circuit  for  operation  of  cell  .... 

Weight  percent  of  the  dissolved  zinc  in  1 M NH^Cl-NH^OH 
within  6 or  7 minutes  contact  time  


Profiles  of  Zn  and  ZnO  for  the  low  porosity  electrode 
discharged  to  8Z  depth  at  50  mA/cm2.  Q and  , mea- 
sured profiles  for  Zn;  histograms,  measured  profiles 
for  ZnO;  smooth  curves,  predicted  profiles  for  Zn  and 
ZnO  


Profiles  of  Zn  (circles  and  hexagons)  and  ZnO  (histo- 
grams) for  the  low  porosity  electrode  discharged  at 
20  mA/cm Q and  histogram  A,  measured  profiles  after 
10Z  depth;  Q and  histogram  B,  measured  profiles  after 
23Z  depth;  smooth  curves,  predicted  profiles  after  10Z 
and  23Z  depths  of  discharge  


Profiles  of  Zn  and  ZnO  for  the  high  porosity  electrode 
discharged  to  10Z  depth  at  50  mA/cm2.  Q and Q , mea- 
sured profiles  for  Zn;  histograms,  measured  profiles 
for  ZnO;  smooth  curves,  predicted  profiles  for  Zn  and 
ZnO  


Profiles  of  Zn  and  ZnO  for  the  high  porosity  electrode 
discharged  to  21Z  depth  at  50  mA/cm2.  O and  Q » mea" 
sured  profiles  for  Zn;  histograms,  measured  profiles 
for  ZnO;  smooth  curves,  predicted  profiles  for  Zn 
and  ZnO  


Profiles  of  Zn  (circles  and  hexagons)  and  ZnO  (histo- 
grams) measured  for  the  high  porosity  electrode  dis- 
charged at  20  mA/cm2.  O and  histogram  A for  10Z  depth 
of  discharge; Q and  histogram  B for  22Z  depth  of 
discharge  


51 

52 

53 
59 
61 

64 


68 


69 


71 


72 


73 


xi 


Profiles  of  Zn  (circles  and  hexagons)  and  ZnO  (histo- 
grams) measured  for  the  low  porosity  electrode  dis- 
charged at  50  mA/cm2.  Q and  histogram  A for  20Z  depth 
of  discharge  for  the  electrode  having  thickness  of 
0.05  cm;  Q and  histogram  B for  8Z  depth  of  discharge 
for  the  electrode  having  thickness  of  0.1  cm  . . . . 74 

Profiles  of  Zn  (top  line)  and  ZnO  (bottom  line) 
measured  for  the  low  porosity  electrode  after  one 
cycle.  At  each  half  cycle,  50  mA/cm2  and  10Z  depth  of 
initial  zinc  were  applied 75 


Overpotentials  at  the  face  of  the  zinc  electrode  (x-L) 
on  discharge  at  50  mA/cm2.  and  , measured  over- 

potentials; O,  predicted  overpotentials  for  L50-ZMD, 

% , predicted  overpotentials  for  L50-ZDM.  L and  H 
represent  low  porosity  and  high  porosity  electrodes, 
respectively,  and  Z,  M,  and  D represent  zinc  electrode, 
membrane,  and  dynel,  respectively  (Figure  notation 
L50-ZDM  represents  a low  porosity  zinc  electrode  haying 
a dynel  beneath  the  membrane  discharged  at  50  mA/cm2)  . . 78 

Overpotentials  at  the  face  of  the  zinc  electrode  (x**L) 
on  discharge  at  20  mA/cm2.  and  , measured  over- 

potentials; O and  #,  predicted  overpotentials  for 
L20-ZDM  without  and  with  corrections  for  swelling,  res- 
pectively. Explanation  of  figure  notation  is  given  in 
Figure  24  caption 80 

Photographs  of  the  electrode  surface.  A,  unused  zinc 
electrode;  B,  L50-ZMD  electrode  sfter  discharge  failure  . 82 


Photographs  of  the  electrode  surface  after  discharge 
failure.  A.  L20-ZMD  (low  porosity  electrode  discharged 
at  20  mA/cm2);  B,  H20-ZDM  (high  porosity  electrode 
discharged  at  20  mA/cm2)  


83 


Overpotentials  measured  at  four  different  positions  on 

discharge  of  L50-ZDM  electrode  at  50  mA/cm2.  Curves 

A,  B,  and  C represent  the  overpotentials  measured  at 

the  backing  plate,  at  the  face  of  the  test  electrode, 

and  above  the  membrane,  respectively.  Curve  D is  the 

potential  difference  between  the  zinc  test  and  the 

zinc  counter  electrodes  85 

Overpotentials  measured  at  four  different  positions  on 
charge  of  L50-ZDM  electrode  at  50  mA/cm2.  Explanation 
of  curves  A,  B,  C,  and  D is  given  in  Figure  28  caption. 

The  gaps  F-C  and  C-E  are  the  overpotentials  across  the 
the  membrane  and  across  the  solution  in  the  counter 
electrode  compartment,  respectively  88 


xii 


30 


33. 

34. 


Overpotentials  measured  at  the  £ace  of  the  electrode 
(L50-ZDM)  on  charge  at  50  mA/cm^.  Explanation  of  figure 
notation  Is  given  In  Figure  24  caption.  The  previous 
discharge  time  prior  to  onset  of  charge  is  39.4  minutes 
for  L50-ZMD,  69.6  minutes  for  L50-ZDM,  84.5  minutes  for 


H50-ZDM,  and  138.5  minutes  for  L50-ZD.  gg 

Photographs  of  the  dynel  cloth  after  one  discharge- 

charge  cycle.  Current  density  was  50  mA/cmt.  Zinc 

deposit  penetrates  the  dynel  towards  the  counter 

electrode.  A for  L50-ZDM  and  B for  L50-ZD 90 


Normalized  average  KOH  concentration  as  a function  of 
operation  time.  Curve  A is  for  one  layer  of  RAI  P2291 
membrane,  porosity  of  0.3,  and  electrode  thickness  of 
0.1  cm;  curve  B Is  for  two  layers  of  RAI  P2291  membrane, 
porosity  of  0.3,  and  electrode  thickness  of  0.1  cm; 
curve  C Is  for  two  layers  of  RAI  P2291  membrane,  poro- 
sity of  0.6,  and  electrode  thickness  of  0.1  cm  . . . .96 

Transference  numbers  of  ions  in  a 0.1  M HC1-KC1  ternary 

solution.  Data  points  are  the  measured  values  (67) 

and  solid  curves  are  the  calculated  values  127 

Diffusion  path  during  charge  near  the  solid  surface 
within  the  pores;  (A)  zinc  deposition  by  charge 
transfer  reaction  and  (B)  dissolution  of  ZnO  by 
chemical  reaction.  Dashed  lines  represent  the  direc- 
tion of  zincate  movement  during  charge  129 


ABSTRACT  OF  THE  DISSERTATION 
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A one  dimensional  mathematical  model,  perpendicular  to  the 

electrode  surface,  of  a porous  zinc  electrode  has  been  developed  on  the 

basis  of  concentrated  ternary  electrolyte  theory.  Numerical  techniques 

were  used  to  predict  the  galvanostatic  behavior  and  thus  the  failure 

mechanisms  of  porous  zinc  electrodes  under  conditions  similar  to  those 

of  secondary  zinc  batteries.  Profiles  of  zinc  and  zinc  oxide  and 

electrode  overpotentials  during  galvanostatic  operation  of  pressed  zinc 

electrodes  having  two  different  porosities  were  measured  at  superficial 

2 

current  densities  of  20  and  50  mA/cm  . A membrane  was  used  to  separate 

the  zinc  test  and  counter  electrode  compartments.  The  experimental 

observations  are  in  good  agreement  with  the  theoretical  predictions. 

2 

Discharge  at  50  mA/cm  caused  electrode  failure  due  to  the  deple- 
tion of  hydroxide  ions  within  the  zinc  electrode  compartment.  At  a 

2 

lover  current  density  of  20  mA/cm  , discharge  continued  until  pore 
plugging  (blockage  of  the  1st  kind)  or  passivation  (blockage  of  the 
2nd  kind)  caused  the  zinc  electrode  to  fail.  In  both  cases,  the  reac- 
reactlon  profiles  are  highly  nonuniform  and  the  reaction  zone,  located 
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near  the  electrode  surface,  is  very  thin.  On  repeated  cycling,  the 

| j 

difference  between  the  anodic  and  cathodic  reaction  distributions 

[ 

(caused  the  movement  of '.zinc  and  zinc  oxide  in  the  direction  perpendi- 
cular to  the  electrode  surface.  This  may  be  a limiting  factor  in  the 
losses  in  cell  capacity  of  secondary  zinc  batteries  which  are  designed 
to  eliminate  dendrite  penetration  and  shape  change. 


I.  INTRODUCTION 


The  zinc  electrode  has  been  commonly  used  as  a negative  electrode 
in  primary  cells  such  as  Leclanche  dry  cells  and  alkaline  zinc  manganese 
dioxide  cells.  It  also  has  received  considerable  attention,  during  the 
last  decade,  as  a negative  plate  in  secondary  batteries,  e.g.,  zinc- 
nickel,  zinc-silver  oxide,  zinc-bromine,  and  zinc-chlorine  cells.  The 
use  of  zinc  electrodes  in  secondary  batteries  is  attractive  because  of 
its  availability,  its  low  cost,  and  its  ability  to  deliver  high  energy 
and  power  density.  However,  secondary  cells  using  zinc  electrodes  have 
poor  cycle  life,  limited  by  zinc  electrode  failure.  Any  improvements 
on  the  cycle  life  of  the  zinc  electrode  can  be  applicable  to  its'use  in 
Zn-Ni  or  Zn-AgO  rechargeable  batteries  which  are  potential  applications 
for  use  in  electric  vehicle  propulsion  or  aerospace. 

Failure  of  the  zinc  electrode  can  be  attributed  to  a variety  of 
causes.  Two  major  causes  of  failure  are  zinc  dendrite  formation  and 
shape  change.  The  zinc  dendrites  which  are  formed  during  charge  pene- 
trate the  separator  and  cause  short  circuits  to  the  counter  electrode. 
Much  of  the  work  to  date  cn  improving  zinc  electrode  behavior  during 
charge  has  been  associated  with  zinc  deposit  morphology,  deposit  adher- 
ency,  and  mechanism  of  dendrite  growth  as  discussed  in  a review  by 
Oswin  and  Blurton  (1).  Zinc  dendrite  growth  appears  to  be  controlled 
by  the  mass  transfer  of  zincate  ions  (2)  and  is  characterized  by  a 
critical  current  density  (3).  The  preferential  dendrite  growth  at  the 
top  of  a porous  zinc  electrode  observed  by  Oxley  (4)  has  been  attributed 
to  zincate  depletion  of  the  electrolytic  solution  in  that  region  (1). 

The  results  of  these  studies  have  suggested  that  the  use  of  low  charge 
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rate,  pulse  charging,  and  metallic  or  organic  additives  have  some  bene- 
ficial effects  on  preventing  the  formation  of  dendritic  zinc  deposits 
(5).  The  detrimental  dendrite  penetration,  however,  has  been  success- 
fully reduced  by  use  of  an  appropriate  membrane  material  (6)  such  as 
commercially  available  grafted  polyethylene  ion  exchange  membrane. 

Another  failure  known  as  shape  change  is  the  redistribution  of 
solid  zinc  species  during  cycling  over  the  electrode  surface,  which 
eventually  decreases  the  effective  surface  area  and  thus  cell  capacity 
of  the  zinc  electrode.  Previous  investigators  (7,8,9)  have  indicated 
that  the  solid  species  moves  from  the  top  and  sides  to  the  bottom  and 
center  of  the  zinc  porous  electrode  during  cycling  and  that  loose  pack- 
ing, compared  to  tight  packing,  of  the  cell  promotes  more  rapid  redis- 
tribution of  zinc. 


. 


A mechanism  based  on  gravitational  effects  has  been  proposed  (10) , 
but  this  mechanism  is  not  consistent  with  the  fact  that  shape  change  is 
independent  of  cell  orientation.  Other  investigators  (7,9,11)  suggested 
concentration  cell  effects  as  a mechanism  of  shape  change.  McBreen  (7) 
has  suggested  that  a concentration  cell  is  generated  by  the  differences 
in  the  current  distribution  during  charge  and  discharge. 

Recently,  Choi  et  al.  ( 12,13)  have  proposed  that  the  main  reason 
for  the  shape  change  is  the  convective  flow  driven  by  membrane  electro- 
osmotic  effect  as  a dominant  factor.  The  shape  change  was  explained 
by  the  concentration  cell  caused  by  the  differences  in  the  convective 
flow  pattern  during  cycling,  i.e.,  flow  of  supersaturated  zincate  solu- 
tion towards  the  zinc  electrode  center  during  discharge  and  flow  of 
under saturated  zincate  solution  towards  the  electrode  edge  (or  solution 
reservoir)  during  charge.  The  agreement  between  their  experimental 
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results  and  theoretical  predictions  demonstrated  that  sealing  or  tight 
packing  of  the  cell,  preventing  convection,  successfully  eliminated  the 
shape  change. 

The  simple  approaches  mentioned  above  have  been  effective  in 
extending  cycle  life,  but  the  resulting  cycle  life  appears  still  to  be 
limited  by  other  factors.  Initial  loss  of  cell  capacity  as  well  as 
rapid  decrease  in  cell  capacity  with  continued  cycling  were  observed 
even  for  the  sealed  cells  for  which  chemical  analysis  for  zinc  species 
confirmed  no  shape  change  (13).  This  indicates  that  there  are  further 
reasons  for  cell  failure  to  be  explained  and  minimized  in  order  to  gain 
further  cycle  life  improvements.  A first  step  toward  disclosing  other 
failure  modes  and  mechanism  is  to  understand  and  be  able  to  predict  the 
zinc  behavior  within  the  porous  negative  electrode. 

A number  of  workers  investigated  extensively  the  anodic  zinc  behav- 
ior in  alkaline  solution,  using  polycrystalline,  flat  plate  electrodes. 
The  major  efforts  have  been  concentrated  on  the  studies  of  anodic  passi- 
vation of  the  planar  zinc  electrode.  It  has  been  generally  accepted 
that  the  applied  current  density  can  be  fitted  to  a linear  relationship 
with  the  inverse  square  root  of  the  passivation  time  (14,15,16,17) . 

A linear  relationship  based  on  the  Sand  equation  (18)  implies  that 
diffusion  of  zincate  ion  away  from  the  electrode  is  the  controlling 
process  for  passivation.  Therefore,  passivation  is  expected  to  occur 
when  the  layer  of  electrolytic  solution  adjacent  to  the  zinc  surface 
reaches  a critical  concentration  of  zincate  (16,17). 

Porous  electrodes  involve  much  more  complicated  and  interrelated 
phenomena  than  the  planar  electrode.  In  addition  to  passivation,  conver- 


sion  of  zinc  into  ZnO  during  discharge  of  the  porous  electrode  decreases 
the  pore  size  and  thus  the  transport  of  electrolyte  into  the  interior 
of  the  porous  electrode.  Some  workers  (19,20)  have  observed  that  utili- 
zation of  zinc  depends  strongly  on  the  initial  porosity  of  the  zinc 
electrode.  Membranes,  which  have  been  used  to  prevent  zinc  dendrite 
penetration,  can  limit  the  transport  of  electrolyte  species  to  and  from 
the  counter  electrode  as  was  discussed  by  Shaw  and  Remanick  (21). 

Porous  electrodes  have  highly  nonuniform  current  distribution  with  an 
effective  reaction  zone  which  is  very  thin,  resulting  in  local,  high 
current  density.  The  remaining  part  of  the  electrode  acts  as  an  inert 
matrix  which  does  not  contribute  to  the  cell  capacity.  Nonuniform  dis- 
tribution can  also  accelerate  the  decrease  in  pore  size  in  the  reaction 
zone.  In  order  to  understand  such  processes  and  further  identify  the 
failure  mechanisms,  it  is  helpful  to  develop  a mathematical  model  which 
describes  properly  the  transient  behavior  of  the  porous  zinc  electrode. 

Various  mathematical  models  have  been  proposed  to  describe  the 
behavior  of  the  porous  electrodes  as  summarized  by  Newman  and  Tiedemann 
(22) . Winsel  (23)  described  battery  discharge  by  considering  a single 
pore  electrode  of  constant  porosity  and  assuming  uniform  concentration 
and  conductivity  in  the  solution.  Newman  and  Tobias  (24)  included  the 
effects  of  a matrix  resistance  and  concentration  changes  in  their 
steady  state  analysis  of  a one  dimensional  homogeneous  mixture  model. 

An  attempt  to  consider  structural  changes  in  the  pores  was  initi- 
ated by  Alkire,  Grens,  and  Tobias  (25).  Using  a straight  pore  geometry, 
they  investigated  the  effects  of  changes  in  porosity,  concentration, 
and  fluid  flow  on  the  dissolution  of  a porous  copper  electrode  in 


sulfuric  acid  solution.  Dunning  et  al.  (26)  proposed  a mathematical 
model  which  predicted  discharge  and  cycling  behavior  of  a Ag/AgCl  elec- 


trode in  NaCl  solution  and  a Cd/Cd(OH)2  electrode  in  KOH  solution.  They 
considered  the  effects  of  complexing  of  the  sparingly  soluble  salt  with 
the  bulk  electrolyte  anion  and  the  detailed  structural  effects  on  mass 
transfer  parameter  on  the  limiting  species.  However,  the  concentration 
derivative  term  was  neglected  and  they  applied  a so-called  pseudo  steady 
state  approximation  to  describe  the  transient  behavior.  Recently,  Gu  et 
al.  (27)  extended  the  model  of  Dunning  et  al.  for  Ag/AgCl  electrode  to 
include  surface  diffusion  of  the  adsorbed  complex  ion  and  more  realistic 
active  surface  area  based  on  their  experimental  measurements.  The  tran- 
sient behavior  of  the  model  was  solved  without  using  the  pseudo  steady 
state  approximation. 

Simonsson  (28,29)  applied  his  model  to  predict  the  failure  mode  of 
of  the  positive  electrode  of  the  lead  acid  battery.  He  showed  that  the 
discharge  of  the  positive  lead  dioxide  electrode  at  high  current  load 
failed  by  the  severe  depletion  of  sulfuric  acid  at  the  pore  mouth,  while 
Gidaspow  and  Baker  (30)  emphasized  the  pore  plugging  as  a major  cause  of 
battery  failure.  A mathematical  model  which  predicts  the  transient 
behavior  of  the  porous  zinc  electrode  has  not  been  reported  as  yet 
except  that  of  Choi  et  al.  (12).  The  results  of  their  model  which  in- 
cludes a convection  term  in  the  y-direction,  parallel  to  the  electrode 
surface,  successfully  predicted  the  mechanism  of  the  failure  called 
"shape  change". 

In  the  present  work,  a mathematical  model  is  developed  to  describe 
the  transient  behavioT  of  the  porous  zinc  electrode  in  the  direction 


r ^ 


perpendicular  to  the  electrode  surface.  It  also  describes  failure 
mechanisms  during  discharge  and  cycling  operation  other  than  shape 
change.  The  porous  electrode  is  treated  as  the  superposition  of  two 
continue,  solution  and  solid  phases,  by  use  of  a macroscopic  averaging 
technique  (22) . This  model  resembles  real  electrodes  when  particle 
and  pore  dimensions  are  on  the  order  of  one  micrometer. 

The  products  produced  during  discharge  of  the  Cd/Cd (0H)2,  Ag/AgCl, 
and  Pb02/PbS04  electrodes  are  believed  to  be  sparingly  soluble  in  the 
corresponding  electrolytic  solutions.  For  the  zinc  electrode,  the  dis- 
charge  product  ZnO  is  highly  soluble  in  KOH  solution.  The  electrolytic 
solution  generally  used  for  secondary  batteries  using  zinc  negative 
electrodes  is  highly  concentrated  KOH  solution  (about  10M)  saturated 
with  ZnO  (about  1M) . Furthermore,  the  electrolytic  solution  in  contact 

[ 

with  a discharging  zinc  electrode  can  become  supersaturated  with  zincate 
(31),  resulting  in  the  concentrated  potassium  hydroxide-potassium 
zincate-water  ternary  mixture.  Therefore,  fundamental  equations  neces- 
sary for  describing  the  concentrated  ternary  solution  are  developed 
based  on  irreversible  thermodynamic  principles  (32) . The  resulting 
equations  can  also  be  applied  to  other  multi-component  systems  such  as 
molten  salt  and  membrane  transport. 

Commercial  zinc-nickel  or  zinc-silver  oxide  batteries,  in  general, 
have  separators  between  the  two  electrode  compartments  which  affect 
significantly  the  concentration  variations  inside  the  porous  electrodes. 
The  effect  of  the  membrane  separator  and  the  electrolytic  solution 
reservoir  on  the  zinc  electrode  performance  are  included  in  the  model 
to  predict  the  mathematical  behavior  under  conditions  similar  to  a 
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practical  battery  system.  The  mathematical  model  presented  here  con- 
sists of  Ohm's  law  in  the  solution  phase.  Ohm's  law  in  the  matrix  phase, 
kinetic  rate  expression,  conservation  of  charge,  and  conservation  equa- 
tions of  hydroxide  ion,  zincate  ion,  water,  and  solid  species.  The 
effect  of  local  mass  transfer  between  the  electrode  active  surface  and 
the  bulk  solution  in  the  pores,  changes  in  porosity,  variations  of 
active  surface  area,  and  convection  arising  from  the  differences  in 
molar  volumes  of  reactants  and  products  are  included.  The  results  of 
the  numerical  solution  are  used  to  suggest  certain  modes  of  failure  and 
performance  limitations  of  the  zinc  electrode. 

Experimental  measurements  of  zinc  electrode  overpotentials  and 

current  distribution  were  carried  out  using  galvanostatic  operation  at 

2 

two  different  current  densities;  20  and  50  mA/cm  . Those  results  were 
compared  with  the  mathematical  model  predictions  for  zinc  electrode 
failure. 

There  are  several  alternate  procedures  for  fabricating  porous  zinc 
electrodes.  They  include  pressing  of  ZnO  powder  mix  (ZnO,  HgO,  and 
some  other  additives)  which  gives  the  zinc  electrode  in  the  discharged 
state  (33) , pressing  of  amalgamated  zinc  powder  (34) , metal  spraying 
process,  and  electrodeposition  process  which  result  in  the  zinc  elec- 
trode in  the  charged  state.  Determination  of  failure  modes  during  dis- 
charge and  cycling  is  the  primary  objective  of  this  work.  The  porosity 
and  volume  of  electrolyte  reservoir  of  the  electrode  should  be  carefully 
controlled  during  preparation  in  order  to  obtain  reproducible  electrode 
performance  and  failure  mechanism.  The  pressing  techniques  developed 
by  Morrell  and  Smith  (34)  satisfy  most  of  the  above  mentioned  require- 


ments  and  yield  an  electrode  in  the  charged  state. 

There  are  three  approaches  for  the  experimental  determination  of 
current  distribution.  The  first  approach  is  the  direct  measurement  of 
potential  distribution  by  using  numerous  reference  electrodes  along  the 
electrode.  This  approach  was  employed  by  Brodd  (35)  in  the  study  of 
current  distribution  in  a manganese  dioxide  electrode.  The  second 
approach  is  the  direct  measurement  of  the  current  distribution  by  using 
a sectioned  electrode  and  measuring  the  current  in  each  section  as  was 
used  by  Coleman  (36)  for  Mn02  electrode  and  by  Gagnon  and  Austin  (37) 
for  Ag/Ag20  electrode.  The  third  approach  is  to  measure  the  reaction 
distribution  of  the  discharged  electrode  by  electrode  sectioning  fol- 
lowed by  the  analysis  of  each  section.  The  charge  state  of  each  section 
can  be  analyzed  using  chemical  analysis  as  shown  by  Nagy  and  Bockris 
(38)  for  the  zinc  electrode,  Bro  and  Kang  (39)  for  the  cadmium  electrode, 
and  Alkire  (40)  for  the  copper  electrode. 

The  first  two  direct  approaches  have  the  advantage  of  allowing 
continuous  measurements  during  operation,  however,  reference  electrode 
probes  or  a discontinuous  sectioned  matrix  may  disturb  the  real  current 
distribution  and  can  lead  to  distorted  results.  The  third  approach  is 
limited  to  one  experimental  data  set  per  electrode,  and  the  sectioning 
procedure  can  result  in  a net  reduction  and  redistribution  of  material. 
The  sectioning  approach  seems  to  be  one  of  the  promising  experimental 
techniques  required  for  battery  development  since  its  modification  (so 
called  freezing  technique)  allows  determining  the  distribution  of  both 
solid  and  liquid  species.  By  freezing  the  electrolyte  in  the  boundary 
layer,  followed  by  sectioning  and  chemical  analysis,  Brenner  (41)  and 
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Flatt  et  al.  (42)  successfully  determined  the  concentration  profiles  in 
solution  at  dissolving  anode  surfaces  and  at  a depositing  cathode  sur- 
face, respectively.  In  the  present  study,  zinc  and  zinc  oxide  profiles 
are  measured  using  the  sectioning  method  with  EDTA  titration.  Freezing 
was  not  tried. 


i 


j 


Zn(OH)"  ♦ OH"  — — 

» Zn  (OH)  ~ ♦ e" 

(G) 

Zn  (OH)  ~ + OH" 

==  Zn(OH)* 

(H) 

where  reaction  step  (G)  is  the  rate  controlling  step. 

In  a concentrated  KOH  solution  saturated  with  ZnO,  it  is  possible 
for  the  following  complexes  of  zinc  species  to  be  present  as  intermedi- 
ate compounds:  Zn++,  Zn(OH)+,  Zn(OH)2,  Zn(OH)j,  Zn(OH)*  ...  Zn(OH)^~n. 
These  complexes  are  formed  according  to  the  following  equation 

Zn+%  nOH~—  — Zn(OH)^~n  (I) 

4 15 

The  equilibrium  constants  reported  by  Butler  (53)  are  10  ‘ for  n»l, 
IqIO.15  £0r  1q14.25  ^ 1015  for  n»4.  Using  these  values, 

Boden  et  al.  (54)  concluded  that  the  only  species  present  to  any  appre- 
ciable extent  in  KOH  solution  was  Zn(OH)*  ion.  Other  workers  (55,56,57) 
also  confirmed  that  the  predominant  species  was  zincate  ion.  Therefore, 
in  the  present  work,  the  electrolytic  solution  is  considered  as  a four 
component  system  containing  K+,  OH",  Zn(OH)*,  and  H20. 

II-2.  Ternary  Electrolyte  Theory 

It  is  shown  that  the  solution  under  consideration  is  a ternary 
mixture  containing  potassium  hydroxide  as  well  as  potassium  zincate. 

The  solution  for  use  in  alkaline  zinc  batteries  is,  in  general,  highly 
concentrated  KOH  solution  (about  40  wt%)  saturated  with  ZnO.  The  ZnO 
is  highly  soluble  in  KOH  solution,  and  it  is  generally  agreed  that 
anodic  discharge  of  zinc  in  KOH  solution  produces  a supersaturated 
zincate  solution  which  is  often  greater  by  a factor  of  two  than  the 
equilibrium  zincate  concentration  (31). 
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As  Che  concentration  of  an  ionic  species  become  large,  interactions 
between  ions  and  solvent  and  ion-ion  interactions  become  important. 
Newman  et  al.  (58)  have  developed  a binary  electrolyte  theory  for  con- 
centrated solution.  The  following  analysis  of  a four  component  system 
will  parallel  that  of  a binary  electrolyte  considered  by  Newman  et  al. 

The  general  concept  of  relating  a driving  force  per  unit  volume 
on  a species  i,  d^C-c^Vu^),  to  a linear  sum  of  frictional  interactions 
exerted  on  species  i by  species  j,  (v^-v^ , was  proposed  by  Maxwell, 
often  referred  to  as  the  Stefan-Maxwell  equation  (59). 

di  ■ ci  71,1  ■ ]l  kij  <VJ  - V (1) 

J 

where  v.^  is  the  velocity  of  species  i,  c^  is  the  concentration  of  spe- 
cies i,  and  is  the  chemical  potential  or,  in  case  of  ionic  species, 
is  to  be  interpreted  as  electrochemical  potentials  as  suggested  by 
by  Guggenheim  (60).  The  friction  coefficients  are,  in  general, 
quite  concentration  dependent.  However,  they  can  be  defined  in  terms 
of  binary  diffusion  coefficients  2^ , which  are  found  to  be  less 
condensation  dependent  (58,59) 


RT  ciej 


cT®ij 


where  c^  is  the  total  concentration  of  the  solution.  From  Newton  s 
Third  Law  of  Motion,  Onsager's  reciprocal  relation  is  satisfied  (61); 


w Therefore,  the  number  of  independent  friction  coefficients 
for  the  n component  system  are  n(n-l)/2. 

For  the  four  component  system  under  consideration,  i.e.,  common 
cation  K+  (species  ”1”),  anions  Zn(OH)^  (species  "2")  and  0H~  (species 
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"3"),  water  (species  "o") , the  multicomponent  transport  equation  (1) 
can  be  represented  as  follows: 


c2  vu2  ■ Sl'vrv2>  * *2S(Vt2)  * K2o(VV2> 


lo'VV 

(3) 

2o(Vv2> 

(4) 

30(VV3> 

(5) 

The  fourth  equation  for  water  is  not  an  independent  equation  due  to  the 
Gibbs -Duhem  equation,  SjCjVu^O. 

The  chemical  potentials  of  electrolytes  A and  B are  defined  by 
the  following  equations. 

% » v*  Viij  ♦ v*  Vy2  (6) 


% ■ A lvl  * A 7u3  (7) 

A 

where  the  electrolyte  A (K2Zn(OH)^)  contains  Vj  number  of  species  "1" 
(K+)  and  v*  number  of  species  "2"  (Zn(OH)*).  Likewise,  v®  and  v®  are 
defined  by  the  number  of  species  "1"  (K+)  and  species  "3"  (OH-),  res- 
pectively, per  molecule  of  electrolyte  B (KQH) . Electroneutrality 
A A B B 

requires  ZjVj+ZjVj^O  and  z^v^*ZjVj»0  (or  zici+z2c2+z3C3*0^  * 

The  current  density  is  written  as 


i « F £ z^  ■ FZjCjVj  ♦ Fz2c2v2  * Fz3c3V3 


If  the  electrical  potential  ♦ is  defined  by  the  potential  of  a 
suitable  reference  electrode  at  a point  in  the  solution  measured  with 
respect  to  a similar  reference  electrode  at  a fixed  point  in  the  solu- 
tion, thermodynamic  principles  yield 


X'  ...»  * 


- nF  **  ■ sx  Vlij  * s2.Vu2  ♦ s3  Vy3  + sq  VyQ 


(9) 


The  coefficients  are  the  stoichiometric  coefficients  of  species  i 

for  the  electrode  reaction  represented  by 

z. 


E.  si  Mi 


ne 


(10) 


where  NT  is  a symbol  representing  the  chemical  formula  of  species  i. 

Equations  (3)  to  (8)  are  six  linear  equations  containing  ten  varia- 
bles; and  i*  rearran9ement  equations 

(3)  to  (8)  in  terms  of  fluxes  N^(or  c.v^)  is  straightforward  but  quite 
involved.  The  above  system  of  four  components  is  very  similar  to  those 
considered  by  Bennion  (membrane,  single  salt,  and  water) (62) , by 
Tiedemann  and  Bennion  (cation,  anion,  neutral  salt,  and  water) (63) , and 
by  Miller  (two  salts  with  common  cation  and  water) (64) . 

Using  the  reference  velocity  as  the  solvent  velocity  vq,  the 
resulting  rearrangement  through  matrix  inversion  is  given  by  the  follow- 
ing equations. 


\ \ t° 

N1  - - (^AA+  l!LAB  H-  ( 7;  1 + C1 

N2  ‘ * V^AA  % ' V2**AB  % * "7  1 * c2Vo 

*2F 

B B t0x 

N3  * • v3LAB  % - v3LBB  % * , 1 * c3vo 

s3F 


N ■ c v 
0 00 


(11) 

(12) 

(13) 

(14) 


The  arrangement  of  Equation  (9)  for  potential  coupled  with  equa- 
tion (11)  through  (14)  results  in  the  following  equation. 
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transport  parameters  are  also  shown  in  Appendix  A- 2.  Appendix  A- 3 

shows  how  the  equations  (11)  to  (15)  can  be  transformed  to  represent 
transport  equations  for  the  other  systems  such  as  single  salt-membrane- 
water  and  single  salt-neutral  salt-water. 

The  flux  equations  containing  the  diffusion  coefficients,  which 
are  usually  reported,  are  expressed  in  terms  of  concentration  gradients. 
The  chemical  potential  gradient  is  related  to  concentration  gradients 
as 


3cA  dx  3cB  dx 


(16) 


The  solvent  velocity  vq  has  been  chosen  as  a reference  velocity 
in  the  above  treatment.  In  some  cases,  such  as  in  porous  electrodes, 
it  is  more  convenient  to  use  volume  average  velocity  v^  as  a reference 
velocity.  The  volume  average  velocity  is  defined  by 


* nVi 


ci7i 


Z 

i 


where  V *s  t*ie  Part^a^  ®°lar  volume  of  species  i. 


IS 


" , > v ' t • 


(17) 


The  flux  equations  referred  to  the  volume  average  velocity,  using 


concentration  gradients,  are  shown  below: 


N2  dcA 

_£  « - D — — 

vA  Mdx 


D i ♦ c vD 


^3  „ ^ n ^ + t3  . □ 

B " " °BA  . " °BB  . B 1 CB 

v3  dx  dx  z3v3F 


n dcA  dcB  * * n 

0 » ♦ V — + V — - Cotn  ~ ♦ c vu 
0A  dx  08  dx  ° 0 F 0 


where 


V*«tV  ! dm-l 

o 


AA  AA  3cA  “AB  3ca 


CA 
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AB  AB  c OB 
o 


,0  CB 


AB  AA  3cg  “AB  3cb 


dba-dba-^doa  ' D 

o 


BA  " UAB  3ca  BB  3ca 


Ddo  - D°  - D ; D°  - L 


BB  BB  c OB 
o 


BB  “AB  3cB  “BB  3cb 


°0A  * Co^ALAA+  VBLAB)  3ca  + Co(VALAB+  VBLBB)  3ca 
°0B  * co^ALAA*  ^BLAB5  3c^  + Co(\LAB+  ^BLBB}  3c^ 


t2  * *2  ‘ *2C2to 
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* Vi  Vl  v< 


16 


For  the  solvent  velocity  reference  frame,  the  parameters  Dq^,  dob* 

* 

and  tQ  disappear  from  the  above  equations.  The  quantities  , t^,  and 
v^ should  be  replaced  with  D?^ , t?,  and  vq,  respectively,  to  represent 
the  flux  equations  referred  to  the  solvent  velocity  vq. 

* 

The  partial  molar  volumes  of  the  ions  shown  in  the  parameter  tQ 
satisfy 


VA-  '’iVl 

* v2  ^2 

(22) 

V 

+ vB  7 

V3  3 

(23) 

These  relations  are  not  sufficient  to  determine  the  partial  molar 

volumes  of  ions  separately.  Consequently,  we  let 
o o 

V J-V-Vb  (24, 

V1  3 


This  choice  is  based  solely  on  convenience  (65),  and  no  physical  signi- 

* 

ficance  should  be  attached  to  it.  With  this  choice,  tQ  becomes  zero, 
and  in  equations  (18)  and  (19)  reduce  to  t°. 

The  expressions  of  Appendix  A-l  enable  us  to  calculate  the  trans- 
port parameters  for  a ternary  mixture  from  the  data  of  binary  electroly- 
tes. Tabulated  data  for  (or  Q ^)  for  1:1  and  1:2  binary  elec- 
trolytes have  been  compiled  by  Chapman  (66) . The  transference  numbers 
and  conductivities  calculated  from  the  two  sets  of  binary  data  showed 
good  agreement  with  the  experimental  data  for  the  ternary  mixtures. 

An  example  is  shown  in  Appendix  A-4  where  measured  transference  numbers 
for  HC1-KC1  solution  (67)  are  compared  with  the  transference  numbers 
estimated  from  the  binary  data  (66) . 


However,  successful  prediction  of  the  diffusion  coefficients,  D. ., 


for  ternary  mixtures  requires  the  knowledge  of  the  ternary  activity 
coefficients  as  indicated  by  equations  (21) . The  activity  coefficients 
for  a mixed  solution  may  be  approximated  by  use  of  Newman's  treatment 
(68).  Some  workers  (64,69,70)  have  used  experimental  ternary  activity 
coefficients  to  show  reasonable  agreements  between  the  measured  and  the 
predicted  diffusion  coefficients.  They  applied  the  multicomponent 
transport  equations  of  Miller  (64),  which  were  derived  based  on  an 
inverse  description  of  the  equation  (1)  of  the  form  : Jj*  Z^l^X^  where 
Xj  is  the  thermodynamic  driving  force,  is  the  molar  flux  of  species 
i,  and  1^  are  the  phenomenological  or  transport  coefficients.  The 
flux  equations  (18) , (19) , and  (20)  and  current  equation  (15)  will  be 
used  to  describe  the  mathematical  model  of  the  porous  zinc  electrode. 

II-3.  Description  of  the  Model 

A one  dimensional  mathematical  model  of  a porous  zinc  electrode 
has  been  proposed  which  predicts  OH~  concentration  c^,  Zn(OH)*  concen- 
tration c2,  superficial  current  density  in  the  solution  i2,  transfer 
current  density  j,  porosity  e , superficial  volume  average  velocity  vR 
potential  in  the  solution  d»2 , and  potential  in  the  matrix  ^ as  a 
function  of  time  and  position  perpendicular  to  the  surface  of  the  elec- 
trode. The  geometric  features  of  the  model  are  shown  in  Figure  1.  A 
backing  plate  is  located  at  x-0  and  the  membrane  or  external  boundary 
layer  is  at  x«L.  The  effects  of  membrane  as  well  as  volume  of  solu- 
tion reservoir  between  the  electrode  surface  and  membrane  on  the  elec- 
trode performance  were  included  to  predict  zinc  electrode  behavior  un- 
der conditions  similar  to  a practical  battery  system.  The  porous 
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electrode  is  treated  as  a homogeneous  mixture  of  the  ionically  conduct- 
ing solution  phase  denoted  by  subscript  "2"  and  an  electronically 
conducting  matrix  phase  denoted  by  subscript  "1".  The  solid  matrix 
might  contain  zinc,  zinc  oxide,  and  inert  conducting  materials  (for 
example,  graphite)  or  any  combinations  of  the  three.  The  equations 
developed  in  the  previous  section  can  be  applied  only  to  the  pure 
homogeneous  solution  phase  (porosity  of  unity) . In  the  porous  electrode, 
effects  of  porosity  and  tortuosity  on  the  transport  equations  must  be 
included. 

The  macroscopic  averaging  technique  performed  over  the  volume  of 
pores  within  the  electrode  has  been  described  in  detail  by  Newman  and 
Tiedemann  (22)  and  Diinning  (71) . This  macroscopic  averaging  treatment 
is  applied  to  the  present  model  by  disregarding  the  actual  geometric 
details  of  the  pores.  For  example,  consider  c^  as  the  concentration  of 
species  i averaged  over  the  volume  of  solution  in  the  pores.  Then  the 
superficial  concentration,  averaged  over  the  unit  volume  of  the  elec- 
trode, is  tc^  where  £ is  the  porosity  of  the  electrode.  If  is 
defined  by  the  superficial  flux  based  on  the  unit  cross  sectional  area 
of  the  electrode,  N^/e  represents  the  flux  through  the  unit  area  of 
solution  phase  in  the  pores.  Similarly,  v^c  is  the  velocity  in  the 
solution  phase  if  v^ is  a superficial  bulk  fluid  velocity. 

The  use  of  averaged  quantities  (Nj/e  , i2A  > and  vR'e  ) into  equa- 
tions (18),  (19),  and  (20)  yields  the  following  equations. 
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(26) 
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where  t?  is  the  transference  number  of  species  i,  t is  the  tortuosity 
factor,  c^  is  the  bulk  concentration  of  species  i averaged  over  the 
volume  of  solution,  and  H^,  i^,  and  v^ refer  to  the  superficial  quanti- 
ties based  on  the  unit  cross  sectional  area  of  both  matrix  and  pores. 

The  diffusion  coefficients  D.^  are  replaced  with  e1+t  to  account  for 
the  tortuosity  correction.  Schofield  and  Dakshinamurtl  (72)  have 
suggested  that  the  effective  diffusion  coefficients  are  proportional 
to  (or  t-0.5). 

The  basic  equations  of  the  model  to  be  discussed  are  Ohm's  law  in 
the  solution  phase,  redox  kinetic  expression,  conservation  of  zincate, 
conservation  of  hydroxide  ion,  overall  conservation,  and  Ohm's  law  in 
the  matrix  phase. 


Ohm's  Law  in  the  Solution  Phase 


A modified  Ohm's  law  applied  to  the  ternary  electrolyte  solution 
in  the  porous  electrode  is  obtained  from  equation  (15) . 
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where  d>2  is  the  potential  of  a reference  electrode  of  the  same  kind  as 
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the  working  electrode  placed  in  the  bulk  solution  within  the  pores  of 
the  porous  electrode;  < is  the  conductivity  of  solution  in  the  pores  at 
concentration  c^;  c^  is  the  concentration  of  species  i where  subscript 
o refers  to  solvent;  ^ is  the  chemical  potential  of  species  i;  and  s^, 
n,  and  z^  have  meanings  as  defined  in  equation  (10) . The  effective 
solution  conductivity  corrected  for  tortuosity  considered  here  is  repre- 
sented by  ^+t.  Meredith  and  Tobias  (73)  and  Gagnon  (74)  suggested 
t as  0.5,  while  t-2.0  was  observed  by  Ksenzhek  et  al.  (75).  Simonsson 
(28)  and  Romanova  and  Selitsky  (76)  observed  an  effective  conductivity 
which  is  only  about  one  tenth  of  the  conductivity  of  the  free  electro- 
lytic solution.  They  attributed  this  great  reduction  to  the  gas  bubbles 
left  in  the  pores. 


Redox  Electrochemical  Rate  Expression 


The  rate  of  the  charge  transfer  reaction  A at  the  solid-solution 


interface  can  be  represented  by  the  following  equation. 


c*  * c*  / “Sf  (V  " ~5T  CV 


J - i°0  C 
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where  or  and  a are  the  effective  anodic  and  cathodic  transfer  coeffi- 

a c 

cients,  respectively;  j is  the  local  transfer  current  per  unit  area  of 
active  zinc  surface  which  acts  as  active  sites  for  reaction  (A) ; ig  is 
the  exchange  current  density  evaluated  at  reference  concentrations  c° 
and  Cjj  y and  C are  the  orders  6f  dependence  of. the  exchange  current 
density  1q  on  zincate  and  hydroxide  concentrations,  respectively;  6^  is 
the  potential  of  the  matrix  phase;  and  6*  is  the  potential  of  the 
reference  electrode  of  the  same  kind  as  the  working  electrode  located 
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just  outside  the  double  layer. 

The  conservation  of  charge  requires 
di, 

a j - (30) 

m dx 

where  a^  is  the  solid-solution  interface  area  per  unit  volume  of  elec- 
trode which  acts  as  active  sites  for  reaction  (A). 

If  the  diffusion  of  a species  to  or  from  the  reaction  sites  is 


assumed  to  proceed  through  a well  defined  characteristic  diffusion 
length,  6 , the  concentration  of  zincate  ion  (Cj)  and  hydroxide  ion  (c^) 
at  the  metal  surface  can  be  determined  from  the  following  equations. 

(31) 


(32) 


di. 

nF 

c 

■ - — 

5 

*VC2 

- c2) 

dx 

2 

di. 

nF 

o 

2 

a - 

S _ 

. k,(=3 

- c3) 

dx 

3 

where  a is  the  solid-solution  interface  area  per  unit  volume  of  the 
electrode,  is  the  mass  transfer  coefficient  (in  cm/sec)  of  species  i 
from  the  bulk  solution  in  the  pores  to  the  reaction  sites  which  are 
active  for  charge  transfer  reaction  (A)  or  vice  versa. 

The  potential  at  the  metal  surface,  <j>2»  is  assumed  to  be  the  same 
as  the  potential  of  the  bulk  solution  in  the  pores,  4>2*  Equations  (29) 
through  (32)  are  combined  to  give 


crF 

- “RT  CV  V 
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(33) 


The  constants  y,  C,  a , and  a .corresponding  to  the  reaction  mechanism 
(E)  to  (H)  of  Bockris  et  al.  (52),  are  0.75,  0.0,  1.5,  and  0.5,  respec- 
tively. The  exchange  current  density,  ip,  reported  by  a number  of 
investigators  (52,77,78,79,80)  has  the  values  in  the  range  of  10  to  300 
mA/cn2. 

Species  Conservation  Equation 

The  use  of  appropriate  averaging  quantities  yields  the  following 
conservation  equation. 

9cc.  d N. 

♦ R.  (34) 

at  dx  1 

where  represents  production  or  consumption  of  species  i due  to  the 
electrode  reaction.  During  discharge,  zincate  ion  is  produced  by  charge 
transfer  reaction  (A)  and  consumed  by  precipitation  reaction  (B) . On 
the  other  hand,  hydroxide  ion  is  consumed  by  reaction  (A)  and  produced 
by  reaction  (B) . The  reaction  source  R^  can  be  written  as 

"*■-  it ' Si*k»  <c*‘c^  l3S) 

♦ 

where  k , during  discharge,  is  defined  by  the  rate  constant  combining 

9 

mass  transfer  of  zincate  ions  from  bulk  solution  in  the  pores  to  the 
reaction  sites  with  chemical  rate  constant  for  precipitation  of  ZnO; 
c^  is  the  equilibrium  (or  saturation)  concentration  of  zincate  ion  at 

the  solid-solution  interface  which  are  active  for  reaction  (B) ; and  s^ 

* 

and  s^  axe  the  stoichiometric  coefficients  of  species  i for  the  reac- 
tion (A)  and  for  the  reaction  (B),  respectively,  defined  by  equation 
(10).  Substitution  of  equation  (35)  into  equation  (34)  yields 
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where  S2m-S2--1,  s^"^,  s^-2t  8q-0,  and  sq--1  have  been  substituted. 
Overall  Conservation  Equation 

The  porosity  of  the  electrode  will  increase  as  zinc  dissolves  and 
decrease  as  ZnO  precipitates.  A material  balance  on  the  solid  phases 
shows  how  the  porosity  changes  with  time  as  shown  below. 

17’  rf  + £ V <*i  - (39) 


solid 

species 


solid  * 
species 


where  is  the  partial  molar  volume  of  species  i. 

Similarly,  a material  balance  on  the  species  in  the  solution 
phase  gives  the  following  electrolyte  continuity  equator. 


9»  f dvD_  / V2  + VBt3  + + VBS3  + Vo \ ^2 

9 1 dx  ' z2V2  z3v3  >4  n v3  n ' F dx 


This  equation  is  obtained  by  adding  equations  (36),  (37),  and  (38) 

A B 

which  have  been  multiplied  by  Vg/ v“,  and  VQ,  respectively,  and 


I. 


using  the  relations  of  c/^V^c^-1 , V^V^-FO^O.  and 

VaDad*V0Ddd-V  D »0.  The  last  two  relations  for  Dn.  and  D__  are 
A AB  B BB  0 OB  OA  OB 

obtained  from  equations  (21) . 

From  equations  (39)  and  (40) , the  convective  bulk  flow  motion  in 
the  porous  electrode  can  be  described  by  the  following  overall  continuity 


equation. 


1 trv.  t^V.  di, 
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The  second  bracket  term  is  the  volume  change  of  reaction  (B) . The 
first  bracket  term  represents  the  volume  change  due  to  reaction  (A)  if 
equation  (24)  is  applied. 

Ohm's  Law  in  the  Solid  Matrix 


Ohm's  law  in  the  solid  matrix  can  be  written  as 


ij  . I - i2  . - <r  e. 


where  ^ is  the  potential  of  the  solid  matrix,  and  is  the  volume  of 
conducting  material  per  unit  volume  of  the  electrode  (not  necessarily 
equal  to  1-e).  The  superficial  current  density  in  the  solid  phase  ij 
and  that  in  the  solution  phase  i^  are  related  to  the  total  applied 
current  density  as  ij+i2«I.  The  conductivity  of  solid  matrix  is  approx- 


imated by 


3 ■ «L»cw**  * "i"/ 


where  a^,  a®n0»  and  °i  are  1116  conductivities  of  the  zinc  metal,  zinc 
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oxide,  and  inert  materials  in  the  pure  solid  state,  respectively,  and 

eZn*  eZnO’  £I  are  t*ie  zinc,  zinc  oxide,  and  inert 

materials,  respectively,  per  unit  volume  of  the  electrode. 

Mass  Transfer  Coefficients  (k^,  kg,  and  k*) 

The  reaction  rate  is  determined  by  the  rate  of  mass  transfer  of  a 
limiting  species  between  bulk  solution  within  the  pores  and  the  active 
reaction  sites.  During  discharge,  hydroxide  ion  is  depleted  at  the 
metal  surface,  producing  zincate  ion.  The  zincate  ions  are  expected  to 
diffuse  into  the  bulk  solution  in  the  pores  resulting  in  a supersatu- 
rated zincate  solution.  At  the  same  time,  the  zincate  ions  in  the 
supersaturated  solution  will  diffuse  to  the  ZnO  surface  to  precipitate 
as  ZnO.  The  limiting  species  at  the  metal  surface  during  charge  is  the 
zincate  ions  which  must  be  supplied  from  the  bulk  in  the  pores.  The 
solution  becomes  undersaturated  with  ZnO  and  solid  ZnO  will  dissolve 
into  the  solution. 


Based  on  the  assumption  that  there  exists  a characteristic  dif- 
fusion length,  6 , which  is  the  distance  between  the  solid-solution 
interface  and  the  bulk  solution  inside  the  pores,  mass  transfer  coeffi- 


* 

cients  k^,  kg,  and  ks  defined  in  equations  (30),  (31),  and  (35)  are 
described  by  the  following  equations  (see  Appendix  B for  derivation) 


1 1 1 


where  and  kg  are  the  mass  transfer  coefficients  of  potassium  zincate 
and  potassium  hydroxide,  respectively,  transferring  from  metal  surface 

to  the  bulk  within  the  pores  or  vice  versa,  6 is  the  characteristic 

* 

diffusion  length,  kg  is  the  rate  constant  combining  mass  transfer  of 
zincate  from  bulk  to  ZnO  surface  with  chemical  rate  constant  for  preci- 
pitation of  ZnO,  and  and  are  the  diffusion  coefficients  of  the 
potassium  zincate  and  potassium  hydroxide,  respectively.  The  area  a^ 
and  ag  represent  the  interfacial  areas  which  are  active  for  electro- 
chemical reaction  (A)  and  for  chemical  reaction  (B) , respectively. 

Surface  Area  (a  , a , and  a) 

■ m s 

The  value  of  the  true  electrochemically  active  area  is  difficult 
to  determine.  To  a first  approximation,  it  is  assumed  that  the  macro- 
scopic averaged  particles  with  uniform  radius  are  arranged  in  face 
centered  cubic  closed  packed  positions.  For  a given  reference  porosity 
e°  and  radius  R°,  the  solid-solution  interface  area  a°  per  unit  volume 
of  the  electrode  can  be  estimated  from  the  relation  of  a°«3(l-e°)/R0. 
The  specific  surface  area  of  the  electrode  at  any  state  having  porosity 
of  e can  be  related  to  the  reference  interface  area  a°  as  a/a°»(l-e)R°/ 
(l-e9)/R.  Since  the  ratio  R/R°  is  equal  to  (a/a°)0'5  , the  specific 
area  a can  be  represented  as  follows: 
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The  specific  interface  area  a is  the  sum  of  the  area  a which  is  active 

m 

for  charge  transfer  reaction  (A)  and  the  area  ag  which  is  the  active 

site  for  dissolution  or  precipitation  of  ZnO. 

In  the  absence  of  the  inert  conducting  material,  the  area  a and  a 

ms 

become  a^n  (zinc  area)  and  aZnQ  (ZnO  area),  respectively.  They  are 
approximated  by  the  following  equation. 


% tor  V * * - zZnO 


<eZn>' 


<'zn>P*  "zn</ 


(46) 


where  z^q  is  the  volume  fraction  of  ZnO,  £Zn  is  the  volume  fraction  of 
active  zinc  which  can  be  converted  into  ZnO,  and  p and  q are  the  cons- 
tants relating  volume  to  surface  area  which  depend  on  morphology  of  the 
2 

crystals  (p«q-  y for  cubical  or  spherical  crystals) . The  volume  frac- 
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tion  of  active  zinc,  EZn>  is  defined  by 
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plus  covx 

where  * and  eZn  represent  the  volume  fractions  of  nonactive  zinc 

caused  by  pore  plugging  (blockage  of  the  1st  kind)  and  complete  coverage 

on  zinc  surface  by  ZnO  (usually  called  passivation  or  blockage  of  the 

2nd  kind),  respectively. 

covr 

The  quantity  EZn  is  approximated  by 
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The  meaning  of  the  various  volume  fractions  is  explained  by  using  the 
following  example. 

Let  us  consider  a local  section  of  the  electrode  having  porosity 
of  0.3,  60%  Zn  by  volume,  and  10%  ZnO  by  volume.  The  partial  molar 

3 

volumes  of  zinc  and  zinc  oxide  are  9.15  and  14.51  cm  /mole,  respectively. 
Due  to  the  differences  in  the  partial  molar  volumes,  85.3%  conversion 
of  Zn  into  ZnO  results  in  pore  plugging,  i.e.,  e»0,  eZn*0.088,  and  eZn0* 
0.912.  The  remaining  Zn  can  not  be  converted  into  ZnO,  and  is  denoted 
by  e|^u®(«0.088) . However,  complete  surface  coverage  of  zinc  by  ZnO 
might  occur  before  pore  plugging.  If  complete  conversion  of  ZnO  into 
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tion.  The  six  equations  can  be  solved  by  use  of  numerical  techniques 
with  the  following  initial  and  boundary  conditions. 


Initial  conditions  at  time»0 
o 

c2  * c2 
o 

C3  " C3 


For  the  other  variables,  the  analytical  solutions  of  the  following 
approximated  equations  are  used  to  estimate  the  initial  profiles  of  i2» 
d>1 , and  $2* 

d<b 

i2  - - <°(e°)Ut 
L dx 

di,  (a  ♦ a )F 

i°  C—  (♦.  - <J>,)  (50) 

dx  " RT 

„ n d *i 

ij  - I - i2  - - a eB  

dx 

where  the  superscript  o refers  to  the  values  evaluated  at  time«0,  and 
I is  the  applied  current  density.  The  solution  of  equation  (50)  can  be 
written  as  follows  (e.g.,  for  current  distribution) 
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where  y is  the  dimensionless  distance  from  the  backing  plate  (x/L),  L is 

the  electrode  thickness,  a is  the  effective  matrix  conductivity  (o° e°) , 
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and  k is  the  effective  solution  conductivity  (tc°(e°) *+t) . 
6 

The  boundary  conditions  at  x*0  are  as  follows: 
d<jL 

— - - 0 (or  i_*0) 
dx 


dx 

dc,  (52) 

— i . 0 

dx 

\P  - 0 

<i>j  * o 

At  the  boundary  x-L,  the  current  density  in  the  solution  phase  (i2) 
is  the  same  as  the  applied  current  density  (I).  For  the  species  bound- 
ary conditions,  the  limited  mass  transfer  of  the  various  species  across 
the  boundary  layer  should  be  accounted  for  to  give  a better  approxima- 
tion to  an  actual  battery  system.  Most  of  the  aqueous  secondary 
batteries  have  a separator  adjacent  to  the  electrode  surface.  Membranes 
have  been  used  as  separators.  In  some  cases,  porous,  inert  materials 
(for  example,  dynel  cloth)  which  contain  excess  solution  reservoir  are 
placed  between  the  electrode  surface  and  the  membrane. 

If  the  solution  reservoir  is  assumed  to  be  completely  mixed,  the 
changes  in  amount  of  species  i with  time  in  the  solution  reservoir  are 
equal  to  the  flux  of  species  i across  the  boundary  at  x*L,  N^,  minus 
the  flux  of  species  i across  the  boundary  layer  (or  membrane) , N^1 . 

Based  on  this  concept,  the  boundary  conditions  at  x«L  can  be  written 
as  follows: 

At  x»L  i2  * I 


1 
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boundary  conditions  simply  represent  the  conservation  of  fluxes  at  the 
boundary  , i.e.,  Nb-  Nbl  and  Nj-  Nb*. 

The  fluxes  Nb  and  are  given  by  equations  (25)  and  (26) , respect- 
ively, modified  to  describe  the  fluxes  at  the  boundary  at  x-L.  The 
fluxes  Nb^  and  are  written  as  follows: 


_ bl  bl 
1 + CA  v 


Dbl  tbl 

d Z2V? 

Dbl  tbl 

Cc  . c«  vbI 

d B “ 2jvBF  B 


(54) 


where  d is  the  thickness  of  the  boundary  layer  (or  membrane) . 

The  six  equations  can  be  solved  by  numerical  techniques  developed 
by  Newman  (81,82)  subject  to  the  initial  conditions  (49)  and  (SO)  and 
the  boundary  conditions  (52)  through  (54) . The  electrode  length  of  L 
was  divided  into  NJ  mesh  points  separated  by  a distance  h between  mesh 
points.  The  boundaries  were  at  J«1  (x«0)  and  J-NJ  (x-L) . Mesh  point  J 
refers  to  a distance  of  (J-l)h  from  the  backing  plate.  In  the  computer 
programming,  the  variables  were  denoted  as  follows: 


C(1,J)  - i2(J) 
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C(2,j)  - <i»2(J) 

C(3,J)  - c2(J)  or  cA(J) 

C(4,J)  - c3(J)  or  cB(J) 

C(5,J)  - vRj) 

C(6,J)  - ^(J) 

In  order  to  apply  Newman's  technique,  the  six  equations  were 

linearized  about  a trial  solution  and  put  into  finite  difference  form. 

Central  difference  form  was  used  for  equations  (36)  and  (37) . Forward 

difference  form  for  equation  (28)  and  backward  difference  form  for 

equations  (33),  (41),  and  (42)  were  employed  with  the  coefficients  of 

terms  being  averaged  between  two  appropriate  mesh  points.  The  result* 

2 

ing  six  finite  difference  equations  are  correct  to  second  order,  h . 

At  the  boundaries  J«1  and  J>NJ,  the  following  finite  difference 
forms  were  used  for  the  first  order  boundary  conditions. 


For  J-l, 

a(J)  » a(J) 

dx 

For  J-NJ 
dC(K,J) 


-3C(K,J)+4C(K,J+l)-C(K,J+2) 


♦ 0(h2) 


C(K,J-2) -4C(K,J-1)+3C(K,J) 


♦ 0(h2) 


where  K refers  to  the  unknown  variable  and  J refers  to  the  mesh  point. 

The  Crank-Nicolson  implicit  method  (83)  was  used  for  averaging 
time  derivatives  in  equations  (36)  and  (37) . For  example,  considering 
the  following  equation, 
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the  finite  difference  form  used  is 

C(3,J)-Z(3,J)  1 C(3,J+1)+C(3, J-l) -2C(3,J) 

- ~ a(J)  

At  2 2 h 

1 Z(3,J+1)+Z(3,J-1)-2Z(3,J) 

♦ ~ zq(J)  

2 2 h 

where  Z and  za  represent  the  values  of  C and  a of  the  previous  time 
step,  respectively. 

The  resulting  matrix  was  solved  using  numerical  techniques 
developed  by  Newman  (81,82).  The  convergence  criterion  used  was  10' 
which  is  defined  by  the  ratio  of  the  difference  between  the  two  suc- 
cessive solutions  to  the  present  solution. 

The  numerical  calculations  require  a knowledge  of  the  various 
parameters  used  in  the  mathematical  model.  The  conductivity  of  the 
electrolytic  solution  (84),  solubility  of  ZnO  in  KOH  solution  (85),  and 
activity  coefficients  of  KOH  and  potassium  zincate  (68,86)  were 
expressed  as  a function  of  concentration  based  on  the  data  taken  from 
various  literature  reports.  Those  formulas  and  references  are 
described  in  Appendix  C.  Other  chemical  and  physical  parameters  perti- 
nent to  the  system  under  consideration  are  described  in  Appendix  D-l. 
The  computer  programs  and  the  control  parameters  are  shown  in  Appendix 
D.  The  numerical  solutions  were  carried  out  in  double  precision  on  the 
UCLA  IBM  360/91  digital  computer. 


III.  ZINC  ELECTRODE  BEHAVIOR  PREDICTED  FROM 
THE  MATHEMATICAL  MODEL 


The  mathematical  model  described  in  the  previous  section  was  solved 
to  predict  the  behavior  and  failure  modes  of  the  porous  zinc  electrode 
having  the  following  initial  conditions.  The  electrolytic  solution 
initially  contained  in  the  pores  is  taken  as  8M  KOH  solution  (about  35 
wt%)  saturated  with  zincate  (c°»0.66M).  The  electrode  under  considera- 
tion is  initially  40%  void  volume,  50%  zinc,  and  10%  ZnO  by  volume  (e°“ 
°*4’  ESn-°-5’  ^O"0-1’  and  e j*0) . The  zinc  electrode  behavior  during 

galvanostatic  operation  with  superficial  current  densities  of  50  and 
2 

20  mA/cm  were  considered.  These  current  densities  are  believed  to  be 
high  enough  to  represent  those  required  during  operation  of  common 
secondary  batteries  and  low  enough  to  eliminate  anodic  passivation,  . 
at  least  in  the  initial  stages.  Therefore,  it  is  assumed  that  complete 
coverage  of  the  zinc  surface  by  precipitated  ZnO  (passivation)  does  not 

COVT 

occur;  e • is  zero, 
in 

A general  membrane  boundary  case,  where  a membrane  is  placed 
between  the  zinc  electrode  surface  and  the  free  solution,  is  considered 
by  taking  the  parameters  shown  in  Table  1 (type  I).  However,  some  of 
the  results  for  the  other  boundaries,  when  they  appear  to  be  useful, 
are  also  included;  for  example,  solution  boundary  (the  third  column) 
and  cation  exchange  RAI  P2291  membrane  boundary  (the  fourth  column  of 
Table  I) . The  other  parameters  used  in  this  calculation  are  shown  in 
Appendix  D-l.  The  predicted  results  during  galvanostatic  discharge  of 
a porous  zinc  electrode  are  shown  in  Figures  2 to  7.  Figures  8 to  13 
illustrate  the  sunnarized  zinc  electrode  behavior  during  cycling. 
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Table  I 


Transport  Parameters  Applied  to  the 
Various  Boundary  Layers 


Type  1 

Membrane 

Boundary 

Solution 

Boundary 

Type  II 
Membrane 
(RAI  P2291) 

'S1 

0.5 

0.23 

0.5 

0.0 

0.05 

0.0 

0.5 

0.72 

0.5 

„bl  , 2,  . 

daa  (cm  /8) 

2 . 5xl0~8 

6 . 0xl0"6 

l.OxlO-8 

_bl  , 2.  . 

dbb  (c$  /8) 

2.5x10"' 

2 . OxlO"5 

7.5xl0-8 

d (cm) 

0.0025 

0.005 

0.0025 

laa 

l.OxlO"5 

0.0012 

4. OxlO-6 

*5  <ca/8> 

l.OxlO-4 

0.004 

3. OxlO-5 

* Transport  parameters  for  RAI  P2291  membrane  were  taken  from 
the  data  of  Slnha  and  Bennlon. (87) . The  parameter  L^}  Is 
defined  by  D^l/d  (refer  to  equation  54).  ~ 


III-l.  Discharge  Behavior 


r 


Figure  2 shows  the  anodic  overpotentials,  <Jl  (o) - ^ (L) , defined  by 

2 

the  potential  of  the  backing  plate  minus  the  potential  in  the  solution 
at  the  electrode  surface,  as  a function  of  time  or  depth  of  discharge. 
The  shape  of  potential -time  curve  resembles  those  reported  by  Elsdale 
et  al.  (88).  In  contrast  to  the  sharp  increase  at  a point  of  passiva- 
tion on  a planar  zinc  electrode  (for  example,  see  references  14  and  15), 
the  potential-time  curve  Is  characterized  by  an  Initial  gradual  rise  in 

overpotential  followed  by  a steeper  rise  at  a later  stage.  The  overpo- 
tential goes  up  substantially  in  the  region  of  40%  depth  of  discharge. 

It  is  to  be  noticed  that  even  in  the  idealized  model,  only  40%  of  ini- 
tial zinc  metal  can  be  converted  into  ZnO  prior  to  electrode  failure. 

The  electrode  failure  at  40%  depth  of  discharge  is  caused  by  high  ohmic 
loss  in  the  solution  due  to  the  decrease  in  pore  size  as  well  as  KOH 
concentration  and  partly  by  high  activation  and  concentration  overpoten- 
tials. 

Figure  3 shows  how  the  concentration  profiles  for  KOH  (top  curves) 

and  potassium  zincate  (bottom  curves)  develop  during  discharge  at  50 
2 

mA/cm  current  density.  Initial  concentration  of  the  solution  was  8M 
KOH  saturated  with  zincate  (0.66M).as  shown  by  dotted  lines.  After  one 
minute  polarization,  KOH  concentration  decreases  and  zincate  becomes 
supersaturated  (curve  A) . On  further  discharge,  KOH  concentration 
decreases  further  to  give  curves  B for  10%,  C for  20%,  and  D for  30% 
depth  of  discharge.  Saturation  concentration  of  zincate  decreases  as 
KOH  concentration  decreases.  Therefore,  zincate  solution  concentration 
also  decreases  as  discharge  continues;  however,  there  is  always  some 
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degree  of  supersaturation.  The  significant  depletion  of  KOH  electrolyte, 
e.g.,  4.5M  at  30%  depth  of  discharge  is  due  to  the  restricted  flux  of 
hydroxide  ions  across  the  membrane.  This  phenomenon  become  more  obvious 
in  Figure  4 where  KOH  concentration  profiles  for  three  types  of  bound- 
aries are  presented.  The  average  KOH  concentration  after  10%  depth  of 
discharge  becomes  7.5M  for  solution  boundary  (curve  A)  but  becomes  2.0M 
for  type  II  membrane  boundary  (curve  C) . Clearly,  the  cell  having  type 

II  membrane  as  a separator  can  fail  due  to  KOH  depletion. 

2 

The  charge  transfer  current  density,  j (mA/cm  ) , is  plotted  as  a 
function  of  distance  in  Figure  5.  Initially,  the  distribution  of  charge 
transfer  reaction  is  shown  as  curve  A,  having  the  maximum  at  the  front 
face  of  the  electrode.  As  discharge  continues  up  to  10%  depth  of  dis- 
charge, KOH  concentration,  porosity,  and  therefore  conductivity  of  solu- 
tion decrease.  Furthermore,  the  increase  in  the  surface  zincate 
concentration  at  the  zinc  surface  results  in  smaller  activation  cverpo- 
tential.  These  two  effects  shift  the  reaction  toward  the  front  face 
and  the  reaction  profiles  become  less  uniform  as  shown  in  curve  B.  If 
the  electrode  is  discharged  further  beyond  the  20%  depth  of  discharge, 
the  reaction  at  the  face  of  the  electrode  starts  decreasing  due  to  the 
depleted  active  zinc  metal  (curve  C).  For  40%  discharge,  all  the  zinc 
near  the  front  face  of  the  electrode  is  used  up  as  shown  in  curve  0. 

The  reaction  profile  is  highly  nonuniform  and  the  effective  reaction 
zone  moves  step  by  step  from  the  face  of  the  electrode  toward  the  back- 
ing plate.  Similar  profiles  are  obtained  for  the  other  two  boundaries. 

The  distribution  of  solid  zinc  and  zinc  oxide,  in  terms  of  volume 
fractions,  is  plotted  in  Figure  6 for  type  I membrane  and  in  Figure  7 


CONCENTRATION  OF  KOH  (MOLES /LITER) 


DISTANCE  FROM  BACKING  PLATE  (cm) 


Figure  4.  Concentration  Profiles  of  KOH  After  10Z  Depth 
of  Discharge  (1*50  mA/cm^  and  A,B,and  C for 
Solution  Boundary,  Type  I Membrane  Boundary, 
and  Type  II  Membrane  Boundary,  Respectively) 
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TYPE  I MEMBRANE  BOUNDARY 


DISTANCE  FROM  BACKING  PLATE  (cm) 


Figure  5.  Distribution  of  Local  Charge  Transfer  Current, j(mA 
/a*2),  for  Type  I Membrane  Boundary  (1-50  wA/cvfi, 
L-0.1  cm,  and  A,B,C,and  D Represent  0Z,10Z,20Z, 
and  40Z  Depth  of  Discharge,  Respectively) 


DISTANCE  FROM  BACKING  PLATE  (cm) 


Distribution  of  Zn  and  ZnO  Plotted  As  Volume 
Fractions  (cm3 /cm’)  for  Type  I Membrane  Boun- 
dary (1-50  mA/cm2,  and  AfB(C(and  D Represent 
10Z,20Z,30Z,and  40Z  Depth  of  Discharge,  Res- 
pectively) 


VOLUME  FRACTIONS  «„  AND  „ (cmJ/cm 

IXi  ZnO 


Figure  7. 


DISTANCE  FROM  BACKING  PLATE  (cm) 


Distribution  of  Zn  and  ZnO  Plotted  ea  Volume  Fractions 
(cm^/car ) for  Solution  Boundary  (1*50  mA/cm%  and  A,B, 
C,and  D Represent  10Z,20Z,30Z,and  40Z  Depth  of  Discharge, 
Respectively) 
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for  solution  boundary.  As  shown  in  Figure  5,  solid  zinc  is  consumed 
mainly  near  the  face  of  the  electrode  and  used  up  after  only  23%  dis- 
charge for  the  membrane  boundary  (Figure  6)  and  10%  discharge  foT  the 
solution  boundary  (Figure  7).  During  discharge,  solid  zinc  dissolves 
into  solution  by  charge  transfer  reaction  (A) , producing  supersaturated 
zincate  solution  from  which  ZnO  is  precipitated  by  chemical  reaction 
step  (B).  Consequently,  the  total  conversion  of  zinc  into  zincate  ion 
at  a specified  time  t,  i.e.,  the  integrated  area  between  the  dotted  line 
(denoted  by  and  the  curve  at  t*me  01151  *he  saoe  f°r  both 
cases.  However,  the  total  amount  of  ZnO  precipitated,  i.e.,  the  inte- 
grated area  of  ZnO  profile,  is  smaller  for  the  solution  case  than  for 
the  membrane  case.  For  example,  after  40%  depth  of  discharge,  the 
amount  of  precipitated  ZnO  obtained  from  Figure  6 and  Figure  7 are  95% 
and  55%,  respectively,  of  the  total  amount  of  the  consumed  zinc  which 
is  calculated  from  the  current  applied  during  the  operation  period  (8400 
seconds).  The  missing  zinc  oxide  is  lost  into  the  counter  electrode 
compartment  as  a soluble  species,  zincate  ion,  before  precipitation 
occurs. 

III-2.  BehavioT  on  Cycling 

The  electrode  initially  composed  of  40%  liquid  filled  pores,  50% 

Zn,  and  10%  ZnO  by  volume  was  cycled  with  constant  current  density  of 
2 

SO  mA/ca  . For  each  cycle,  20%  depth  of  discharge  followed  by  20%  depth 
of  charge  are  applied.  At  the  beginning  of  each  discharge  or  charge 
step,  the  concentrations  of  zincate  and  hydroxide  ions  are  assumed  to  be 
c°(0.66M)  and  c^(8M),  respectively.  This  assumption  is  not  significant 


since  the  time  required  to  reach  or  c^  is  negligible  compared  to  the 
total  operation  time  at  each  half  cycle.  For  example,  at  the  end  of  20% 
discharge  of  the  first  cycle,  the  solution  is  5M  KOH  for  type  I membrane 
boundary  case  and  7.4M  KOH  for  the  solution  boundary  case.  The  time 
required  after  cycle  reversal  for  the  solution  to  reach  8M  KOH  is  about 
3 minutes  for  type  I membrane  boundary  and  30  seconds  for  the  solution 
boundary  case. 

Figure  8 shows  the  concentration  profiles  of  hydroxide  (top  curves) 
and  zincate  (bottom  curves)  during  charge  of  the  first  cycle.  In  Figure 
9,  transfer  current  distribution  on  charging  is  compared  to  that  on  dis- 
charging for  type  I membrane  boundary  case.  Current  distribution  during 
charge  becomes  more  uniform  since  precipitated  ZnO  initially  reduces 
the  active  zinc  surface  area  and  shifts  the  reaction  towards  the  backing 
plate.  The  overpotential  is  plotted  as  a function  of  time  in  Figure  10. 
The  anodic  overpotentials  increase  to  give  20  mV,  24mV,  and  32mV  at  the 
end  of  discharge  of  the  1st,  2nd,  and  3rd  cycle,  respectively.  Similar 
but  less  significant  increase  in  overpotential  is  observed  during 
charge.  This  increase  in  overpotential,  namely,  the  decrease  in  cell 
capacity  with  cycling  is  due  to  the  redistribution  of  zinc  and  zinc 
oxide  which  causes  pore  plugging  near  the  front  face  of  the  electrode. 

The  volume  fractions  of  zinc  and  zinc  oxide  during  two  cycles  are 
plotted  as  a function  of  distance  in  Figure  11  for  the  type  I membrane 
boundary  and  Figure  12  for  the  solution  boundary.  The  distribution  of 
zinc  and  zinc  oxide  becomes  less  uniform  as  cycling  continues.  This 
nonuniformity  is  more  serious  for  the  solution  boundary  (Figure  12)  than 
for  the  membrane  boundary  (Figure  11).  The  volume  fractions  of  Zn  and 


DISTANCE  FROM  BACKING  PLATE  (cm) 


Figure  8.  Concentration  Profile  During  Constant  Current 
Charge  for  Type  I Membrane  Boundary  (I“50  mA/ 
cm2,  l-0.1  cm,  and  A for  One  Minute  Polariza- 
tion and  B and  C for  10Z  and  20Z  Depth  of 
Charge,  Respectively) 


VOLUME  FRACTIONS  eZn  AND  eZnQ 


FIRST  CYCLE 


CHARGE 


DISCHARGE 


SECOND  CYCLE 


DISTANCE  FROM  BACKING  PLATE  (cm2 


Distribution  of  Zn  and  ZnO  Plotted  As  Volusia 
Fractions  (cm3/cm3)  During  Two  Cycles  for 
Type  I Membrane  Boundary  (I«50  mA/cn2,L-0.1 
ca,«2  -O.S.a®  .-0.1) 


FIRST  CYCLE 


MEMBRANE  BOUNDARY  (1-50  tnA/cnV*) 


MEMBRANE  BOUNDARY  (1-20  mA/cm2) 


DISTANCE  FROM  BACKING  PLATE  (cm) 


Distribution  of  Zn  and  ZnO  Plotted  As  Volume 
Fractions  At  Two  Current  Densities  for  Type 
I Membrane  Boundary.  —.After  One  Cycle; 
After  Two  Cycles. 
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ZnO  for  the  low  current  density  of  20  mA/cm  are  compared  with  those  for 
the  high  current  density  of  SO  mA/cm  in  Figure  13.  For  each  cycle,  20% 
depth  of  discharge  followed  by  20%  depth  of  charge  was  applied.  Pro-' 
files  of  Zn  and  ZnO  are  clearly  more  uniform  for  the  low  current  density 
compared  to  those  for  the  high  current  density. 

III-3.  Zinc  Electrode  Failure  Modes 

Zinc  electrode  behavior  occurring  during  cycling  has  been  predicted 

by  a mathematical  model.  The  model  predicts  that  discharge  of  the  zinc 

2 

electrode  having  porosity  of  0.4  with  current  density  of  SO  mA/cm 
failed  after  only  40%  depth  of  discharge  (see  Figure  2)  due  to  pore 
plugging  (blockage  of  the  1st  kind)  as  shown  in  Figure  6.  The  model 
also  predicts  that  the  reaction  distribution  is  highly  nonuniform,  and 
the  effective  reaction  zone  is  very  thin  (Figure  5) . This  nonuniform 
reaction  distribution  will  accelerate  the  failure  due  to  pore  plugging 
at  the  face  of  the  electrode.  Furthermore,  the  zinc  metal  in  the  region 
outside  the  effective  reaction  zone  serves  simply  as  an  inert  matrix  and 
does  not  contribute  to  the  cell  capacity.  Therefore,  thin  electrodes 
are  required  for  design  of  industrial  cells  in  order  to  have  high  cell 
capacity  per  unit  weight  of  electrode. 

The  loss  of  ZnO  species  into  the  counter  electrode  compartment 
during  discharge  is  very  large  for  cells  without  a membrane  (see  Figure 
7),  while  for  the  cell  with  a membrane  (Figure  6),  mass  transfer  of 
zincate  ion  into  the  counter  electrode  compartment  was  restricted. 
During  charge,  for  the  cell  without  a membrane  or  with  a punctured 
membrane,  the  zincate  ions  which  have  moved  towards  the  counter  elec- 
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trode  on  discharging  are  expected  to  move  back  towards  the  zinc  test 
electrode.  This  type  of  mass  transport  leads  to  zinc  dendrite  formation 
which  eventually  causes  short  circuits  to  the  counter  electrode.  This 
is,  probably,  the  main  reason  along  with  preventing  zinc  metal  nuclea- 
tion  within  the  membrane  why  industrial  cells  having  appropriate  mem- 
branes successfully  reduce  dendrite  formation  and  short  circuits. 

However,  for  cells  having  a membrane  such  as  RAI  P2291  cation 
exchange  membrane,  KOH  concentration  was  depleted  significantly  during 
discharge,  which  can  be  a cause  of  electrode  failure.  Discharge  failure 
due  to  KOH  depletion  was  predicted  in  Figure  4 (curve  C) . 

Shape  change,  i.e.,  the  redistribution  of  solid  zinc  species  in  the 

. 

y-direction,  has  been  successfully  explained  by  the  differences  in  con- 
vective  flow  and  the  degree  of  super-  (or  under-)  saturation  of  zincate 

j 

ion  during  cycling  (12,13).  Similarly,  redistribution  of  solid  zinc 
species  in  the  x-direction  is  predicted  as  shown  in  Figures  11  to  13. 

ZnO  precipitated  during  discharge  was  not  completely  dissolved  during 
charge  but  accumulated  at  the  electrode  surface  as  cycling  continues. 

This  movement  may  be  a limiting  factor  in  the  loss  of  cell  capacity 
after  many  cycles. 

Based  on  the  zinc  electrode  behavior  predicted  by  the  present 
mathematical  model,  failure  modes  occurring  in  the  zinc  electrode  are 
summarized  as  follows; 

Failure  during  discharge 

* 

(1)  Blockage  of  the  1st  kind  (pore  plugging) 

* 

(2)  Blockage  of  the  2nd  kind  (passivation  ot  surface  coverage  on  zinc 
surface  by  the  precipitated  ZnO) 


55 


* 

(3)  Depletion  of  KOH  solution 

(4)  Gas  evolution 

Failure  during  charge 

* 

(1)  Depletion  of  zincate 

(2)  Zinc  dendrite  formation 

(3)  Gas  evolution 

Failure  during  cycling 

(1)  Zinc  dendrite  formation  and  short  circuit 

(2)  Shape  change  in  the  y-direction 
# 

(3)  Redistribution  of  Zn  species  in  the  x-direction 

(4)  Accumulated  gas  evolution 

where  symbol  denotes  the  failure  modes  to  be  considered  in  the 
present  work. 
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IV.  EXPERIMENTS 


IV-1.  Preparation  of  the  Pressed  Zinc  Electrode 

The  test  electrodes  were  prepared  from  pure  zinc  powder  (Plasmadyne 

Co.)  with  particle  sizes  between  250  and  325  mesh  (60-44  microns).  The 

. 

zinc  powder  was  amalgamated  with  mercury  in  a mercuric  acetate  solution 
containing  an  amount  of  mercury  equal  to  2.5  wt%  of  the  zinc  weight. 

The  amalgamation  solution  also  contains  one  percent  glacial  acetic  acid. 
Complete  amalgamation  was  confirmed  by  adding  a mixture  of  NH4C1  and 
NH^OH  to  a small  sample  of  the  solution,  where  any  dissolved  mercury 
would  form  a white  precipitate.  After  amalgamation  is  complete,  the 
zinc  powder  is  rinsed,  dried  in  a vacuum  oven,  and  stored  in  a desicator. 

The  test  zinc  electrodes  were  prepared  by  pressing  amalgamated  zinc 
powder  onto  a backing  plate.  The  backing  plate  consisted  of  a titanium 
foil  (0.001  inches  thick)  upon  which  a silver  screen  (0.004  inches  thick) 
had  been  spot  welded.  Electrodes  with  two  different  porosities  were 
tested;  low  porosity  and  high  porosity  electrodes.  The  low  porosity 
electrodes  were  obtained  by  pressing  the  Zn(Hg)  powder  in  a stainless 
steel  mold  using  13000  psi.  The  resultant  electrode  contained  about  0.5 
grams  of  zinc  amalgam  and  had  a porosity  of  about  0.3  to  0.35.  The  high 
porosity  electrodes  were  prepared  by  the  procedure  given  by  Morrell  and 
Smith  (34).  A mixture  of. NaCl  powder  and  Zn(Hg)  powder  having  the 
weight  ratio  of  40  NaCl-60  Zn(Hg)  was  pressed  onto  a backing  plate  using 
20000  psi.  The  NaCl  powder  used  as  a filler  had  a particle  size  of  250 
to  325  mesh.  After  pressing,  the  resultant  electrode  was  soaked  in 
distilled  water  for  24  hours.  Complete  dissolution  of  NaCl  was  con- 
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finned  during  this  period.  This  gave  a highly  porous,  stable  structure 
having  a porosity  of  about  0.6.  After  the  electrode  was  dried,  a white 
precipitate  formed,  which  may  be  zinc-oxychloride,  zinc-hydroxy  chloride, 
or  zinc-hydroxy  carbonate.  This  compound  was  easily  dissolved  by  dip- 
ping it  in  dilute  HC1  solution.  The  resultant  electrode  contained  about 
0.3  g of  zinc. 

The  electrodes  made  by  the  above  procedures  were  disks  having  a 

2 

cross  sectional  area  of  1 cm  and  thickness  of  about  0.1  to  0.11  cm. 

2 

Small  cross  sectional  area,  1 cm  , was  chosen  in  order  to  minimize  the 
effect  of  nonuniform  current  distribution  in  the  direction  parallel  to 
the  electrode  surface. 

IV-2.  Potential  Measurements 

A schematic  drawing  of  the  test  cell  is  shown  in  Figure  14.  The 
main  body  of  the  cell  was  made  of  Plexiglas.  The  test  zinc  electrode 
was  completely  sealed  by  two  0-rings  when  the  cell  is  assembled  by 
tightening  screws.  A silver  screen  having  a zinc  deposit  was  used  as  a 
counter  electrode.  Two  layers  of  RAI  P2291  (40/20)  cation  exchange 
membrane  (thickness  of  0.001  inches)  and  one  layer  of  dynel  (thickness 
of  0.005  inches),  used  as  separators,  were  placed  on  top  of  the  test 
electrode.  The  RAI  membrane  is  a beta  radiation  cross  linked  low  den- 
sity polyethylene  film  which  has  been  gamma-radiation  grafted  with 
methacrylic  acid.  The  dynel  is  a porous  non-woven  acrylonitrile  which 
is  resistant  to  oxidation  in  the  cell  environment.  The  sealed  test  zinc 
electrode  contained  limited  electrolytic  solution  (about  0.03  to  0.11 
cm"*),  while  the  counter  electrode  compartment  contained  effectively  a 
large  amount  of  solution  (about  1.5  cm^). 
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Figure  1*.  Test  Cell  Assembly 


The  teflon  capillaries  lead  to  the  Hg/HgO  reference  electrodes. 

The  Hg/HgO  reference  electrodes  were  made  by  the  procedure  given  by  Falk 
and  Salkind  (89) . The  H-shaped  glass  tube  contained  40  wt%  KOH  solution 
which  contacted  with  a mercury  oxide -mercury-  platinum  wire  couple.  A 
dynel  fiber  plug  was  inserted  into  the  end  of  the  capillary  tube  to 
minimize  the  flow  of  electrolytic  solution  from  the  zinc  electrode  com- 
partment into  the  reference  electrode  or  vice  versa. 

After  evacuation,  filling  the  cell  with  solution,  and  assembling 
the  cell,  the  test  electrodes  were  discharged  galvanostatically  in  a 
40  wt%  KOH  solution  saturated  with  ZnO.  Figure  15  is  the  basic  circuit 
used  in  operation  of  the  cell.  The  overpotentials  of  the  test  zinc 
electrode  were  measured  with  respect  to  Hg/HgO  reference  electrodes 
using  a Keithley  Model  602  electrometer  and  one  channel  of  a Varian 
Associate  G2000  two-channel  recorder.  Three  Hg/HgO  reference  elec- 
trodes were  used  to  measure  the  zinc  electrode  potentials  at  three  dif- 
ferent positions;  one  above  the  membrane,  one  on  the  electrode  surface, 
and  one  on  the  backing  plate.  An  automatic  multipole  switch  and  a 
Flexopulse  HG  108A6  timer  were  used  to  allow  periodic  measurements  at 
three  positions  with  intervals  of  20  seconds.  An  external  zinc -mercury 
cell  was  connected  in  series  in  the  electrometer  circuit  to  cancel  out 
the  equilibrium  potential  of  the  zinc  electrode  with  respect  to  Hg/HgO 
reference  electrode.  This  arrangement  allowed  zero  potential  at  open 
circuit  and  the  sensitivity  of  the  potential  data  was  improved.  Poten- 
tial differences  between  the  zinc  test  and  the  zinc  counter  electrodes 
were  recorded  continuously  using  the  other  channel  of  the  recorder. 

The  integrated  current  passed  was  counted  by  a coulometer,  Model  541 


Figure  15.  Electrical  Circuit  for  Operation  of  Cell. 

A,  ammeter;  C,  counter  electrode;  El,  electrometer  for 
monitoring  teat  to  counter  electrode  potential  difference; 

E 2,  electrometer  for  meaaurlng  overpotentials  at  three  dif- 
ferent positions;  H,  Hg/HgO  reference  electrodes  (40  vtZ 
KOH  solution);  M,  separators;  0,  resistor;  PS,  constant  cur- 
rant power  supply;  R,  two  channel  recorder;  Z and  Za,  sine 
test  and  sine  auxiliary  electrodes  (40  wtZ  KOH-saturated 
slncate) .respectively. 
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Some  electrodes  were  sectioned  using  the  technique  developed  by 
Bro  and  Kang  (39).  The  discharged  electrode  was  fitted  into  a fixture 
having  a microadvancing  screw.  This  allowed  the  electrode  to  be  ad- 
vanced 0.025  inches  per  revolution.  The  electrode  having  thickness  of 
0.1  to  0.11  cm  was  then  machined  off,  by  use  of  a milling  machine,  in 
the  direction  parallel  to  the  electrode  surface.  This  sectioning  opera- 
tion resulted  in  about  5 to  10  sections.  The  machined  chips  per  each 
section  were  collected  and  analyzed  to  determine  the  amount  of  zinc  and 
zinc  oxide. 

The  standard  technique  for  analyzing  zinc  and  zinc  oxide  with  EDTA 
titration  has  been  described  in  detail  by  Welcher  (90) . The  sectioned 
chips  were  dissolved  and  stirred  vigorously  in  10  ml  of  1M  NH4C1-NH40H 
buffer  solution  for  five  minutes.  Previous  investigators  (38,91)  have 
suggested  that  this  solution  dissolves  ZnO  completely  but  attacks  zinc 
negligibly  during  this  period.  The  solution  was  vacuum  filtered  and 
washed  six  times  with  distilled  water.  About  5 ml  of  1M  NH^Cl-NH^OH  was 
added  to  the  filtrate  and  the  resulting  solution  was  titrated  with  0.1M 
disodium  EDTA.  Near  the  end  point,  0.01M  disodium  EDTA  was  used  to  give 
a more  accurate  titration.  The  indicator  used  was  0.2%  solution  of 
Eriochrome  Black  T in  triethanolamine.  The  color  of  indicator  changes 
from  red  to  blue  at  the  end  point.  The  weight  of  ZnO  equivalent  to  1 ml 
of  0.1M  disodium  EDTA  is  8.138  mg.  Following  the  analysis  of  the  ZnO, 
the  filtered  zinc  powder  was  dissolved  in  1 to  2 ml  of  10M  HC1  solution 
for  several  hours.  The  solution  was  adjusted  to  ptW  by  adding  1M  KOH 
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solution  and  then  pH«10  by  adding  1M  NH^Cl-NH^OH  solution,  and  titrated 
with  0.1M  (and  0.01M)  disodium  EDTA.  The  equivalent  weight  of  zinc  for 
1 ml  of  0.1M  disodium  EDTA  is  6.S38  mg. 

The  error  associated  with  the  above  procedure  was  estimated  from 
the  following  experiments.  Samples  containing  known  weights  of  Zn  or 
ZnO  were  titrated  without  filtration.  For  ZnO  the  results  of  the  titra- 
tion agreed  with  the  known  weight  of  the  ZnO  sample  within  0.5%.  For 
the  sampled  zinc  powder,  titration  error  was  within  1%.  However,  when 

j 

a sample  containing  a mixture  of  known  amounts  of  Zn  and  ZnO  is  analyzed 

by  the  above  procedure,  it  appears  that  the  vigorous  stirring,  used  to 

ensure  the  complete  dissolution  of  ZnO,  increases  the  dissolution  rate 

of  Zn.  In  order  to  quantify  the  error  caused  by  the  dissolution  of  Zn, 

samples  of  amalgamated  zinc  powder  were  vigorously  stirred  in  10  mi  of 

1M  NH.C1-NH.0H  solution  for  five  minutes  and  vacuum  filtered  for  one  or 
4 4 

two  minutes.  The  filtered  solution  was  titrated  with  disodium  EDTA. 

It  was  observed  that  the  dissolution  of  zinc  in  1M  NH.C1-NH.0H  solution 

4 4 

was  about  4 wt%  of  the  initial  zinc  within  6 or  7 minutes  contact  time 
as  shown  in  Figure  16.  Thus,  when  the  amount  of  dissolved  zinc  during 
selective  dissolution  of  ZnO  from  a Zn-ZnO  mixture  is  accounted  for, 
titration  errors  remain  within  1%  unless  the  sampled  chips  require  less 
than  about  0.02  ml  of  0.1M  disodium  EDTA  for  titration. 

The  reaction  profiles  are  represented  by  the  absolute  volume  frac- 
tions of  Zn  and  ZnO  as  a function  of  distance.  During  discharge,  the 
electrode  thickness  generally  increased  as  much  as  10%  depending  on  the 
porosity  and  depth  of  discharge.  Therefore,  the  fractional  thickness  of 
the  sectioned  electrode  was  taken  as  the  moles  of  Zn  and  ZnO  in  the  chip 
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Figure  16.  Weight  Percent  of  the  Dissolved  Zinc  In  1M 

NH/C1-NH/0H  Within  6 or  7 Minutes  Contact  Tine 


sample  divided  by  the  total  moles  of  Zn  and  ZnO. 
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where  n?n  and  n^”°  are  the  number  of  moles  of  Zn  and  ZnO  in  the  i th 
sectioned  chip  sample,  respectively,  N is  the  total  number  of  sections, 
and  is  the  dimensionless  thickness  of  the  i th  sectioned  sample. 


V.  RESULTS 


V-l.  Reaction  Profiles 

The  reaction  distribution,  in  terms  of  the  volume  fractions  of  Zn 
and  ZnO,  is  shown  in  Figures  17  to  23.  In  all  of  these  experiments,  a 
dynel  cloth  was  placed  between  the  test  electrode  surface  and  the  mem- 
brane. This  arrangement  is  designated  by  "ZDM"  configuration  through- 
out this  report.  The  measured  volume  fractions  of  ZnO  are  represented 
by  the  histograms  and  those  of  zinc  are  plotted  as  circles  and  hexagons, 
which  are  located  at  the  center  of  each  level  of  the  ZnO  histogram. 

Some  of  the  measured  profiles  are  compared  with  the  profiles  pre- 
dicted from  the  mathematical  model.  The  data  of  Appendix  D-l  are  used 

for  the  numerical  calculations.  Some  of  the  important  parameters  are 

o 2 3 

as  follows:  the  initial  specific  surface  area,  a =300  cm  /cm  ; initial 
KOH  concentration,  c°»10  M;  initial  porosity,  e°=0.33  for  the  low  poro- 
sity electrode  and  e°»0.6  for  the  high  porosity  electrode;  initial 
volume  fractions  of  ZnO  and  the  inert  conducting  material,  E°nQ«0  and 
e°-0;  exchange  current  density,  i°«0. 1 amp/cm2  (taken  from  reference 
52);  tortuosity  factor,  t*0.5;  transport  parameters  of  hydroxide  ions 
across  the  two  layers  of  RAI  P2291  membrane,  t®*0.5  and  1.5x10  “ 
cm/sec  (taken  from  reference  87) . 

All  the  parameters  used  are  taken  from  experimental  data  reported 
in  the  literature  or  are  estimated  from  theoretical  considerations. 

One  exception  is  the  effective  conductivity  of  the  electrolytic  solu- 
tion which  was  chosen  to  give  the  best  fit  between  experimental  and 
theoretical  results.  The  effective  conductivities  are  taken  as  27% 
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(for  the  low  porosity  electrode)  and  12%  (for  the  high  porosity  elec- 
trode) of  the  conductivity  of  the  free  electrolytic  solution.  This 
choice  may  not  be  unreasonable  considering  the  observations  of  Simon- 
sson  (28)  and  Romanova  and  Selitsky  (76) . They  observed  an  effective 
conductivity  which  is  only  about  one  tenth  of  the  conductivity  of  the 
free  solution.  They  attributed  this  large  reduction  in  conductivity 
to  the  gas  generaged  in  the  pores. 

Figure  17  shows  the  volume  fractions  of  zinc  (circles  and  hexagons) 

and  zinc  oxide  (histograms)  measured  for  low  porosity  electrode  (ZDM 

2 

configuration)  discharged  to  8%  depth  at  50  mA/cm  . The  predicted 

profiles  are  shown  as  the  smooth  curves.  The  reaction  profiles  are 

highly  nonuniform  resulting  in  a very  thin  reaction  zone.  The  observed 

as  well  as  predicted  reaction  profiles,  discharged  at  a low  current 

2 

density  of  20  mA/cm  , are  shown  in  Figure  18.  Circles  and  solid  line  A 

are  measured  after  10%  depth  and  hexagons  and  solid  line  B are  measured 

2 

after  23%  depth  of  discharge.  Even  for  the  20  mA/cm  discharge,  a 
thin  reaction  zone  is  still  observed.  After  23%  depth  of  discharge 
(4.51  hour),  all  the  zinc  at  the  electrode  face  is  used  up  and  reaches 
the  region  for  pore  plugging.  Discharge  failure  occurs  at  about  22% 
depth  in  the  theoretical  calculation  but  at  about  32%  depth  in  our 
experiments.  Cracking  or  swelling  of  the  electrode  observed  during 
operation  allows  further  discharge  to  32%  depth  (6.3  hour  discharge) 
before  the  electrode  fails  due  to  pore  plugging.  It  is  to  be  noticed 
that  only  one  fifth  of  the  electrode  is  effectively  utilized  prior  to 
electrode  failure. 

The  zinc  and  zinc  oxide  profiles  for  the  high  porosity  electrode 
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Volume  Fractions 


Fractional  Distance  From  Backing  Plate 


Figure  17.  Profiles  of  Zu  end  ZnO  for  Che  Lov  Porosity 
Electrode  Discharged  to  8Z  Depth  at  50  mA/c»  .Q  andQ, 
Measured  Profiles  for  Zn;  Histograms,  Measured  Profiles  for 
ZnO;  Smooth  Curves,  Predicted  Profiles  for  Zn  and  ZnO. 
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Figure  18.  Profiles  of  Zn  (Circles  and  Hexagons)  and  ZnO 
(Histograms)  for  the  Low  Porosity  Electrode  Discharged  at  20 
mA/cm^.Qand  Histogram  A,  Measured  Profiles  After  10Z  Depth; 

and  Histogram  B,  Measured  Profiles  After  23Z  Depth;  Smooth 
Curves,  Predicted  Profiles  After  10Z  and  23Z  Depths  of  Dis- 
charge. 

■ I 


£0 


I 


2 

to  10%  depth  at  SO  mA/cm  are  shown  in  Figure  19.  Further  discharge  up 
to  21%  depth  (one  hour  discharge)  results  in  the  profiles  in  Figure  20. 
Two  sets  of  data  as  well  as  the  predicted  profiles  are  presented.  For 
the  high  porosity  electrode,  ohmic  loss  as  well  as  concentration  over- 
potentials through  the  solution  are  smaller  than  those  for  the  low  poro- 
sity electrode.  This  small  solution  resistance  shifts  the  reaction 
profile  towards  the  backing  plate  to  produce  a more  uniform  profile . 
Agreement  between  the  measured  and  the  predicted  profiles  becomes  poor 
compared  to  that  for  the  low  porosity  electrode.  This  is  due  to  the 
severe  swelling  of  the  high  porosity  electrode  which  is  weak  in  mecha- 
nical strength.  Evidently,  porosities  for  this  electrode  are  ?arge 
enough  to  eliminate  the  failure  due  to  pore  plugging,  at  j east  in  the 
short  cycle  range. 

2 

In  Figure  21,  the  measured  profiles  for  20  mA/cm  discharge  are 

shown.  Curve  A is  for  10%  depth  of  discharge  and  curve  B is  for  22% 

depth  of  discharge  (2.6  hour  operation).  Since  the  mechanical  strength 

of  this  electrode  discharged  to  22%  depth  was  very  weak,  the  electrode 

was  broken  into  pieces  during  sectioning  process  of  the  missing  portion 

2 

of  solid  line  B.  Discharge  of  this  electrode  at  20  mA/cm  continued 

to  more  than  68%  depth  (8  hour  operation) . 

In  Figure  22,  the  measured  profiles  of  zinc  and  zinc  oxide  for 

the  low  porosity  electrode  having  thickness  of  0.05  cm  is  compared  to 

that  for  the  electrode  having  thickness  of  0.1  cm. 

The  profiles  shown  in  Figure  23  are  measured  after  one  cycle  for 

2 

a low  porosity  electrode  with  an  applied  current  density  of  50  mA/cm  . 
The  total  charge  passed  at  each  half  cycle  was  equivalent  to  10%  conver- 
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Figure  20.  Profiles  of  Zn  end  ZnO  for  the  High  Porosity 
Electrode  Discharged  to  21Z  Depth  at  50  mA/cnr . Q and 
Measured  Profiles  for  Zn;  Histograms,  Measured  Profiles  for 
Zr.O;  Smooth  Curves*  Predicted  Profiles  for  Zn  and  ZnO. 
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Fractional  Distance  From  Backing  Plate 


Figure  21.  Profiles  of  Zn  (Circles  end  Hexagons)  and  ZnO 
(Histograms)  Measured  for  the  High  Porosity  Electrode  Dis- 
charged At  20  mA/cm^. Q and  Histogram  A for  10Z  Depth  of  Dis- 
charge;^ and  Histogram  B for  22Z  Depth  of  Discharge. 
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Figure  22.  Profile*  of  Za  (Circles  sad  Hexagons)  and  ZnO 
(Histograms)  Measured  for  the  Low  Porosity  Electrode  Discharged 
At  50  aA/ca*.Qaad  Histogram  A for  20Z  Depth  of  Discharge  for 
the  Electrode  Saving  Thickness  of  0.05  and  Histogram  B 

for  8Z  Depth  of  Discharge  for  the  Electrode  Having  Thickness 


FRACTIONAL  DISTANCE  FROM  BACKING  PLATE 


Figure  23.  Prof 11m  of  Zn  (Top  Lin«)  and  ZnO  (Bottom 
Lina)  Measurad  for  eha  Low  Porosity  Elaetroda  After  One 
Cycle.  At  Each  Half  Cycle,  30  mA/cm2  and  10%  Depth  of 
Initial  Zinc  ware  Applied. 
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sion  of  initial  zinc.  Some  of  the  ZnO  precipitated  on  discharging  was 
not  dissolved  completely  on  charging  but  remained  near  the  face  of  the 
electrode.  This  type  of  ZnO  movement  on  cycling  might  cause  long  time 
failures. 

V-2.  Overpotentials  of  the  Zinc  Electrode 

The  potential  data  taken  during  the  preliminary  test  period  were 
not  reproducible.  It  was  found  later  that  the  reproducibility  of  the 
data  was  successfully  improved  when  the  amount  of  electrolytic  solution 
accessible  to  the  electrode  was  carefully  controlled.  Consequently, 
several  cell  configurations  having  different  amounts  of  electrolytic 
solution  have  been  used  in  the  potential  measurements. 

Figure  24  illustrates  the  anodic  overpotentials  during  discharge 

2 

at  SO  mA/cm  for  four  different  configurations.  They  are  a low  poro- 
sity electrode  having  dynel  on  top  of  the  membrane  (designated  by  L50- 
ZMD) , a low  porosity  electrode  having  dynel  beneath  the  membrane 
(L50-ZDM) , a high  porosity  electrode  having  dynel  beneath  the  membrane 
(H50-ZDM) , and  a low  porosity  electrode  having  dynel  but  without  a mem- 
brane (LSO-ZD).  For  all  the  cases,  dynel  and  membrane  refer  to  the 
use  of  one  layer  of  dynel  cloth  and  two  layers  of  RAI  P2291  membrane, 
respectively.  Placing  dynel  cloth  below  the  membrane  in  ZDM  configura- 
tion provides  an  additional  solution  reservoir  of  about  0.03  cm3,  but 
the  dynel  cloth  used  in  an  arrangement  ZMD  simply  serves  as  a support- 
ing material  to  minimize  the  difficulties  associated  with  swelling  of 
the  electrode. 

The  zinc  electrode  compartment  is  about  0.12  cm  thickness  (see 
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Figure  14),  while  the  zinc  test  electrode  has  thickness  of  about  0.1  to 

0.11  cm.  For  0.1  cm  thickness  electrodes,  the  configurations  ZMD  and 

ZDM  using  low  porosity  electrode  (e»0.3)  have  electrolytic  solution  of 
3 3 

about  0.05  cm  and  0.08  cm  , respectively.  The  total  capacity  for  the 

the  low  porosity  electrode  is  23.6  A-min  (or  471  minutes  operation  at 
2 

50  mA/cm  ) and  that  for  the  high  porosity  electrode  is  14.05  A-min  (or 

2 

281  minutes  operation  time  at  50  mA/cm  ).  Discharge  failure  occurs 
after  33  minutes  for  L50-ZMD  (curve  A) , after  62  minutes  for  L50-ZDM 
(curve  B) , and  afteT  70  minutes  for  H50-ZDM  (curve  C) . These  observa- 
tions indicate  that  the  electrode  failure  is  caused  by  KOH  depletion 
since  the  failure  time  is  approximately  doubled  by  increasing  the 
electrolyte  volume  from  0.05  cm3  for  ZMD  to  0.08  cm3  (for  L50-ZDM  con- 
figuration) or  0.11  cm3  (for  H50-ZDM  configuration).  The  theoretical 
result  shown  in  Figure  4 (curve  C)  also  indicates  that  the  discharge 
to  10%  depth  of  the  cell  configuration  corresponding  to  ZMD  decreases 
the  KOH  concentration  from  8 M to  2 M. 

The  predicted  overpotentials  using  the  data  of  Appendix  D-l  were 
plotted  as  white  circles  for  L50-ZMD  configuration  and  as  black  circles 
for  L50-ZDM  configuration.  The  volume  of  the  external  solution  reser- 
voir, which  is  located  between  the  electrode  surface  and  the  membrane, 
is  taken  as  0.01  cm3  for  ZM)  and  0.05  cm3  for  ZDM  configurations  in 
the  numerical  calculations.  Discharge  for  the  cell  without  a membrane 
(L50-ZD) (curve  D in  Figure  24)  was  continued  to  137  minutes.  Since 
there  will  be  a large  supply  of  KOH  solution  from  the  counter  electrode 
compartment,  the  failure  observed  for  L50-ZD  electrode  might  be  associ- 
ated with  the  pore  plugging. 
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Figure  25  shows  the  anodic  overpotentials  measured  at  the  face  of 

2 

the  zinc  test  electrode  during  discharge  with  20  mA/cm  current  density. 

2 

Those  for  50  mA/cm  discharge  are  shown  for  comparison.  Curve  C (L20- 
ZDM)  and  curve  A and  B (two  sets  of  data  for  L20-ZMD)  give  about  6 to 
7 hours  failure  time  (31  to  36%  depth  of  discharge) . Discharge  of  the 
high  porosity  electrode  (H20-ZDM)  can  continue  further,  beyond  8 hours 
operation  (or  68%  depth  of  discharge) . It  is  to  be  noticed  that  the 
addition  of  dynel  as  a solution  reservoir  (curve  C)  does  not  increase 
significantly  the  failure  time  compared  to  L20-ZMD  cell  (curves  A and 
B)  On  the  other  hand,  the  use  of  a high  porosity  electrode  increases 
substantially  the  failure  time.  Furthermore,  significant  drop  in 
overpotential  instead  of  continuous  increase  was  frequently  observed 
as  shown  in  curve  A at  1.25  hour  discharge  and  also  in  other  curves 
as  smooth  humps,  which  are  caused  by  swelling  and  cracking  of  the  elec- 
trode. These  observations  favor  the  pore  plugging  model,  attributing 
the  failure  during  a 20  mA/cm2  discharge  to  pore  plugging. 

It  was  mentioned  that  the  calculated  discharge  used  in  Figure  18 
failed  at  about  22%  depth  (4.32  hour  operation).  The  overpotential 
obtained  from  the  numerical  calculations  for  Figure  18  are  plotted  as 
white  circles.  The  predicted  discharge  failure  time  is  much  shorter 
than  the  observed  failure  time.  When  the  electrode  is  allowed  to 
retain  its  porosity  of  0.3  during  discharge  in  the  numerical  calcula- 
tions, the  resulting  overpotentials  resemble  the  observed  overpoten- 
tials as  shown  in  the  black  circles. 

The  numerical  calculations  modified  to  include  the  swelling  effect 
(black  circles)  can  be  described  as  follows.  The  reaction  zone  of  the 
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Figure  25.  Overnotentlals  at  the  Face  of  the  Zinc  Electrode  (x-L)  On 

Discharge  at  20  nA/cm^.  and , Measured  Overpotentials;  O and  • > 

Predicted  Overpotentials  for  L20-ZDM  Without  and  With  Corrections  for 
Swelling,  Respectively.  Explanation  of  Figure  Notation  is  Given  in  Figure 


low  porosity  electrode  (initial  porosity  of  0.33)  is  allowed  to  retain 
a constant  porosity  of  0.3.  The  increase  of  the  solid  volume  due  to 
the  precipitated  ZnO  is  transferred  to  the  reservoir  region  (volume  of 
0.01  ~ 0.02  cm3)  between  the  electrode  surface  and  the  dynel  cloth. 

When  the  porosity  of  this  reservoir  region  reaches  0.3,  the  porosity 
of  the  reaction  zone  is  allowed  to  decrease  due  to  the  precipitated 
ZnO,  causing  the  pore  plugging.  The  black  circles  obtained,  based  or 
these  numerical  procedures,  agree  well  with  the  observed  overpotentials. 

Several  photographs  of  the  discharged  electrode  surface  are  shown 
in  Figures  26  and  27.  In  Figure  26,  the  surface  of  the  unused  zinc 
electrode  (picture  A)  is  compared  with  that  of  the  electrode  L50-ZMD 
which  failed  due  to  K0H  depletion.  The  porous  structure  is  still 
visible  in  picture  B.  In  Figure  27,  picture  A shows  the  surface  of  the 
electrode  L20-ZMD  after  6.33  hour  discharge  to  failure.  The  ZnO  preci- 
pitated as  a continuous  coherent  film  causing  the  pore  plugging. 

Picture  B is  the  surface  of  H20-ZDM  electrode  taken  after  8 hour  dis- 
charge.  The  electrode  did  not  fail  during  this  period  and  the  result- 
ing electrode  still  has  an  open  structure.  The  cracks  and  swelling, 
clearly  visible  in  Figure  27,  are  responsible  for  the  irregular  poten- 
tial behavior  as  described  previously  in  the  discussion  of  Figure  25. 

Figure  28  shows  the  overpotentials  measured  at  three  different 
positions  (curves  A,  B,  and  C)  and  the  potential  differences  between 
the  zinc  test  (L50-ZDM  configuration)  and  the  zinc  counter  electrode 
(curve  D).  The  curves  A,  B,  and  C are  the  overpotentials  measured  at 
the  backing  plate,  at  the  face  of  the  electrode,  and  on  top  of  the 
membrane,  respectively.  The  potential  difference  of  80  mV  between 
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gure  26.  Photographs  of  the  Electrode  Surface.  A,  Unusei 
Zinc  Electrode;  B,  L50-ZHD  Electrode  After  Dis- 
charge Failure. 
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lines  B and  C,  which  is  the  potential  drop  across  the  two  layers  of 
RAI  P2291  membrane,  is  comparable  to  the  ohmic  loss  of  about  90  mV  cal- 
culated based  on  the  measured  conductivity  across  this  membrane  (87) . 
The  first  discharge  failed  after  67.5  minutes.  By  allowing  the  cell  to 
rest  for  8 minutes  at  open  circuit,  an  additional  discharge  time  of 
12  minutes  is  obtained  prior  to  failure.  Increase  in  the  rest  period 
leads  to  a longer  time  tc  failure  for  the  second  discharge.  This 
implies  that  KOH  depletion  is  the  failure  mechanism  under  these  condi- 
tions. 

Figure  29  is  the  corresponding  potential -time  curves  during 

continued  charge  of  the  test  electrode.  The  applied  current  density 
2 

was  50  mA/cm  . During  continued  charge,  about  3.9  A-minute  was  * 
recovered,  prior  to  significant  potential  drop,  out  of  3.98  A-minute 
applied  during  discharge.  Curves  A,  B,  and  C are  the  overpotentials  - 
measured  at  the  backing  plate,  at  the  electrode  surface,  and  on  top  of 
the  membrane.  Curve  D is  the  potential  difference  between  the  zinc 
test  and  the  zinc  counter  electrodes.  Curve  E is  drawn  based  on  the 
solution  ohmic  loss  in  the  counter  electrode  compartment.  The  dis- 
tance between  the  membrane  and  the  counter  electrode  surface  (1.1  cm) 
and  the  conductivity  of  10  M KOH  solution  (about  0.51  flf^cm"1)  (84) 
results  in  the  ohmic  loss  of  about  110  mV  which  is  the  gap  C-E.  The 
gap  D-F,  represents  the  anodic  overpotentials  of  the  counter  zinc 
electrode. 

Comparison  of  Figure  28  with  Figure  29  provides  some  useful 
information  of  the  potential  behavior.  The  test  zinc  electrode  has 
0.08  cm3  of  electrolytic  solution  and  porosity  of  0.33,  while  the 
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electrodeposited  counter  zinc  electrode  contains  l.S  cm5  of  electroly- 
tic solution  with  porosity  of  about  0.6  to  0.8.  Notice  that  the  anodic 
overpotential  of  the  test  zinc  electrode  (curve  B of  Figure  28)  is  much 
larger  than  that  of  the  counter  zinc  electrode  (gap  D-E  of  Figure  29) . 
This  consideration  indicates  that  the  high  anodic  overpotentials  of 
the  test  zinc  electrode  consists  mainly  of  concentration  overpotential 
due  to  the  KOH  depleted  solution. 

Another  interesting  point  to  note  is  the  potential  differences 
between  the  electrode  surface  and  the  top  of  the  membrane  (line 
spacings  B-C  of  Figures  28  and  29) . This  gap  B-C  is  the  potential 
drop  across  the  solution  reservoir  located  between  the  electrode  sur- 
face  and  the  dynel,  across  the  dynel,  and  across  the  membrane.  The 
gap  B-C  of  Figure  29  rises  with  time  and  is  much  larger  compared  to 
the  gap  B-C  of  Figure  28  which  is  nearly  constant  throughout  the 
operation.  The  potential  drop  across  the  membrane  of  about  80  mV, 
i.e. , the  gap  B-C  of  Figure  28,  is  used  to  produce  the  gap  C-F  of 
Figure  29.  The  significant  increase  of  the  gap  B-F,  potential  dif- 
ferences across  the  dynel  and  the  solution  reservoir,  with  time 
indicates  that  zinc  deposition  reaction  takes  place  above  the  elec- 
trode surface  and  grows  through  the  dynel. 

The  overpotentials  during  continued  charge  of  the  electrodes 
shown  in  Figure  24  are  plotted  with  time  in  Figure  30.  The  over- 

potentials  measured  at  the  electrode  surface  are  presented.  Charge 
failure  time  observed  was  37  minutes  for  L50-ZMD  electrode,  78  minutes 
for  L50-ZDM  electrode  which  has  excess  solution  reservoir  in  the  dynel, 
and  112  minutes  for  the  high  porosity  electrode.  Most  of  the  coulom- 


blc  charge  applied  during  discharge  Is  recovered  during  the  charge  step 
prior  Co  significant  overpotential  rise.  For  the  cell  without  a mem- 
brane (curve  D) , very  small  cathodic  overpotentials  were  observed 
throughout  the  operation  period.  Zinc  deposits  penetrated  the  dynel 
cloth  and  continued  to  grow  towards  the  counter  electrode  at  a fairly 
constant  overpotential.  No  failure  was  observed  during  charge  of  this 
electrode.  The  dendrite  growth  of  this  type  will  short  out  an  actual 
cell  whenever  the  membrane  is  bridged  during  operation. 

After  disassembling  the  cell,  a photograph  of  the  dynel  was  taken 
from  the  side  opposite  of  the  test  electrode.  Picture  31-A  was  taken 
after  one  discharge-charge  cycle  for  LSO-ZDM  electrode.  Zinc  dendrites 
penetrating  the  dynel  cloth  are  clearly  shown.  Charge  failed  at  this 
stage  due  to  the  depletion  of  zincate  ions.  A photograph  of  the  dynel 
separator  used  for  a discharge-charge  cycle  of  a L50-ZD  electrode  was 
taken  from  the  same  direction  as  shown  in  Picture  31-B.  Zinc  dendri- 
tes, for  the  cell  having  no  membrane,  seem  to  grow  in  the  direction 
towards  the  counter  electrode  until  no  zincate  ions  are  available  in 
the  counter  electrode  compartment.  In  one  experiment,  charge  was 
allowed  to  proceed  until  failure,  but  addition  of  drops  of  fresh 
electrolytic  solution  restored  original  performance. 


} 
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VI.  DISCUSSION 


VI-1.  Failure  During  Discharge 

In  general,  discharge  failure  of  planar  zinc  electrodes  has  been 
interpreted  as  a phenomenon  called  "passivation".  However,  discharge 
failure  of  the  porous  zinc  electrode  involves  interrelated  phenomena, 
and  therefore,  it  appears  to  be  valuable  to  divide  these  interrelated 
mechanisms  into  several  sub-mechanisms  having  more  precise  meaning,  i.e., 
passivation,  pore  plugging,  and  solution  depletion.  ' Each  failure  mecha- 
nism is  defined  as  follows.  Consider  a test  electrode  compartment  which 
contains  solution  inside  the  porous  electrode  (designated  by  "internal 
solution") , and  solution  outside  the  porous  electrode  but  within  that 
compartment  ("solution  reservoir") . Solution  outside  the  test  electrode 
compartment  on  the  other  side  of  the  membrane  separator  which  is  usually 
contained  in  the  counter  electrode  compartment  is  denoted  by  "external 
solution".  Discharge  of  the  porous  electrode  may  fail  in  three  ways. 

(1)  If  the  concentration  of  hydroxide  ions  in  the  solution  reser- 
voir is  severely  depleted  due  to  the  limited  OH"  flux  across  the  separa- 
tor, then  the  internal  solution  is  also  depleted  of  electrolyte  KOH. 

This  causes  passivation  due  to  a decrease  in  solubility  of  ZnO  and  a 
rapid  increase  in  the  rate  of  ZnO  precipitation  on  the  zinc  metal,  in 
addition  to  increased  resistance  losses.  This  phenomenon  is  defined  as 
concentrat ion  depletion. 

(2)  If  the  concentration  of  the  solution  reservoir  is  reasonably 
high,  but  that  of  the  internal  solution  is  severely  depleted  due  to  the 
limited  mass  transfer  through  the  porous  structure  near  either  side  of 
the  electrode  surfaces,  then  the  internal  solution  becomes  depleted  and 
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causes  passivation.  This  phenomenon  is  defined  as  pore  plugging. 

(3)  Even  if  the  KOH  concentration  of  the  internal  solution  as  well 
as  the  solution  reservoir  is  appreciably  high,  passivation  can  occur  by 
the  physical  blocking  of  the  active  zinc  surface  by  the  precipitated  ZnO 
or  by  the  depleted  hydroxide  ion  or  excess  zincate  at  the  zinc  surface 
caused  by  high  local  current  density.  This  phenomenon  is  defined  as 
passivation,  having  the  same  meaning  as  was  used  for  planar  electrodes. 

Failure  causes  or  parameters  which  are  expected  to  eliminate  fail- 
ure now  become  clear.  For  example,  membrane  parameters  can  be  adjusted 
to  allow  faster  0H~  transport  or  volume  of  the  electrolytic  solution 
reservoir  can  be  increased  to  reduce  depletion  of  0H~  concentration. 

High  porosity  electrodes  allow  for  swelling  and  eliminate  pore  plugging, 
and  high  surface  area  reduces  local  transfer  current  density  reducing 
the  chance  for  passivation. 

Vl-2.  Depletion  of  Electrolytes 

The  overall  discharge  reaction  of  the  zinc  electrode  consumes 
hydroxide  ions  and  produces  water.  Hydroxide  ion  concentration  is 
expected  to  decrease  in  the  zinc  electrode  compartment  during  discharge 
if  the  supply  of  hydroxide  ion  from  the  counter  electrode  is  not  suffi- 
cient to  compensate  for  the  consumed  hydroxide  ion.  Depletion  of  ICOK 
becomes  even  more  serious  for  sealed  or  tight  packed  cells  having  mem- 
branes designed  to  minimize  failure  due  to  zinc  dendrite  penetration 
and  shape  change.  Tight  packing  decreases  the  volume  of  electrolytic 
solution  reservoir  and  generally  membranes  used  to  prevent  zinc  dendrite 
formation  also  limit  the  mass  transfer  of  hydroxide  ions  between  the 
two  electrode  compartments.  A.  restricted  amount  of  electrolytic  solu- 


tion  and/or  limited  flux  of  hydroxide  ions  across  the  membrane  can  lead 
to  a situation  in  which  the  zinc  electrode  fails  on  discharge  earlier 
than  if  a large  amount  of  electrolyte  is  available. 

Failure  observed  during  discharge  of  the  zinc  test  electrode  having 

2 

two  layers  of  RAI  P2291  membrane  with  current  density  of  SO  aA/ca  was 
attributed  to  the  KOH  depletion  (see  Figure  24-A,B,and  C) . Theoretical 
prediction  shown  in  Figure  4 (curve  C)  confirmed  these  observations. 
Depletion  of  KOH  strongly  depends  on  the  type  of  the  membrane,  volume  of 
electrolytic  solution  reservoir,  and  porosity  of  the  electrode. 

In  order  to  make  such  information  more  useful  for  industrial 
battery  design,  an . approximate  relationship  between  the  electrode  design 
parameters  and  the  change  in  KOH  concentration  with  time  was  obtained 
as  follows.  It  is  assumed  that  all  the  zincate  ions  produced  by  zinc 
dissolution  are  converted  into  ZnO  and  that  0H~  concentration  inside  the 
electrode  is  completely  mixed.  Under  these  conditions,  the  decrease  in 
KOH  concentration  inside  the  electrode,  allowing  for  hydroxide  ion  flux 
from  the  boundary  layer,  can  be  written  by  the  following  simplified 


conservation  equation. 
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where  V is  the  volume  of  the  electrode  compartment,  A is  the  cross 
sectional  area  of  the  electrode,  e is  the  porosity,  I is  the  superficial 
applied  current  density,  c3  is  the  average  (mixed)  concentration  of  the 
hydroxide  ion  inside  the  electrode,  c°  is  the  hydroxide  concentration  in 
the  counter  electrode  compartment,  and  and  t^1  are  the  mass  transfer 
coefficient  of  potassium  hydroxide  and  transference  number  of  the 
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(87),  electrode  thickness  of  0.1  cm,  and  initial  KOH  concentration  of 


10M  are  used  in  the  calculations.  Curve  A is  for  one  layer  of  membrane 


0.5)  and  porosity  of  0.3.  Curve  B is  for  two 


layers  of  membrane  (L 


and  curve  C is  for  two  layers  of  membrane,  and  porosity  of  0.6. 

2 

The  average  concentration  of  KOH  for  50  mA/cm  discharge  decreases 

to  zero  after  30  minutes  for  the  low  porosity  electrode  (curve  B)  and 

after  one  hour  for  the  high  porosity  electrode  (curve  C) . The  two  fold 

increase  in  porosity  (or  amount  of  KOH  contained  in  the  zinc  electrode) 

doubles  the  concentration  limited  operation  time.  On  the  other  hand, 

2 

for  20  mA/cm  discharge,  the  average  KOH  concentration  inside  the  low 
porosity  electrode  remains  3.14M  KOH  at  infinite  time  of  operation 
(curve  B),  which  means  that  the  flux  of  KOH  from  the  counter  electrode 


compartment  through  the  membrane  can  supply  the  KOH  solution  consumed  in 

the  zinc  test  electrode.  Consequently,  the  predicted  (Figure  32)  and 

the  observed  (Figures  24  and  25)  results  lead  to  the  conclusion  that  at 

2 

high  rate  discharge  (50  mA/cm  ),  KOH  depletion  causes  the  zinc  electrode 

2 

to  fail  prior  to  pore  plugging,  while  discharge  at  low  rate  (20  mA/cm  ) 
continues  until  pore  plugging  causes  zinc  electrode  failure. 


Discharge  of  most  battery  electrodes  is  accompanied  by  changes  in 
pore  sizes  because  of  the  differences  in  the  molar  volumes  of  solid 
reactant  and  solid  product.  Since  the  molar  volume  of  ZnO  (14.51  cm3/ 

T 

mole)  is  larger  than  that  of  zinc  (9.15  cm  /mole),  the  conversion  of 
zinc  into  ZnO  during  discharge  decreases  the  pore  size  and  eventually 
causes  plugging  of  the  porous  structure  unless  sufficient  space  is 


provided  for  the  swelling.  A nonuniform  current  distribution  will 
accentuate  this  problem  at  the  pore  mout  This  phenomenon  has  been 
given  some  attention  recently  by  Gidaspow  and  Baker  (30)  who  emphasize 
porosity  changes  and  pore  plugging  as  a mechanism  of  battery  failure. 

Discharge  failure  due  to  pore  plugging  was  predicted  from  the 
mathematical  model  as  shown  in  the  reaction  profiles  of  Figures  6 and  7 
(see  also  Figure  18).  If  it  is  assumed  that  all  the  zincate  ions  result- 
ing from  zinc  dissolution  precipitate  as  ZnO  on  discharge,  the  zinc 
electrode  having  a relatively  low  porosity  (less  than  0.4)  can  be  easily 
blocked  by  the  precipitated  ZnO.  Failure  due  to  pore  plugging  is  ex- 
pected to  occur  after  only  10%  depth  of  discharge  when  a penetration 
depth  of  0.1  is  considered.  This  pore  blocking  was  confirmed  by  experi- 
mental and  theoretical  reaction  profiles  obtained  for  the  low  porosity 

2 

electrode  discharged  at  20  mA/cm  current  density,  as  shown  in  Figure  18. 

Potential  measurements  support  this  phenomenon.  Failure  time  observed 

2 

during  low  rate  discharge  (20  mA/cm  ) was  about  6 to  7 hours  which  is 
equivalent  to  31  to  36%  depth  of  discharge. 

The  discrepancy  between  the  observed  and  expected  failure  time  is 
primarily  due  to  swelling  and  cracking  of  the  electrode  and  partly  by 
incomplete  precipitation  of  ZnO.  The  test  electrode  used  in  our  study 
was  supported  by  plexiglas  except  the  front  surface  which  had  dynel  cloth 
as  a supporting  material.  Therefore,  the  stress  arising  from  the 
increase  in  solid  volume  during  discharge  causes  swelling  and  sometimes 
cracking  of  the  frontal  surface  which  increases  the  time  required  for 

plugging  by  several  times.  This  phenomenon,  so  called  pore  plugging, 

2 

is  also  observed  during  high  rate  discharge  (at  50  mA/cm  ) of  L50-ZD 
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electrode  having  no  membrane  (see  curve  D of  Figure  24) . 

Complete  pore  plugging  leading  to  zero  porosity  is  not  necessary  to 
cause  failure.  The  electrode  face  having  pore  size  smaller  than  some 
certain  value  can  limit  significantly  the  KOH  supply  into  the  interior 
of  the  electrode,  resulting  in  KOH  depletion  and  then  passivation. 
Simonsson  (28,29)  has  proposed  this  phenomenon,  i.e.,  blocking  of  the 
electrode  surface  in  the  outer  layer  of  the  electrode  by  a discharge 
product  followed  by  acid  depletion  in  the  inner  part  of  the  electrode, 
as  a failure  mechanism  during  high  rate  discharge  of  the  porous  Pt^ 
positive  electrode  in  a lead  acid  battery. 

The  effect  of  pore  plugging  on  concentration  of  electrolyte  in  the 
interior  of  the  electrode  can  be  roughly  estimated  by  use  of  equation 
(56) . Consider  a situation  which  might  occur  after  20%  depth  of  dis- 
charge of  an  L20-ZDM  electrode.  It  is  assumed  that  an  electrode  with  a 
unit  cross  sectional  area  (A*»l  cm?)  consists  of  an  outer  layer  of  ZnO 
which  has  a porosity  of  0.05  and  a thickness  of  0.02  cm  and  an  inner 
layer  of  unattacked  zinc  having  a porosity  of  0.3  and  a thickness  of 

0. 08  cm.  The  following  parameters  are  applied  to  equation  (56):  A«1  cm2, 

V (volume  of  electrolytic  solution  of  the  inner  part) *>0.024  cm3,  I«0.02 
amp/cm2,  t^-O.?,  and  c°  (initial  concentration  of  KOH  outside  the  elec- 
trode) «0.  005  mole/ cm3.  The  mass  transfer  coefficient  across  the  outer 

(Lgg)  is  evaluated  using  the  values  of  D°B»2xlO"5  cm2/sec  and  d>0.02  cm, 

1. e.,  Lgg  ■DgBe1*5/d»1.118xlO“5  cm/sec.  The  OH"  concentration,  c 3, 
calculated  from  equation  (56)  becomes  3.95  M after  8.3  minutes  and  1.63  M 
after.  33  ..minutes.  . 


VI-4.  Passivation 

A number  of  workers  (14,15,16,17)  have  reported  discharge  failure 
associated  with  the  passivation  of  a flat  plate  zinc  electrode  in  alka- 
line solution  as  was  summarized  by  Jolas  (92)  and  Powers  and  Breiter 
(93).  It  is  generally  agreed  that  the  passivated  electrodes  are 
physically  blocked  from  further  dissolution  by  the  presence  of  a passi- 
vating film  of  ZnO  (14). 

Two  possibilities  exist  for  the  formation  of  the  passive  film; 

(a)  direct  formation  on  the  electrode  surface  through  a solid  state 
reaction  and  (b)  indirect  formation  from  a supersaturated  zlncate  solu- 
tion through  a dissolution-precipitation  mechanism  (51,93,94,95,96). 
Further  studies  (94,95,96)  on  planar  zinc  electrodes  have  shown  that 
a white  porous  film  of  ZnO  (type  I)  forms  by  precipitation  from  super- 
saturated solution  before  passivation  occurs,  while  a coherent  compact 
film  of  ZnO  (type  II)  with  color  varying  from  light  grey  to  black, 
which  is  responsible  for  passivation,  appears  directly  on  the  surface 
beneath  the  type  I film.  It  was  also  indicated  that  a type  I film  and 
then  a type  II  film  could  be  formed  in  the  absence  of  convection,  but 
only  the  type  II  film  was  observed  when  the  conditions  for  supersatura- 
tion  were  largely  removed  by  the  presence  of  convection  or  stirring. 

Only  a few  works  (88,94,97)  have  been  reported  on  the  passivation 
behavior  of  the  porous  zinc  electrode.  Breiter  (97)  has  found  that 
only  a small  fraction  of  the  interior  of  the  porous  electrodes  partici- 
pate in  the  electrochemical  process.  Nagy  and  Bockris  (38)  explained 
their  measured  reaction  distribution  by  introducing  the  effects  of 
the  passive  film  on  the  zinc  surface.  Elsdale  et  al.  (88,98)  observed 
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a gradual  increase  in  potential  in  the  vicinity  of  transition  from  the 
active  to  the  passive  region,  in  contrast  to  the  abrupt  increase  in 
potential  for  the  planar  zinc  electrode  (for  example,  references  14  and 
15).  Consequently,  they  defined  a time,  tQ  3,  as  the  time  required 
for  attainment  of  a potential  rise  of  0.3  V from  the  initial  potential 
and  provided  a linear  plot  of  current  density  versus  (tg  3)“®'5. 

It  is  assumed  that  the  two  types  of  passive  film  (i.e.,  porous 
film  and  compact  film)  observed  by  Powers  (95)  are  formed  by  the  follow- 
ing mechanisms  (99) . 


Zn (OH) 2 

— ZnO  ♦ H20 

(a) 

Zn(OH)*  

— ZnO  ♦ H20  + 20H" 

(b) 

Previous  investigators  have  explained  passivation  phenomenon  using 
transport  process  of  zincate  ions  away  from  the  metal  surface.  A very 
useful  interpretation  of  the  passivation  mechanism  can  be  made  if 
hydroxide  concentration  at  the  metal  surface  is  included  in  the  passi- 
vation mechanism. 

If  the  concentration  of  OH*  at  the  zinc  surface  is  severely 
depleted  or  if  the  ratio  of  zincate  to  0H~  exceeds  a certain  critical 
value,  the  intermediate  species,  Zn(OH)2,  of  the  zinc  dissolution  reac- 
tions can  not  dissolve  in  the  KOH  solution  but  produces  a compact  film 
of  ZnO  (type  II)  directly  on  the  solid  surface  by  mechanism  (a).  This 
phenomenon  is  expected  with  very  high  anodic  current  operation.  On  the 
other  hand,  when  the  concentration  of  OH~  at  the  zinc  surface  remains 
reasonably  high  as  expected  with  low  current  dissolution,  the  inter- 
mediate «pecies  Zn(0H)2  dissolves  in  the  KOH  solution  to  produce  a 
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supersaturated  zincate  solution  from  which  a porous  type  I film  of  ZnO 
forms  by  mechanism  (b) . If  the  porous  t>pe  I film  grows  to  reach  a 
certain  thickness,  transport  of  zincate  and  hydroxide  ions  between  the 
zinc  surface  and  the  bulk  in  the  pores  will  be  limited,  causing  forma- 
tion of  the  type  II  film,  which  is  the  passivating  film. 

This  mechanism  explains  the  comson  feature  of  most  previous 
observations  (99)  that  linear  relation  (between  the  current  density 
and  inverse  of  square  root  of  the  passivation  time)  holds  for  short 
passivation  time,  while  for  longer  passivation  time  deviation  from  the 
linear  relationship  with  less  reproducibility  in  data  occurs. 

Popova  et  al.  (100,101)  concluded  that  passivation  occurs  by 
either  (a)  or  (b)  when  zincate  concentration  in  the  vicinity  of  the 
zinc  surface  exceeds  a critical  concentration.  They  suggested  that  the 
molar  ratio  of  zincate  to  hydroxide  ions  at  the  electrode  surface  was 
0.16  at  the  instant  of  passivation.  Other  workers  reported  that  the 
critical  concentration  of  zincate  ion  was  approximately  half  the  hyd- 
roxide ion  concentration  in  the  electrolytic  solution  (16,17). 

The  passivation  effects  on  electrode  failure  included  in  the  pre- 
sent mathematical  model  (see  equations  (46)  to  (48))  are  surface 
coverage  of  zinc  surface  by  the  precipitated  ZnO  (blockage  of  the  2nd 
kind)  which  corresponds  to  type  I film.  Our  model  does  not  include 
passivation  due  to  type  II  film.  Therefore,  predicted  potentials 
increase  gradually  until  KOH  solution  in  the  bulk  of  the  pores  de- 
creases to  zero.  The  abrupt  increase  in  overpotentials  observed  in  the 
present  experiments  (see  Figure  24  and  25)  are  believed  due  to  this 
type  II  passivating  film  when  the  zinc  electrode  compartment  becomes 


severely  depleted  with  KOH. 
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It  was  not  possible  to  observe  precisely  when  a porous  type  I film 
or  compact  type  II  film  started  to  form  since  the  test  electrode  was 
not  visible  and  covered  by  the  membrane  and  dynel  cloth  during  opera- 
tion. Immediately  after  disassembling  the  cell,  the  color  of  the  dis- 
charged electrode  was  grey  in  the  wet  state*  but  it  changed  to  white 
when  the  electrode  was  washed  and  completely  dried.  The  grey  color 
observed  during  discharge  implies  that  the  film  formed  under  KOH  de- 
pleted conditions  is  very  similar  to  the  type  II  film  of  ZnO. 
Considering  that  a separator  was  excluded  from  the  experiments  by 
Elsdale  et  al.  (98)  and  Coates  et  al.  (88),  their  results,  i.e.,  gra- 
dual increase  in  potential  in  the  vicinity  of  passivation,  are  in 
agreement  with  the  present  interpretation  since  the  KOH  concentration 
in  their  experiments  may  not  have  decreased  to  below  the  critical  value 
to  form  a type  II  film. 

The  failure  observed  for  the  H20-ZDM  electrode  (see  curve  D of 
Figure  25)  is  attributed  to  passivation.  Porosity  of  0.6  for  the  H20- 
ZDM  electrode  is  large  enough  to  prevent,  at  least  in  the  initial  first 
cycle,  pore  plugging  (blockage  of  the  1st  kind)  as  shown  in  Figures  19 
to  21.  Complete  conversion  of  zinc  into  ZnO  results  in  about  0.4 
porosity.  Potential  measurement  (Figure  25),  however,  seems  to  indi- 
cate that  the  high  porosity  electrode  also  fails  after  about  8 to  10 
hour  operation  (68  to  85%  depth  of  discharge)  with  a current  density 
of  20  mA/cm  . Discharge  up  to  the  range  of  this  operation  time  might 
cause  type  I film  formation,  i.e.,  surface  coverage  of  zinc  metal  by 
the  precipitated  ZnO.  Type  I film  having  a certain  thickness  can 
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cause  a significant  decrease  in  0H~  ion  concentration,  an  increase  in 
zincate  ion  concentration  at  the  active  zinc  surface,  and  finally  type 
II  passivation. 


Most  of  the  work  concerning  the  failure  initiated  during  charge 
has  been  associated  with  zinc  dendrite  growth  and  short  circuits  which 
have  been  one  of  the  major  problems  after  many  cycles.  Comparison  of 
Figure  30  with  Figure  24  indicates  that  charge  continues  without  signifi- 
cant overpotential  rise  until  most  of  the  coulombic  charge  applied 
during  discharge  has  been  recovered  on  charging.  This  observation  leads 
to  the  conclusion  that  discharge  products  in  any  form  of  zincate,  ZnO, 
or  Zn(0H)2  can  be  converted  easily  into  zinc  on  charging.  However, 
zincate  depletion  as  a failure  mechanism  can  not  be  neglected  simply 
because  of  this  observation. 

During  charge,  zincate  becomes  undersaturated,  and  under  certain 
circumstances  depletion  of  zincate  might  play  an  important  role  as  a 
limiting  factor.  If  solid  ZnO  can  not  dissolve  as  fast  as  zinc  deposits, 
zincate  is  expected  to  be  depleted  inside  the  porous  electrode.  Even 
if  solid  ZnO  dissolves  fast  enough  to  replace  the  depositing  zincate, 
zincate  ions  at  the  zinc  surface  can  be  depleted  due  to  the  limited  mass 
transfer  of  zincate,  especially  at  high  charge  rate.  The  increase  in 
the  zincate  diffusion  length  from  ZnO  source  to  zinc  sink,  resulting 
from  the  high  molar  volume  ratio  of  ZnO  to  Zn  (1.5),  might  accelerate 
this  phenomenon.  The  latter  effect,  namely,  zincate  depletion  due  to 
limited  mass  transfer,  was  not  a dominant  factor  during  charge  of  the 
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present  electrode  with  current  density  up  to  50  mA/cm  . 

Picture  31 -A  implies  that  zinc  deposit  growth,  instead  of  being 
limited  by  depleted  zincate  ion  at  the  zinc  surface,  continues  to  grow 
toward  the  zincate  source  in  the  dynel  cloth.  Similar  behavior  has  been 
observed  by  Hamby  (102)  who  found  the  zinc  dendrite  penetration  through 
the  porous  polyethylene  film  (1/8"  thickness)  which  was  placed,  as  a 
solution  reservoir,  on  the  opposite  side  from  the  counter  electrode  (on 
the  side  of  the  backing  plate) . These  observations  leave  only  one  possi- 
ble reason  for  depletion  of  zincate,  namely,  slow  dissolution  of  ZnO. 

There  is  some  indication  that  the  dissolution  direction  of  chemi- 
cal reaction  (B)  is  slow  (103,104,105)  on  charge;  however,  the  kinetic 
data  has  not  been  established.  The  fact  that  it  takes  one  or  two  days 
even  with  vigorous  stirring  to  obtain  a saturated  solution  from  ZnO 
powder  in  10M  KOH  solution  gives  a rough  idea  how  slow~the  dissolution 
reaction  is.  The  zincate  depletion  of  this  type  may  cause  the  initial 
low  capacity  during  charge  of  an  electrode  made  of  pure  ZnO,  commonly 
used  as  negative  zinc  plates  in  the  construction  of  secondary  batteries. 

Choi  et  al.  (13)  observed  cell  capacity  less  than  30%  of  ZnO  pre- 
sent during  initial  charge  of  a ZnO  electrode.  They  attributed  the 
resulting  low  capacity  mainly  to  the  failure  of  the  positive  silver 
oxide  electrode.  However,  zincate  depletion  seems  to  have  some  effect 
on  the  initial  low  capacity  since  the  electrolytic  solution  contained 
inside  their  electrode  corresponds  to  only  0.37  A-hour  (3.1%  depth  of 

charge)  based  on  the  following  parameters;  cross  sectional  area  of  66 
2 

cm  , electrode  thickness  of  0.16  cm,  porosity  of  0.65,  and  zincate 


concentration  of  1 M. 


On  charge  of  L50-ZD  electrode  (Picture  3L-B),  without  using  a mem- 
brane, the  cell  accepted  charge  without  rapid  increase  in  potential 
until  all  the  zincate  ions  in  the  counter  electrode  compartment  was 
depleted.  The  zinc  reaction  front,  in  dendritic  form,  then  moves  step 
by  step  towards  the  counter  electrode.  Complete  recovery  of  the  dis- 
charge product  on  charge  (compare  Figures  24  and  30)  does  not  prove 
zincate  depletion  due  to  slow  dissolution  of  ZnO.  Consequently,  the 
success  in  interpretation  of  the  charge  failure  outlined  above  depends 
on  further  knowledge  of  whether  the  rate  of  dissolution  of  the  precipi- 
tated ZnO  or,  if  it  exists,  Zn (OH) 2 is  faster  than  that  of  ZnO  powder 
used  in  forming  the  electrode. 

VI-6.  Failure  During  Cycling 

Zinc  electrode  failure  occurring  after  many  cycles  can  be  referred 
to  as  a phenomenon  called  "redistribution  of  solid  zinc  species  (Zn  and 
ZnO)".  Inability  of  reforming  the  solid  zinc  species  in  their  original 
position  during  cycling  causes  zinc  dendrite  penetration  and  shape 
change.  High  solubility  of  ZnO  in  KOH  solution  may  accentuate  the  re- 
distribution of  zinc  species  from  one  part  of  the  electrode  to  the  other. 

For  the  cell  without  a membrane.  Figure  7 shows  that  only  half 
(55%)  of  the  discharge  product  (zincate  ion)  precipitates  as  ZnO  and  the 
remaining  half  diffuses  into  the  counter  electrode  compartment  before 
precipitation.  On  charge,  once  the  zincate  source  in  the  zinc  electrode 
is  depleted,  zinc  deposits  will  grow  towards  additional  zincate  which 
is  transferring  from  the  counter  electrode  compartment.  Zinc  deposits 
formed  under  zincate  limiting  conditions  similar  to  this  are  known  to 
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be  dendritic,  as  shown  in  Figure  31-B  (see  also  references  1 and  2) . 

On  the  other  hand,  if  a membrane  is  used  as  a separator,  95%  of  the  dis- 
charge product  precipitates  as  ZnO  as  shown  in  Figure  6.  Some  of  the 
remaining  5%  of  the  zincate  ion  stays  in  the  zinc  electrode  compartment 
and  the  remainder  is  lost  into  the  counter  electrode.  Use  of  a membrane 
restricts  zincate  movement  away  from  the  test  electrode.  Charge  of 
this  electrode  is  less  likely  to  cause  dendrite  growth  through  the  sepa- 
rator because  the  source  of  zincate  ions  remains  in  the  zinc  test  elec- 
trode. Comparison  of  these  two  examples  shows  an  effect  of  the  membrane 
on  zinc  dendrite  penetration. 

Redistribution  of  solid  zinc  and  zinc  oxide  from  one  part  of  the 
electrode  to  the  other  can  occur  through  the  movement  of  zincate  ions. 
The  most  important  processes  involving  zincate  ions  include  the  charge 
transfer  reaction,  mass  transfer  of  zincate,  and  chemical  precipitation 
or  dissolution  of  ZnO.  Therefore,  current  distribution  which  character- 
izes the  charge  transfer  reaction,  degree  of  super-  (or  under-)  satura- 
tion which  determines  precipitation  or  dissolution  of  ZnO,  and  mass 
transfer  modes  (e.g.,  diffusion,  migration,  and  convection)  will  govern 
the  redistribution  of  solid  zinc  species. 

Let  us  consider  shape  change  (redistribution  of  the  so lie  zinc 
species  in  the  y-direction  parallel  to  the  electrode  surface)  of  a 
vented  porous  electrode  having  a solution  reservoir  on  top  and  a mem- 
brane directly  attached  on  one  side.  Choi  et  al.  (12)  has  proposed 
that  shape  change  is  caused  by  convective  flows  driven  primarily  by 
membrane  electroosmotic  effects.  Discharge  induced  convective  flow  of 
electrolytic  solution  towards  the  counter  electrode  results  in  a 


de.rease  of  the  total  volume  of  electrolytic  solution  inside  the  zinc 
electrode.  If  the  electrode  is  vented,  the  solution  in  the  reservoir 
will  flow  downwards.  The  solution  picks  up  zincate  ions  along  the  flow 
path  causing  greater  precipitation  rates  in  the  direction  of  flow.  On 
the  other  hand,  on  charge,  the  increase  of  total  volume  of  electrolytic 
solution  inside  the  zinc  electrode  will  push  the  undersaturated  solution 
towards  the  reservoir.  Notice  that  the  solution  in  the  reservoir  now 
becomes  undersaturated  and  greater  dissolution  is  expected  near  the  vent 
or  reservoir.  Repeated  cycling  causes  movement  of  solid  zinc  species 
from  near  the  reservoir  portion  of  the  electrode  to  positions  remote 
from  the  reservoir.  If  the  electrode  is  sealed,  convective  flow,  and 
thus  shape  change  will  be  limited. 

Similar  treatment  can  be  applied  to  the  redistribution  of  solid 
zinc  species  in  the  x-direction  perpendicular  to  the  electrode  surface. 
The  present  work,  in  both  experimental  and  theoretical  aspects,  deals 
with  a sealed  electrode  to  minimize  the  convective  flow  driven  by  the 
membrane.  Under  this  condition,  the  convective  flow  arising  from  the 
volume  changes  due  to  the  overall  electrode  reaction  is  on  the  order 
of  10~7  cm/sec  and  is  considered  negligible. 

The  reaction  or  current  distribution  in  the  x-direction  during 
discharge  is  highly  nonuniform  as  confirmed  by  Figure  5 and  Figures  17 
to  22.  On  the  other  hand,  current  distribution  on  charging  is  more  uni- 
form compared  to  that  during  discharge,  since  precipitated  ZnO  reduces 
the  effective  zinc  surface  area  and  shifts  the  reaction  towards  the 
backing  plate  (refer  to  Figure  9 for  detailed  comparison) . On  repeated 
cycling,  this  difference  of  anodic  and  cathodic  current  distributions 


together  with  differences  in  the  degree  of  super-  (or  under-)  satura- 
tion causes  the  redistribution  of  solid  zinc  species  and  thus  signifi- 
cant distortion  of  the  electrode  structure  occurs  as  cycling  progresses. 
The  distorted  structure  is  shown  in  Figures  11  to  13  (theoretical)  and 
also  in  Figure  23  (experimental) . The  results  are  shown  as  an  accumula- 
tion of  ZnO  st  the  electrode  surface.  If  convective  flow  due  to  mem- 
brane pumping  Is  allowed,  more  serious  redistribution  is  expected. 
Similarly,  allowing  free  movement  of  zincate  ion  into  the  counter  elec- 
trode, by  eliminating  the  membrane,  may  accelerate  this  effect  as 
shown  in  Figures  11  and  12. 

The  result  of  Figure  10,  namely,  substantial  overpotential  rise 
with  repeated  cycling,  predicts  that  the  redistribution  of  zinc  and  zinc 
oxide  in  the  x-direction  is  a limiting  factor  in  the  loss  of  cell  capa-„ 
city  observed  for  cells  designed  to  minimize  shape  change  and  dendrite 
penetration. 

Comparison  of  Figure  28  with  Figure  29  shows  that  the  present  cell 
exhibits  much  higher  overpotentials  during  discharge  compared  to  those 
during  charge.  Similar  behavior  has  been  observed  by  Hamby  (102).  He 
has  concluded  that  the  overpotentials  during  charge  remain  essentially 
constant  but  those  on  discharging  increase  significantly  with  large 
differences  in  overpotential  between  different  positions.  These  observa- 
tions imply  that  failure  during  anodic  discharging  may  control  the 
overall  zinc  electrode  failure.  This  means  that  the  redistribution  of 
zinc  species  may  cause  overpotential  rise  through  discharge  failure 
mechanisms  such  as  pore  plugging,  KOH  depletion,  and  passivation.  The 
ZnO  accumulated  for  many  cycles  can  decrease  pore  size,  at  least  locally, 
resulting  in  high  ohmic  overpotential  and  eventually  deplete  the  elec- 


trolytes  in  the  interior  of  the  electrode  to  increase  concentration 
overpotential.  It  also  reduces  the  active  zinc  surface  area,  promoting 
high  activation  overpotential.  Similar  adverse  effects  are  expected 
due  to  the  gas  trapped  inside  the  porous  electrode. 
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VII.  CONCLUSION 


A mathematical  model  has  been  developed  using  the  concentrated 
ternary  electrolyte  theory  to  predict  the  local  charge  transfer  current, 
current  density  in  the  solution,  potential  in  the  solution  and  in  the 
matrix,  concentrations  of  hydroxide  ions  and  zincate  ions,  porosity, 
volume  fractions  of  Zn  and  ZnO,  and  bulk  solution  velocity  as  a function 
of  time  and  position  perpendicular  to  the  surface  of  the  electrode. 

The  predicted  behaviors  have  been  confirmed  by  galvanostatic  experi- 
ments on  overpotentials  and  reaction  distributions.  The  agreement  be- 
tween the  theoretical  predictions  and  the  experimental  observations 
provides  a good  interpretation  of  the  failure  mechanisms  of  negative 
zinc  electrodes  for  use  in  secondary  batteries. 

2 

During  discharge  at  high  rate  (at  50  mA/cm  ) for  both  low  porosity 
and  high  porosity  electrodes,  the  utilization  of  zinc  is  s f.r  ;ely 
limited  by  depletion  of  electrolyte  within  the  zinc  electrode  compart- 
ment, depending  on  the  volume  of  solution  reservoir  and  the  type  of  mem- 
brane used. 

2 

At  lower  current  densities  (e.g.,  20  mA/cm  ),  the  utilization  of 
zinc  is  much  higher  than  that  at  high  current  densities.  However,  pore 
plugging  plays  an  important  role  in  limiting  the  discharge  capacity. 

The  predicted  and  the  observed  reaction  profiles  on  discharge  are  highly 
nonuniform  and  the  reaction  zone,  located  near  the  electrode  surface, 
is  very  thin,  typically  0.2  mm.  This  highly  nonuniform  reaction  distri- 
bution accentuates  the  failure  due  to  pore  plugging  at  the  pore  mouth. 

Low  current  density  discharge  of  the  high  porosity  electrode,  designed 
to  remove  failure  due  to  pore  plugging,  continued  until  the  electrode 
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failed  due  to  passivation. 

The  predicted  reaction  distribution  on  charge  become  more  uniform 

] 

compared  to  that  on  discharge.  The  difference  of  anodic  and  cathodic 
reaction  distribution  caused  the  redistribution  of  solid  zinc  and  zinc 
oxide  species  on  repeated  cycling.  This  movement  of  solid  species  is 
predicted  to  contribute  to  the  loss  in  cell  capacity  for  long  cycled 
cells. 

The  results  of  this  work  can  be  applied  to  the  design  of  industrial 
zinc  electrodes  to  improve  their  performance  and  cycle  life.  The  re- 
sults also  point  out  areas  where  further  research  should  be  emphasizea 
in  order  to  gain  cycle  life  improvements. 
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APPENDIX  A 


A-l.  Ternary  Transport  Parameters  as  Functions  of  Friction  Coefficients 
The  transport  parameters  used  in  equations  (11)  through  (IS)  are 

t 

related  to  the  friction  coefficients  10  ^ s in  the  following  forms 

LM  ■ 7 ( C?H3  * 2Wl3  * C3M1  ) 

LAB  * LBA  ■ ^■(ClCK13*K23^ClC3C2K13*KloJ-C3M1) 

4b  * 7 ( Cfaj*^o^>*lC3<S3^.>4Ck  ) 

*!  ’ * C3«13K2o  * K12K3o>  ) 

A’  ^ ( Wlo'Vu'  * Wlo  - KI2K30)  ) 
-^T-Cf{Wlo'hcW  * VK13Wlo>) 


where 


D - Cl'M3Mo-KL’*2ClC3<Hol'l3*,tloK3o>*C>iMo-,tU> 

E - M0CK12Ku*K12K23.IC13E23)  ♦ ho**>lhi*2s'  * K23K3o(,:12*1C13) 
* ltloK3o«12*,t23)  * hoVio 


Mi  * H Kii 


ci-  2ici 


; i,j»  1,2,3,  and  o 


; i«l,2,  and  3 


(for  example,  Mj  « Ki0+Ki2+K13^ 


It  is  to  be  noticed  that  one  of  the  transference  numbers  is  a de 


pendent  parameter  because  of  the  relation  t°+t°+t°=l. 

A-2.  Friction  Coefficients  as  Functions  of  Six  Transport  Parameters 
"2  C,c. 


K12  ’ - C1C2  T + 


'i3  ~i~3 


-2  c.c. 


K23  - * C2C3  k + 


‘ A(LBBtA  - LABtB  - cl  H) 

- C'C-  ~ * (Wl  ' LABtA  - | H) 

a!b/l  t°.^L  t°  + — L t°  - — H \ 
L I AB  1 CB  BB  A CA  M B CACB  / 

C / 

lo  " “ y ^CALBB_CBLAB^tA  " ^B^’^AB5  tB  ) 

20  " “ ( CcALBB_CBLAB)  (1“  " c“  (cBLAA'CALAB)t:B) 

3o  * ~ ( CCBLAA"CALAB5 C1*  c^  t3)  ' c^  (CALBB'CBLAB5tA  ) 


where 

L * LAALBB  ' LABLAB 


A-3.  Transformation  of  the  Transport  Equations  into  the  Other 
Four  Component  Systems. 

The  transport  equations  (12)  to  (15),  which  are  expressed  based  on 
the  solvent  reference  velocity,  are  rewritten  in  more  convenient  forms 
as  follows: 

N2  tA 

NA  - -X  - - laa%  - LAB%  ♦ 7 1 ♦ Vo 
v2  F 


N t 

nb  ■ -T  ■ - lab%  - lbb7“b  * 7 1 * Vo 

v3  F 


(A. 3.1) 


i > - k 


nF  ' 


(VntA-  r V % * ‘V"*? 

o 


- 7 v %) 


where  new  parameters  defined  by  CZ2V2^  * tB*t3/,^z3v3^  ’ VV'V 

B 

and  sg*s3/v3  include  the  effects  of  charge  of  the  individual  species. 

The  equations  (A. 3.1)  are  transformed  into  several  forms  to  describe 
the  transport  equations  of  similar  four  component  systems. 

Membrane  System 

A system  consisting  of  salt  A (cation  "1"  and  anion  "2") , solvent 

"o",  and  membrane  "m"  is  considered.  It  appears  to  be  convenient  to 

select  the  reference  velocity  as  the  membrane  velocity  v which  is  zero 

m 

rather  than  solvent  velocity  v . Therefore,  the  symbol  "o"  of  equations 
(A. 3.1)  is  replaced  with  the  symbol  "m".  The  remaining  species  "B"  of 
equations  (A. 3.1)  corresponds  to  the  solvent  species  "o"  of  the  membrane 
system.  The  parameter  used  in  the  relations  of  Appendix  A-l  and  A- 2 
Is  taken  as  zero,  l.e.,  C^-z^c^O.  The  resulting  equations  are 
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.<i  - ■ 


..  . r>\  . 


la  ?va  - + — r~  1 + c»vIB 

A A Ao  o ..A_  A m 

v? 


L®  VWi  - L®  Vy  ♦ t®  ♦ c V 

Ao  A o o o p o m 


nF  \ v. 


(A.3.2) 


y ‘“'9  _ 

A * —A5  % * (S0  * "V  7»= 
2 *2V2 


where 


lA  ■ CA  'Klo*  K2o+  K»o)/D'  > D“-<Klo+,t2«'Ht™,><,Cl,«2.«o.)-KL 

L‘o  ■ Vo  ("20*  V/D" 

Lo  ’ % (Klo*  K2o*  Kl»*  S.)/D" 


tm  - -2-  (K.  K.  - K.  K.  )/D® 

o z^c  v lo  2m  2o  lm 

t®  - (K,  (K,  ♦ K,  ♦ K ) + K.  K )/D® 

2 lmv  lo  2o  onr  lo  onr 

nP2nm 

■ (K  .♦  K _♦  K )(Kf,K,  + K,-K-  K.  K ) ♦ K,  K.  (K,  +K,  ) 
Km  ml  m2  mo  12  lo  12  2o  lo  2o  lm  2m  lo  2o 

♦ K,  K (K,_+  K,  ) ♦ K K (K, K_  ) ♦ K.  K_K 
2o  om  12  lo'  im  om  12  2o  lm  2m  om 


The  transference  njaber  t®  is  related  to  t®  by  the  relation  t“+t®=l. 
These  equations  are  the  same  as  those  developed  by  Bennion  (62) . 

Single  Salt-Neutral  Species-Water  System 


A similar  treatment  can  be  applied  to  obtain  the  transport  equa- 
tions of  a system  which  contains  single  salt  (cation  "1"  and  anion  "2"), 
neutral  species  "N",  and  water  "o".  The  velocity  of  the  solvent  is 
chosen  as  the  reference  velocity.  If  the  symbols  "m"  and  "o"  are 
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replaced  with  the  symbols  "o"  and  "N",  respectively,  equations  (A. 3. 2) 
reduce  to  the  following  equations  (for  example,  for  the  current  equa- 
tion) : 


i ■ - k 


ie  / s,  nt°  $ 0 s \ 

' Trie  + ~~A  ' C~  CA)VyA  * (SN+  ntN"  c~  CN)%) 
*\*2  Z2V2  0 ° / 


where 


(A. 3. 3) 


*1  ■ 1 • " <K2o(1CN1*  ^2*  "no^  * K2nW/D 

«S  ■ Tfi  0^20  - Vlo>'D 


D • (lSu*  '.KN2*  W"ol*  K02*  KoN>  - “NO 

The  transference  numbers  t°  and  t°  are  positive,  but  t°  can  be  positive 
or  negative.  Equations  (A. 3. 3)  together  with  the  corresponding  flux 
equations  are  the  same  as  those  developed  by  Tiedemann  (63) . 


A-4.  Ternary  Transport  Parameters  Calculated  from  Two  Sets  of  Binary 
Data 

The  transport  parameters  Q)^  for  binary  electrolytic  solution 
have  been  compiled  by  Chapman  (66) . Some  of  the  data  for  HC1  and  KC1 
electrolytes  are  shown  in  Table  II.  If  Cl”  ion,  H*  ion,  and  K+  ion  are 
denoted  by  species  "1",  "2",  and  "3",  respectively,  a mixture  contain- 
ing 0.025  M HC1  (salt  "A")  and  0.075  M KC1  (salt  "B")  have  the  follow- 
ing values  for 

- 1. 984x1 0”5  , 2>2o-9.258xl0’5,  ®12-1.178x10-6 

2fJ0  - 2.063xl0-5  , a>3o-1.978xl0”5,  ®13-3.311xl0'7 
Cj  » O.lxlO-3  ,c2  • 0.025xl0*3,  c3  - 0.075xl0”3 
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and 


, p = p * A,c.+  A„c_ 
o 1 A 2 B 


Co  " <■  cama-  cbmb)/M0 


where  p is  the  density  of  the  ternary  solution  (pq  is  the  density  of 
pure  solvent,  and  A^  and  are  the  best  fit  parameters  calculated  from 

the  densities  of  two  binary  solutions),  NT  is  the  molecular  weight  of 
the  species  i,  and  the  subscripts  "o",  "A",  and  "B"  represent  water, 
HC1,  and  KC1,  respectively.  The  parameter  is  missing  but  parameter 
2>io  has  two  different  values.  For  the  following  mixing  rule  is 
employed.  . 


cAo  * cAo 


CA*CB 


The  parameter  is  approximated  by  the  following  form. 


x . / Q)  c&i 
23  J 2o  3o 


2.0431x10 


The  minus  sign  is  chosen  due  to  the  repulsion  forces  arising  between  two 
different  cations. 

* » 

The  friction  coefficients  K. . s,  obtained  from  these  six  2>.  . s by 
use  of  equation  (2),  are  substituted  into  the  relations  of  Appendix  A-l. 
The  calculated  transference  numbers  of  individual  ions  are  compared 
to  those  reported  by  Maclnnes  et  al.  (67)  in  Table  11-2  and  Figure  33. 
Calculated  transference  numbers  based  on  the  arbitrarily  chosen  value 
of  ®23*  “*-00  are  a3-so  shown  in  Table  II-2. 
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Table  II.  Measured  and  Calculated  Transference  Numbers 


. * 
Transport  Parameters  2>.  . ' s of  Binary  Electrolytes 


HC1 

KC1 

2)  xlO5 
-o 

xlO5 

+o 

2>,  xlO^ 

T — 

(mole/1) 

(mole/1) 

(cm2/s) 

(cm2/s) 

(cm2/ s) 

0.025 

— 

1.984 

9.258 

1.178 

0.05 

- 

1.964 

9.270 

1.702 

0.075 

- 

1.951 

9.289 

2.058 

- 

0.025 

2.051 

1.969 

0.1776 

- 

0.05 

2.058 

1.974 

0.2653 

- 

0.075 

2.063 

1.978 

0.3311 

* These  data  are  provided  by  Professor  John  Newman,  University 
of  California,  Berkeley. 


Measured  and  Calculated  Transference  Numbers  for 
0.1  M HC1-KC1  Mixture 


(mol/1)  (njol/l) 


Measured  (67) 

Calculated 

Calculated  ^ 

to 

H+ 

K+ 

V 

K+ 

V 

K+ 

0.7456 

0.0503 

0.7579 

0.0404 

0.7577 

0.0405 

0.6198 

0.1242 

0.6286 

0.1164 

0.6281 

0.1167 

0.4109 

0.2477 

0.4051 

0.2486 

0.4045 

0.2490 

where  ®23*  " -J  ®2o  ^3o  and  23*  ” ^0  are  used  for  a and  b 
respectively. 


FERENCE  NUMBERS 


Transference  Numbers  of  Ions  in  a 0.1M  HC1-RC1 
Ternary  Solution.  Data  Points  are  the  Measured 
Values  (67)  and  Solid  Curves  are  the  Calculated 
Values . 


17.7 
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APPENDIX  B 


MASS  TRANSFER  COEFFICIENTS  (kA>  kg,  and  kg) 


It  is  assumed  that  the  diffusion  of  a species  proceeds  through 
a well  defined  diffusion  layer,  <5,  which  is  the  distance  between  the 
solid-solution  interface  and  the  bulk  solution  inside  the  pores.  The 
geometry  of  the  situation  is  shown  in  Figure  34.  The  dimensionless 
distance,  £,  is  defined  by  the  distance  from  the  solid  surface,  y, 
divided  by  the  diffusion  layer,  5.  As  an  example,  mass  transfer  of 
zincate  ions  occurring  during  charge  will  be  considered. 

Zinc  Deposition  by  Charge  Transfer  Reaction 


Electrode  reaction  is  given  by 


Zn  + 40H“ 


Zn(OH)  ^ + 2e 


where  the  stoichiometric  coefficients  have  been  defined  in  equation 
(10)  to  give  the  values  of  s2— 1,  s^-4,  and  n«2.  During  charge,  zin- 
cate ions  in  the  solution  diffuse  into  the  area,  a^,  to  deposit  as 
zinc.  Within  the  diffusion  layer,  the  flux  of  a species  is  given  by 


Pick's  lav 


N2  “-°AA 


6 d£ 


(B.l) 


Species  conservation  equation  within  the  control  volume  at  Z can  be 


written  as 


1((  V(.-V)  > H2  ) - 


(B.  2) 


with  boundary  conditions; 


c2-c2  at  C-0  and  c2-c2  at  5-1. 


tion  of  ZnO  by  Chemical  Reaction.  Dashed  Lines  Represent  the  Direction  of 
Zlncate  Movement  During  Charge. 


Integrating  equation  (B.2)  yields  the  following  concentration 


profile. 


s 

C2  ~ C2 
b s 
c2  ’ C2 


In  1+  ( 


r*>) 


ln  ( ) 

a 

m 


(B.  3) 


The  flux  of  zincate  ion  at  5-0  across  the  metal  area,  a , or  across 

V 

the  interface  area,  a,  becomes 

D..  dc„  D..  a-a  . 

a N (5-0)  - a - ( ) - f—  (c!-cb 

" 2 '1  < « lw  * i »<-r> 


The  charge  transfer  current  per  unit  area  of  metal  surface, J,  is 


related  to  the  zincate  flux  as  follows: 

di„  nF  nF  . 

— • V>  ■ ( ~2  )amS2(?-0)  ' - 71  * kA  (C2-C2> 


where  ^ ^m 

, , o a 

k.  - k,  

A A , , a . 

in(  — ) 

m 


.0  AA 

kA  ’ 

A 6 


(B.4) 


(B.  5) 


Equation  (B.4)  and  equation  (B.S)  are  equal  to  equations  (31)  and 
(44),  respectively. 

Similarly,  mass  transfer  coefficient  fdr  hydroxide  ion  can  be 


obtained  to  give 


a 

i-  -a  D 

k-k?  *—  : 

^ ^ ln(  f-  ) * « 

m 


Dissolution  of  ZnO  by  Chemical  Reaction 


On  charging,  solution  becomes  undersaturated  with  ZnO  due  to 
the  charge  transfer  reaction  (A).  Under  this  condition,  ZnO  will 


dissolve  into  the  solution  by  chemical  reaction.  Phase  change  from 


ZnO  crystal  Into  the  zlncate  Ion  at  the  ZnO  surface  can  be  described  by 


N. 


“xn.  *■?  - 


XTL 


(B.  7) 


where  c®**  Is  the  equilibrium  concentration  of  zlncate  ion  and 
represents  the  rate  constant  of  the  chemical  reaction.  The  surface 

g* 

concentration  is  not  necessarily  the  same  as  the  concentration 

g 

at  the  zinc  surface,  . The  zlncate  species  is  expected  to  diffuse 

g* 

from  the  ZnO  surface  (c^  ) to  bulk  in  the  pores  (c^) . 

Zlncate  conservation  equation  at  5 in  the  boundary  layer  becomes 

(B.8) 


1 [(V  c)n2J  -0 


where  the  flux  of  zlncate,  N^,  refers  to  equation  (B.l). 

The  boundary  conditions  are 

g* 

c2  - c2  at  €-0 

b r , 

c2  - c2  at  e-1 

The  integration  of  equation  (B.8)  with  (B.l)  results  in  the 
following  Concentration  profile. 


C2  " C2 
b a* 
c2  " c2 


(a-a  \ 

1 + ( T2-  )C  ) 


in  (|-  ) 


The  flux  of  zlncate  at  4-0  transferring  away  from  the  ZnO  area 
(a#)  or  Interface  area  (a)  becomes 


(D. . dc.  \ D. . 

•f  tL-t 


a-a. 


in(-~) 

* a 


(C2’C2*)  (B*9) 


It  appears  to  be  convenient  to  express  the  rate  of  dissolution 
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C-l.  Solubility  of  ZnO  in  KOH  Solution  (Reference  85) 


The  notation  c represents  the  saturation  concentration  of  ZnO 
in  KOH  solution. 


-3  3 

(a)  c^oh  — 2.0x10  mole/co 


c - 0 


(b)  2.0xl0-3  < c^Qn  < 9.0xl0'3  mole/cm3 


c8  - -2.1xl0‘4  + 9.75xl0~2  ^ +1.25  (c^)2 


3c 


.-2 


■^—-9.75x10  +2.5ckqh 


"KOH 


-3  3 

(c)  c^qh  i 9.0x10  mole/cm 


c8  - l.OxlO'4  + 3.0xl0~2  cK0H  + 5.0  (Ctrrt„)2 


"KOH 


3c6 


,-2 


JcKOH 


3.0x10  +10.0 


CK0H 
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C-2.  Conductivity  of  KOH  Electrolytic  Solution  (Reference  84) 


(a)  c^n  1 1.5xl0~3 (mole/cm3) 


-4. 


k - 0.03655  + 195.34  (c^  - 1.513x10  ) mho/cm 


9c 


— — - 195.34  mho-cm2 /mole 


KOH 


(b)  1.5xl0~3  < c^  <_  3.7xl0-3  (mole/cm3) 

k - 0.3  + 134.10  (ckoH  " 1,5x10"3^  nho/cm 


9k 


9c, 


134 . 10  mho-cm  /mole 


KOH 


(c)  3. 7xl0~3  < cR(JH  < 9.0xl0-3  (mole/cm3) 

c -0.00645 

k ■ 0.6572  cos  {y  ( — q Q105 H mho/cm 


9k 


9c, 


KOH 


. _ . r7T  , KOH  

- - 31.295  v sin  (2  ( 0.0105 


-0.00645 


) } mho-cm  /mole 


(d)  c^  >_  9.0x10 


-3 


(mole/cm  ) 


k ■ 0.61  - 57.805(croh  - 0.009)  mho/cm 


9k 


3ckoh 


■ - 57.805  mho-cm  /mole 


v»A 


C-3.  Activity  Coefficients  of  Potassium  Hydroxide  and  Potassium 
Zincate 


The  activity  coefficients  of  electrolytes  A(K2Zn(OH)^)  and 
B (KOH)  are  approximated  by  the  following  equation  (reference  68) . 


In  y. 


In  y. 


2 a ‘JT’  m +2m  4m, 

— + 2 e12(-^— ^)  + e13  — 

1 + Ba'yl  3 3 


i + tJH 


+ 612m2  + B13(”l  + *y> 


where  m^  is  the  molality  (mole/Kg  of  water)  of  species  i,  and  sub- 
scripts "1",  "2",  and  "3"  denote  K+  ion,  Zn(OH)^  ion,  and  OH  ion, 
respectively.  The  ionic  strength  I is  defined  by 


- £ 

2 i 11 


1 1 


r t L, 


'o  o 2 


The  coefficients  a and  8^  f°r  potassium  hydroxide  are  taken  from 
the  reported  data  (68,106),  while  the  coefficients  Ba  and  8^2  f°r 

potassium  zincate  were  estimated  from  the  data  of  activity  coeffi- 

) 

cients  of  K^CrO^  (86) . 

<*  - 1.1762 

Ba  - 1.158 


812  - - 0.01425 


0.130 


APPENDIX  D 


D-l.  Input  Data 

ALPHA  Kinetic  parameter  used  in  equation  (33),  (a  -1.5) 

£L 

ALPHC  Kinetic  parameter  used  in  equation  (33),  (ac-0.5) 

o 2 

AO  Initial  effective  Interface  area  (a  , taken  as  50  cm  / 

3 2 3 

cm  for  Figures  2 to  13  and  as  300  cm  /cm  for  Figures 

17  to  20,  24,  and  25) 

2 

CITOT  Applied  superficial  current  density,  amp/cm 

COB  Initial  concentration  of  K0H,  moles /cm\  (c°,  taken  as 

B 

0.008  for  Figures  2 to  13  and  0.01  for  Figures  17  to  20, 
24,  and  25) 

CONMEM  Effective  conductivity  of  the  solution  in  the  membrane 
(tc™-  0.5243  mho/cm) 

DIFA  Diffusion  coefficient  of  potassium  zincate  in  the  solu- 

—6  2 

tion  (D. .,  taken  as  6.86x10  cm  /sec) 

AA 

DIFB  Diffusion  coefficients  of  KOH  in  the  solution  (D__, 

BB 

2x10  ^ cm^/sec  at  10  M KOH  solution) 

DJS  Parameter  characterizing  complete  coverage  of  the  zinc 

surface  by  ZnO  ( 1,  taken  as  zero) 

EPI  Initial  porosity  of  the  electrode  ( £°,  taken  as  0.4 

for  Figures  2 to  13.  For  Figures  17  to  20,  24,  arid  25, 
o o 

e -0.33  and  e -0.6  are  used  for  the  low  porosity  and 
the  high  porosity  electrodes,  respectively) 

EPGI  Initial  volume  fraction  of  the  Inert  conducting 

material  (e°,  taken  as  zero) 
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2 

for  Figures  2 to  13  and  0.1  amp/cm  for  Figures  17  to 
20,  24,  and  25) 

FL  Thickness  of  the  Zinc  electrode  (L-0. 1 cm) 

FLCOUN  Distance  between  the  membrane  arid  the  counter  elec- 
trode (taken  as  1.1  cm) 

GAM  Kinetic  parameter  used  in  equation  (33)  (y  >0.75) 

QI  Initial  volume  fraction  of  ZnO  (e?  taken  as  0.1  for 

£nu 

Figures  2 to  13  and  zero  for  Figures  17  to  20,  24, 
and  25) 

RMO  Mass  transfer  coefficient  of  zlncate  ion  defined  by 

aquation  (44)  (k°  , taken  as  0.001  cm/sec) 

ROHO  Mass  transfer  coefficient  of  hydroxide  ion  defined  by 

equation  (44)  (k^,  taken  as  0.003  cm/sec) 

RXT  Rate  constant  of  the  chemical  reaction  defined  by 

equation  (44)  taken  as  0.005  cm/sec) 

SGMO  Conductivity  of  pure  zinc  metal  (a°n>  2x10^  0 1cm  *) 

SGSO  Conductivity  of  ZnO  (o^nO*  ^cm  *) 

SGGO  Conductivity  of  inert  conducting  material  (Oj,  taken 

as  100  0 *cm  *) 

Si  Stoichiometric  coefficients  of  species  1 for  the  charge 

transfer  reaction,  defined  by  equation  (10) ; SI  fpr  K+ 
ion  (Sj“0) , S2  for  zlncate  ion  (Sg^-l)  > S3  for  OH  ion 
(Sj«  4) , and  SO  for  water  (sq-0) 

SSi  Stoichiometric  coefficients  of  species  1 for  chemical 

precipitation  or  dissolution  reaction;  SSI  for  K+  ion 
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1 


* * _ 

(Sj-O),  SS2  for  zincate  ion  (s2*l),  SS3  for  OH  ion 

„ * * 

(s  *-2),  and  SSO  for  water  (s  * -1). 

J o 

T2R  Trnnsference  number  of  zincate  ion  (t°,  taken  as  0.05 

at  10  M KOH-saturated  ZnO  solution) 

T1R  Transference  number  of  the  potassium  ion  (t^,  taken  as 

0.23  for  10  M KOH-saturated  ZnO  solution) 

TOR  Tortuosity  factor  for  diffusion  coefficients  (TOR)  and 

TOT 

that  for  conductivity  of  the  solution  (TOT)  (t*0,5) 

TR2M  Transference  number  of  zincate  ion  in  the  membrane 

(t“-0) 

TR3M  Transference  number  of  OH  ion  in  the  membrane  (t™*0.5) 

_ 3 

VA  Molar  volume  of  potassium  zincate  (V  -67.0  cm  /mole) 

— 3 

VB  Molar  volume  of  potassium  hydroxide  (V  “17.8  cm  /mole) 

B 

— 3 

VM  Molar  volume  of  metallic  zinc  (V_  * 9.15  cm  /mole) 

zn 

VS  Molar  volume  of  ZnO  (V_  _ *14.51  cm'Vmole) 

ZnO 

— 3 

VO  Molar  volume  of  water  (V  ■ 18.07  cm  /mole) 

V0L1  Partial  molar  volume  of  potassium  ion  (V^+  , taken  as 

3 

7.2  cm  /mole) 

VOLUME  Volume  of  the  electrolyte  reservoir  between  the  zinc 
electrode  surface  and  the  membrane  (V,  taken  as  zero 
for  Figures  2 to  13.  For  Figures  17  to  20,  24,  and  25, 
V Is  taken  as  0.01  cm^  (for  ZMD)  and  0.05  cm^  (for  ZDM)) 
ZETA  Kinetic  parameter  used  in  equation  (33)  (c*0.0) 
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D-2.  Programming  Parameters 


Control  Parameter 

N Number  of  unknowns  (N»6) 

NJ  Number  of  mesh  points  used  (NJ«31) 

g 

JP  Has  the  value  of  0 if  <J>2  is  assumed  to  be  the  same  as 

g 

$2  (see  equation  (29)),  and  1 or  2 if  correction  for  $2 
Is  Included. 

JPLG  Has  the  value  of  1 or  2 if  swelling  effects  are  included. 
Otherwise,  JPLG  is  zero. 

L Specifies  the  initial  mode  of  operation.  Equal  to  1 for 

discharge  and  2 for  charge. 

LMX  Has  the  value  of  1 without  inert  conducting  matrix  and 

2 with  inert  conducting  material. 

LCB  Specifies  the  initial  concentration  of  the  continued 

half  cycle.  Equal  to  zero  if  c^  and  c°  are  used,  but 
equal  to  1 if  average  concentrations  of  salts  A and  B 
resulting  from  the  previous  half  cycle  were  used. 

LPUN  Has  the  value  of  1 if  punched  cards  for  output  are 

desired.  Equal  to  zero  if  punched  data  cards  are  not 
necessary. 

LREP  Activates  reading  of  the  punched  output  cards  (obtained 
from  LPUN*1)  and  continues  calculation 

NACT  Use  of  binary  activity  coefficients  for  the  value  of 

zero  and  ternary  activity  coefficients  for  the  value,  2. 

NBC  Specifies  types  of  boundary  layer.  Equal  to  zero  for 
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NCYCLM 


NT 


TDEL 

TLIMIT 


TSTOP 

Subroutine 


BAND 

MATINV 

EQUIL 


RKAY 


SACTIV 


SCOND 


SDIF 


solution  boundary  and  1 for  the  membrane  boundary 
Specifies  the  number  of  half  cycle.  Equal  to  0 or  1 for 
the  initial  half  cycle.  The  value  of  4 activates  the 
run  for  the  two  complete  cycles. 

Has  the  values  greater  than  zero.  The  first  time  step 
of  each  half  cycle  (onset  of  each  half  cycle)  is  di- 
vided by  the  NT  number  of  intervals. 

Time  step  duration,  seconds 

Fixes  maximum  of  operation  time  of  each  half  cycle, 
seconds 

Specifies  the  rest  period  between  two  half  cycles,  seconds 


Called  at  each  mesh  point  to  solve  the  sot  of  linearized 
equations.  These  are  developed  by  Newman. 

Equilibrium  or  saturation  concentration  of  ZnO  in  KOH 
sclution  and  its  derivatives  with  respect  to  concentra- 
tions , C£  and  Cj 

Active  surface  area  (a  , a , and  a) , mass  transfer 

m s 

coefficients  (k^,  1^,  and  k*)  and  their  derivatives 
with  respect  to  concentrations  , C£  and  c. 

Activity  coefficients  of  KOH  and  potassium  zincate 
and  their  derivatives  with  respect  to  concentrations. 
Conductivity  of  solution  and  its  derivatives  with 
respect  to  concentrations  , and  c^ 

Diffusion  coefficients  of  potassium  hydroxide  and 
potassium  zincate  and  their  derivatives  with  respect  to 
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SDIFP 


concentrations,  and 

Second  derivatives  of  the  diffusion  coefficients  with 
respect  to  concentrations,  and  c^ 

SIGMA  Effective  conductivity  of  the  solid  matrix  and  its  deri- 
vatives with  respect  to  volume  fractions  of  zinc,  zinc 
oxide,  and  inert  conducting  materials. 

Dimensional  Parameter 


A,B,D,G 
X.and  Y 

C(I,J) 


AMET 


AREA 


CAS 


CRATIO 


CONVRN 


EQUL 


Coefficients  of  the  linearized  equations  which  were 
defined  by  Newman 

Variables  I (1-1  to  6)  at  mesh  point  J of  the  present 
iteration  at  the  present  time  step.  CI(I,J)  refer 
to  those  of  the  previous  interatlon  at  the  present  time 
step.  Z(I,J)  represent  the  variables  at  the  previous 
time  step. 

Active  surface  area,  a^,  which  is  active  for  charge 
transfer  reaction  (at  the  present  time  step) 

Effective  solid-solution  Interface  area,  a,  of  the  pre- 
sent time  step. 

Ratio  of  the  concentration  of  zincate  ion  at  the  metal 
surface  to  that  in  the  bulk  of  the  pores. 

Ratio  of  zincate  concentration  to  hydroxide  concentra- 
tion at  the  metal  surface 

Fractional  conversion  of  Zn  Into  ZnO,  l.e. , moles  of  the 
precipitated  ZnO  divided  by  the  moles  of  the  Initial  Zn. 
Zincate  concentration  of  the  bulk  solution  in  the  pores 


1 


minus  the  saturation  concentration  of  zlncate 
Porosity  of  the  electrode  at  the  present  time  step.  ZEP 
(J)  refers  to  that  at  the  previous  time  step. 
Overpotentials  with  respect  to  the  Hg/HgO  reference 
electrode. 

Volume  fractions  of  ZnO  at  the  present  time  step.  ZQ(J) 
are  the  volume  fractions  of  ZnO  at  the  previous  time  step. 
Volume  fractions  of  Zn  at  the  present  time  step.  ZQM(J) 
are  the  volume  fractions  of  Zn  at  the  previous  time  step. 
Volume  fractions  of  Zn  which  is  nonactive  due  to  pore 
plugging. 

Product  of  the  specific  surface  area,  a,  and  the  rate 

constants  for  chemical  reaction,  k*.  ZRS(J)  are  the 
* 

values  of  ak  at  the  previous  time  step, 
s 

Transfer  current  density,  j,  at  the  present  time  step. 

D-3 . Computer  Program 

Computer  program  used  in  this  work  is  presented. 
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EP 

OVERP 

Q 

QM 

QMNON 

RS 

TRNCUR 


SK5£*"*",“ 


IMPLICIT  REAL*0  (A-N.3-Z)  " 

DIMENSION  A (6.6) .at  4.6) .CI6.93) ,0(6.13) .G(6) .X(6 .6) .Y  (6.6) . 

1 Z( 6.93 >, TRNCUR(93) . =»  ( 33).ZEP(93) .0(93). Cl  ( 6. 93 ) . QMNON( 9 3) . 


_2_A9f  A (.93)  .ZQ{  ?3  ) .CAS LIU  ._A9£T (93  ) , RS ( 93  ) ,_S»S ( 93J  . ON ( 93  ) . ZOM ( 9 3 ) 
jTEQULt  93  ) . C0NV9N  ( 93 ) . Ovf9»(  93)  .CRA  TI0  ( 93  ) 

COMMQN/WORKA/A.S. O.S.X.Y 
COMMON/ WORKB/N.NJ 

..  COMMON/PORK C/AO. T t 9E. 99 O.R SO .RXT .CP ( .EPGl  ,E PNUCS.gPNUCM  .EQUI  AS. 
t EQUMAX.EPL1M.N0ZN. MM, J 
COMMON /PORK 0/C 

COMMON/ WORXE/S  GMO  . SCSO  . SGGO  , EMI.  t M 


C0MM0N/90RKP/NACT 

common/ MORKS7o'M,6.SI>.W.vaL'/M£.ei»(:(CC.j»L5 

COMMON/MORKH/TOT.L 
COMMON/ ttOR X t /SPEC  ND 
98  READ  101. NJ 

read- 1 0 a . l . LMX  r NC  V CCV7C R S p.  0>  unTn  tfCV  n TTL  C BTJ  P~ 
REAO  10S.N0ZN.NACT, JPLG.N3T0P 
R EAO  107. N JOEL  T . N J09N , N J09X , EPLGMN . ESPLG 


0 

o"a 

109 

100 


PORMAT( 112) 

“ MATrtiina.au 


r 


CRM) 

PORNAT (All ) 

PORMATI/l 1 I A/) 

107  PORMATt  312. AX. 2PI 0.5) 

i o«  porm at'( sx  .ta-.pzoys  , a r*7^2a'3T • 

IP(NJ.EQ.O)  STOP 
N*6 

— : 

REAO  102 . T ZR.T1R.S2, S3 . SO, SS2. SS3 • SSO 
REAO  102.VA  .VB .VO.SGMO.SGSO.SGGO.RMO.RXT 

REAO  102.P10.AO.EPt.3t.EPGI.OJS  .VCL1.RCH0  

REAO'  10 2 iCALML  . CAl.9R  . POTL  I 9 “.  EPNUCS'.  EPNUCMTEOiJMA  X TOCI E*CTH 

REAO  t02.TRaM,TR3M.(H.AAM.OLS8M.  VELM. VOLUME. EPPCNO 
102  PORM AT (8010. A) 

PITOT»CtTOT 

TPTCTE 0T21  P[TdT*-ctTQT  " - 


tP(LMX.EO.l)  EPGImO.O 

EPG lCC»EPg  I 

QMJ ■! ,0-?P I-QI-EPGI 
TM( -TOEL/200.0 
NJM1»NJ-1 
NJM2*NJ>2 

"P-06A07 . 0 

R-8.31A3 
T»298 .1 6 

EN«2.0 


MMs~l 

TR1 *TIR 
TR2«T2R 

TR3*1.0-TR1-TR2 

"A  ST  A “iT  vM  /EN/P  

BST AR«(  -2  • OA VOU  ♦ S2* VAA  S3* VBaSO* VO  ) /2  .0 /P 
C STAR*  VS 


TIME*0.0 

CALL  EOUIL  < CO B. COB.  CO  A . E'JQ  XB.  ECBQP  ) 

.CALL  SCONOt COA, COB. CONGO. C 3N0AA .CONORS) 

00  l l J-l.NJ  ~ 

0(J)»QI 
EP( J)*EPt 
QM( J)*OMt 

talc  t o-jnrrcdATaifTgigTrrrrr  > — 

EOUL(J)  * EQUtAS 
C(3 • J)*COA 


„ C ( A . J ) — C OB 

TI  CONTINUE 


QMMAX*OMt ♦Ot*VM/VS 
OMPLG-OMJ -EPI • VM/( VS- VM) 

IP(OMPUG.LE.O.O)  Q9PLG*  9.0 

“oweoviojs*  rcMMAir=jriPcn 

0PN0NI aOMPLCAOMCOV 
OMOaQMt-OMNONt 

CALL  St  GMA( ERG  I . EP 1.390, 31.SI COM , S t GOS . S IGM ( ) 

“"CALL  RK  AY ( TOEL  .EPt  .090 .01 . AMO , AMCUR . AMC A, PM  t .PSI . PMCUR » RSCU9. 

1 RMCA.RSCA) 

PRINT  1 06  *NJ,|.,  LMX.NCVCL9.L9EP.LPUN.NM.NBC.NT.LCB.jp 
PRINT  1 06 .NGZN • NAC  X . J PL3 . NSTOP 
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ttis  PA®  IS  BBS1  quality  mctiCABia  _ ; 

1WB  OOPY  FUBA1SHBD  TO  DDC.  *" 

PHI  NT  1 09*  N JOEL  T.  S»LCMN.N  JOHN  . NJ  D_M  XT?  5TOJ 

PRINT  201 . PITOT. PLtC04.C39.TOEL.TLlM. ALPHA. ALPHC. CAM. 2ETA, 

1 TOR.VM.VS. VA.V8.va.TR2.TR3.S2.Sl. 

2 ».SS2.SS3.SS0.FI3.SaI..>MI,0t.EPGI.A0.  

"yRMO.RXT.CONOO.SCMO.SCS3.SCGO.  SI  CM  I .DJSTVCCTTCXLml  . 

4 CALMH.POTL IM. TR2M « TR  )M ,0L AAM .0LS9M. 6PL  t M .EPNUCS • EPNUCM. EQUMAX 
9 RM I ,RS  I.AMO.OMNONf , OMPLC.VELM.VOLUME.ROHC. TOT 


« .PLCOUN.CONME  m ,E  ML I M . gPPCMO 
201  ^ P OHM  ATd  0/  fSPSOTOT) 


iwnnsep'.Ecr^nnMrTTT 

DO  12  J*1.NJ 

REAO  t03.C( l.J) .CC  2.J) .CC3. J)  .CC4.J) .CCS. J) 
READ  103.CC6.J>.EPtJ).0M< J).Q<J).RS< J) 

T53 — roRMTT T3D l 4 J f 

12  CONTINUE 

REAP  103. COMNONt J). J»1 »NJ) 


REAO  103.1  SOUL ( J t • J»1 . N J 1 
NCYC-1 

REAO  204. TIME. OSLT.TLIM.NTI 

REAO  1 03.SUM4.SUM9 

204  PQRMATC 3P1 0 .3. 13) 

IP< LREP  «CE .2 } CO  TO  4 
GO  TO  6 

31 CONX1MUC 

NTI»0 

OELT»TOEL 

NCVC«l 


SUM4*C0A 

suM3«coe 
4 CONTINUE 

00  7 J-l.NJ 

OMMAXaQMt J l*Ot  liaWM/V« 

OMPLC«OM( J ) -EPt  J ) *VM/ ( VS- VM ) 
IP(QMPLC.LE.O.O)  OMPLC*0.0 
QMCOV»OJS*< CMMAX-OMPLG) 

QMNCNC  J J jtQM  PLC >.QMCQ3t 

7 CONTINUE 

« CONTINUE 

« IP( TIME.GT.TMt ) GO  TO  31 

IPlLga-ll  21.21.?? 

21  C AI NT«OPLO AT1LCB ) *SUM4 ♦OPLOAT ( 1 -LCB ) *C0 A 

CBINT»OPLOATC  LCB ) •SUM3+0FL0AT ( 1 — LCB ) 4 CO 8 
GO  TO  24 

—22 CBI NTaOFLOAJ CL.CB— 2 1 4SUMS.TPFL3  AT  C 3— LCB  1 aCQB 

CALL  EOUIL(COA.CBINT.CAINT.EOUAB.ECBOP) 

24  CONTINUE 

C APPROXIMATE  SOLUTION  FOR  TIMC«0.0 


DO  1 J«1.NJ 
C(3. JI-CAINT 
CCA. J)«C8INT 
QMlaOMl J) -OMNONt J I 

"CALL  EOU I Lt  CC3.J)  .C  <4  , J » 346  .fiou  f AB,  EOCBP) 

CALL  SIGMACEPCI «EP<  J» .0M1 .3(  J) . SGM, SGS. SI GMj) 

CALL  RXAY CTCEL.EPC  J1.0M1 , O(J) , AMJ . AMCUR .AMC4 .RMJ .RS ( 31 . RMCUR • 

1 BSCUH.BHCA.flSCA) 

ametc J)aAMJ 
CASC J)«C<3. Jl/COA 

IFCCASC J) .LE.CALML)  CAS ( Jl *CALML 

IILCCAS  (J)  . 6T  .CALMM  1 CASCJlagALMH 

EPJ-EPC J ) 

CALL  SCONO (C( 3 • Jl ,C<4 . J> .CONO J. CNOA .CNOS) 

RPOL*AMJ*FI  0*F*<  ALPH4«-4LP  MC  ) /R/T4C  A SC  J » **GA  M 

BiM<»apoL«M.B/<icij/(  t -a -»t>  1 1*1  .n/fPNnj/n»j««(i  .fltrnrn 

RAMaOSORT CRAMS) 

RKAPaRPOL*PlTOT/S I CM J /RAMS/ C 1 .O-EPJ ) 

EXPRaOEXPC  RAM4FL  >— OEXPC -RAMaPL ) 

_PA-iF.ITOT-RKAP*C  1,  J-Oexat-RAMAPLJJ  LZ*XPJ 

RBa  C -F  ITOT-RK  AP4C  OEXPC  RAM4FL  1-1.0))  /EXPR 
EXaOFLOATCJ-l) /OPLOAT CM J-l )*FL 
CC1 • J)aRA«0exP<RAM*EX)*R3«0exPC -R AM*E X ) 4RKAP 
_ C.C6  . J)  a CR AFOEXPI  R4M4EX ) / RAM  — RBaQFXPI  —RAM  ag X I /RAMaRXAPREX 
1 -F I TOT  *E  X— CRA— RB) /R AM ) /S  t CM  J/( I .O-EPJ) 

CC2.  J)aCC6.  J)-CRAaR4MOEX9CRAM*CX)-R84RAM*0EXPC-RAN»EX)  )/RPOL 
CCS.Jla  -CASTA  ♦BSTAR)4CCt .J) 
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-IBIS 


fia£  IS  BBS!  »««  «>«““ 
JODD.Q  — ^ 


CONTINUE 
CONTINUE 
MINT  217 


OQK  *ik*ishje 


III PORNATC /20X, •OWNON»i/J 

MINT  209.  ( QMNQNt  JJ  . 3«l  .NJJ 


IP(L.EQ.l)  MINT  211 
IP(L.E0.2>  PRINT  212 

n KHtimbi 

IP(LMX.EO.l)  PRINT  213. 53G I 
tPn.WX.EO. 2)  PRINT  21*. Pool 

_213  PORWATC//IM  « 8X.  * POR  ZINC  MATRI X • . 20X  . «EPGI  ■♦.P9.21 
21*  'PORWATC//lH'.ax.*PaR  INERT  WATR  IX»‘.  20'xT*  tfPOl  P3.2] 
PRINT  222 . T tME 

222  PORNATI//  20X. •TIME** .P10.3//1 

8StnTJTo4  7c(  17V  > 7g(f7j>  .£(  1757  .£(*'.  j>.g(S  .j>  7C  TSTj  i 1 

209  PCRMAT ( A020 .8 I 

S CONTINUE 

PRINT  29 

If PORN AT ( /l  SX • 

1 •CAS**/) 

00  29  J»1.NJ 

-JJ 3 ) «QM( J ) «0(J1«RS<3I . AMET  C J I. CAS ( J 1 

2 JCOUNT-O 

DELTIa  ASTA  aOELT/2.0 

OGLaCST  AR  *0ELT/2.  0 

IP(WN.EO.O)  OELaO  .6 
IP(MM.EO.O)  OELTI *0 • 0 
00  3 J«1.NJ 

Z0( J)»0<J 


1'BIBU 


DO  13  I«1.N 
00  13  K«1.N 


IPI3.NE.1I  GO  TO  *0 

C 

t PROGRAM  POP  J«l 

C 

e<2.n«i.o 

ots.si*  1.0 

9(6 .A) a 1.0 

CAP  ao.S*(C(3. J)*C(3.JM> ) 

COP  aO.S*(C(*.J)*C(*. J*l » » 

CURPaO.9* (C < 1 » 3 1 *C( 1 . J* 1 1 ) 

tuRppa  < crrry*i  »^n,h>Mr 

CAPP  a(C(3.3*l »-C( 3.3) >/H 
CEPP  a(C(*.3>l )-C(*. JJJ/M 
PPIPPa (C(2.J*1 )*C  ( 2.3) >/H 

CAUL  9C0N0(  CAP  .C9P  .CONO  .COfWA.'CCNOOl 

IP( TIWC.GT .TMI ) GO  TO  32 
CPPa0.9*(SPC3)  *EP(  ■»♦!  I) 

O taCONO«E  PP»(1.0*TOT1 

o ( i . i y*  a .« — : 

0(1 .1 laO.S 
9(  1 .2) a-OI /M 

0(1.2»aQl/H  


I 


IBIS  FAOB  IS  BSST  QUALITY  FJUCTICABL1 
mm  Qgg  MBH.SHBD  TO DftS  — — 

‘8(3*31*1.0 
8(4*4) a 1*0 
6(3>*CAINT 
8(4) *CBtNT 
~ 60“T0j 4 
CONT 1NUE 

CALL  SACTIV(CAF,CeF,ACTA. 4CT8.YA.YB .Y AA.Y AB» VBA. V 09 ) 
CALL  SHATR<  T1R.T2R .COA .C4F.CSF.PA.PB. PAA.PAB *PBA . PBB) 

QFa ( Q ( J ) 4Q  (J  + l ) )/2  .0  

QMP1*(0M( JJ.QMt JMD/2.0 
OMNONF*(OMNON(  JH-QVINONI  J+l> )/2.0 
OMFaOMFl-QMNONF 

■“Eprrrepi  J7>epf3-n  n/3.0 

CALL  EaulL(CAF.CBF.EJCS.SOUIAB.EaCBP) 

CALL  PKAr (OELT.EPP.OMP.QP. AMF.AMCUP.AMCA.RMP.RSP.RMCUR. 
t RSCU8.RMCA.RSCA) 

“DELTAa-OELaRSP-  

0ELTBa0ELaRSF*SQC8P 
Ol*CONO*EPF**< I .04T3T) 

02a  V A4PA/CAF 


OAAa  C0N0A4Y4A4PAA 
OBAa  CONOA+YBA+PBA 
OASa  CONOB ♦YAB4PA8 

“OBBS  CQNOB  .Y8B+P9(J 

02a  (1.04rOT)«01*(PH(pP402«C4PP4034CePP)/EPP 
O Aa  01*{C3N0A*PHtPF432*(044— 1 • O/CAF )«CAPF .03 *Q8AaCSPF ) 
8«phipfv]2«-34S*CAPF403*(0RB-1  .OZCSF)  *CBPF) 


LM  JB/J  Lllr 


OII.Da  0 . S+OE *06LTI  /H- 
0(  t .2 )■  OX/H 
B(  1.2)  a-OO.2) 


0(1  .31*01 «02/H>  O.S*(OS«QELrA+OA) 
8(1 .4) a-Ol *G3/H+  ( O04OS*OSLT3) /2 .0 
0(  1 «*)»  0 1 *C3/H4(  09-0*  aQEL  TB  ) /2. 0 

Oi 


B(3.3>*-3.0 

O (3*3 1*4.0  

X(3.3>«-1.0 

8(4*4 )«— 3 .0 
0 (4 ,4 ) a4 .0 

JU  

CONTINUE 
CALL  BANO(J) 

80  TO  IS 


PROCRAM  FOR  J.6T.1  AMO  J.LT.NJ 


T.W4.7WI'’|I[j 


CAF  *O.Sa (C(3.J)+C(3»J4I  I 1 
C8F  aO.S*(C(4, J).C(4, JM) ) 
CURF«0.S«(C(1. J)4C( 1.J41) ) 


rjrij  XI  I 


CAPF  a(C( 3* J>1 )-C( 3. J) )/H 
CBPF  *(C(4.J*1 >-C(4.J ))/H 
PMIPF»(C(2. J*1 )-C(2 . J) )/M 
_CALL_3C0N0(CA(L»JC8F*.C0N0»£fli)0/ 
IF< TIME.ST .TMI ) 80  TO  42 
EPP*0.S«(EP(J>4EP( J4l  ) ) 
0l*C0N0*EPF4*( 1 .04T0T ) 


0(1.1) aO. S 
B(1  .2)«-01 /M 
0(1  .2) *01 /H 

BO.  31 « 1.0 

8(4 .4 )al .0 
6(3 ) *CA (NT 
6(4 )*CB 1NT 


42 

CONTINUE 

CALL  SACT I V (CAF.CBF, ACTA. ACTS .V A. YB .YAA • YAB • VBA. V8B ) 
CALL  S«ATR(TlR.T2R,CdA.CAF.CB«.PA.PB.PAA.PAB.PBA.PBB) 
8MFla(QM( J)»OM( jm ) 1/2.0 

OMNONFa  ( OMNON ( J 1 40  4S0N  ( J4  ll  ) /2^ 

OMFaOMFl-QMNONF 

0Fa(0(  J )4Q(  JM  D/2.0 

i«j  jpBKTjjpnm 


THIS  PIGS  IS  BEST  QUALITY  PRACXI£ABfj| 

rao*  oan  nmsHSD  to  ddo 


CALL  EQUI L ( CAP  « C9F  • E')Cs3  .E  Qvj  I Ad  « EOCBP) 

CALL  RK  AY  { DEL  T ,£PF  , O'AF  . QP  , AMP  . AMCUP  . AMCA . BMP.  RSP  * PMCUR • 

I RSCUR.BMCA.RSCA) 

_ DELTA *-OEL*HSP  _ 

OELTB*  DEL  *BSP*E0C3° 

01 »C0N0*6  PPa* (1*0 ♦TOT ) 

02*VA*PA/CAF 

_XL3P3fBPPB/CBP : 

0 AA*CONOA ♦ Y AA4PAA 
ORAaCONOA  4VBA4PSA 
0A8a  CON03  AY  A8 APAB 

_0BB*CON08  ♦YBB^PBB __ 

oe*  < 1.0»T0T)a0l*(»HtoPO2*CAPP>03*CBPP)/EPP 

OAa  01* (CONOAAPHt P«*0?»(TAA-l . O/CAP )*CAPP403 *08A *CSPP 1 

08a  0la(C0N0a*PHtPPO2*')AB«CAPP«-03A(0BB-l  .0 /C3P)  *C8PP » 

an  . t taa.p-cganpt  ri/n 

Oil •! laOiS ♦CE*OELT l/H 

011.2)*  01/M 

Bit  «21»-0(  1 *2J^ 

8(1 .3)*-Cl «Q2/M+0.SB~( CEaOELT A+O A 1 
0(1 .31*01 *02/M*  0.5A(OS«OELTA>OA) 

8(1  .*l*-Ot  *C3/MA.(OBaOS*OELTS  1/2.0 

__0( I . *JL»  O l *03/M4<  090E*0*L T3 1/2 . 0 

0(1  ) ■ OEaOELTI*CUR»P«>OEAOSLTAaCAP*(JAaCIPA.OB*CBPA-ae*OeLTS 
CUPP  «(C(  1 .J*l  >-C<  1 .J-l  1 1 /2.0/M 
PMIP  ■(C(2«JM)*C(2«J*1II/1.9/H 


VBLP  *(C(9.J4l 1-C( 9. J-l 1 1 /2.0/H 
ZVELP*(Z(9.J»1  )-Z(S.  J-l  11/Z.O/M 

PMHP*«C(6.  J*l  l-CI.4.  J — l ) > /2.0ZH 

~ZCUPP*( Z( 1 . J*l 1-Z( 1. J-l 1 > /2.0/M 
ZPM X P*CZ( 2. J*l )-Z(2.J-t 11/2. 0/H 
ZCAP  •( Z<  3 . J + l 1— Z ( 3 • J- 1 ) I /2. 0/M 

■X£BP,*jLZ.<_»«J»1  Ir.ll 

CLPPP*( C(  1 . J*1  )*C< 1 . J-l  >-2.0*C< I * J 1 l/M/M 
CAPP  *(C( 3. JM 1*C(3.  J-l l-2.0*C( 3.J1 l/M/M 
CBPP  *(C(4. J*1 1>C(4. J-l  I -2 .0  *C(  * » J ) l/M/M 

_ZCURPPa(Z(  1.  J41  lt,Z(.».A_Jr»  Lr2._0*ZU  « J-U/H/M 

^CCAPP  *(Z(3. JM »4Z( 3. J-t 1—2. 0*Z( 3. J) l/M/M 
ZC8PP  • ( Z ( 4 « JM >4Z(4. J-l 1-2.0 *Z (4. J) l/M/M 
BPPCaCEPI J ♦ 1 1 -EP( J— l 1 ) /2.0/M 

-|«PPC*(jEPUtIir.«P(j-t  

^Xd«c(  3.j) 

CBC*C(4.J I 
ZAC-ZO.J1 

_ZBC*zr4.j; 

VELCaCIS.  J) 

ZVELaZ(S.J) 

OCaO(J) 



CPCaEPI J) 

ZEPC»ZEP(  J 1 
OMCaQMCl-OMNONC 

"CALL  EOUIL  ( CACTCrCTBOCaVEOiTlXBTraCBin 
CALL  eau(L( ZAC.ZBC . ZSOCB. ZE0UA8 . ZEOCBPl 

CALL  RK AY (OELT.EPC.QPC.aC.AMC. AMCUP. AMCA.RMC.PSC. PMCUR. RSCUR. 


1 BMC  A.  R 3CJL1 

TBTta*-oel4 


n5lLTA-i-oeL*Tfl2 

OELT8«OEL4RSC4EOCB“ 

CMO*CSTAR— OSTAR 

APB* AST  A 4BSTAR  

— CALL-  SO  (PC  CAC".  C BCT^TF r.  BTFST5TFTS7BT FTHTToTPl  a . Cfl  PftcfJ 
call  so  ip ( zac.zbc. zo i PA.zotPB.ro tPA a. zoi*ab. zot p e a , zo i pee  1 

CALL  30  IPP( CAC.CBC. 0 IPAAA. 01  PAR A .01 PABB .0 (PBAA. 0 1 PBAB .Ot PBBB 1 
PPTaPPC** TOR 

~ SfePT.ZEPC*  B-fW 

TB2*T2R*CAC/C0A 
TR3*l • O-TR l -TP 2 

OlaOIPA*EPT4CAPP 

""02*  EPT4CAP4CAP 
03*  EPT4CAP4C8P 
04*  ( I .04T0R1*EPT*0IPA*CAP/EPC 

OS*  ( — TR2 ♦ S2  1 /2  .0/P/PPC-CAC49STAR/P.PC 

“Wat  SS2-tAC*0STApr*(  vPt>4LSB*"Cufcfr»7CMff71P^ 
Z01*Z0(PA*ZEPT4ZCAP9 

ZCa*ZEPT*ZCAP4ZCAP  • 


1*7 


jSlSPAflE 


IODDQ  — 


f 0*2 ( f *0*  TORI* ZEPT*  ZO t P A*  ZC  AP/ Z6PC 

ZOS*  < -TR2*  S2 ) /2.0/P/ZE»C-ZAC  *BSTAP/ZEPC 

>g*i  <ss2-ZAC«0STAR  ) •{  zv»UP*AP-MZOJRP)  /CMO/ZEPC  . ..  — 

0T*VELC*CAP/EPC*T0B*ai*'tV'3,V-l  .0  )*EPPC*QA*CURP*aS*  66 

t ♦TOP *0 2*0  I PA*0 IPA  A ♦TOa*QJ*OlPA*OIPAB 

I0Ca«^eLT8»OT/ePC*9lPAl*<Ql*ipPC*O4 1 

X*Q2*0lPABA*03*0tPA83  — 

lQVCY*ISS2*C«C*OiTAft)/CAO/!PC  " 

ZOVE Y*C  332-Z/  C *0  STAR ) /CM9/ZSPC 
AC37ll«-<QS*UVEV*AP9-vll«l/4.0/H 

g (J}_,  1 1 » ♦(  05 »OVgV* *PS~0 l ° J/A . 0 /H 


ai5.3>*  -o.s*otPA*ePT/H/H*as/*«o/H^o.s/DeuT-ocA/A.o 

0(3,3).  -0.3*0  IPA  •EPT/M/'i-OS/A.  0/M  *O.S/DELT-aCA/* .0 

B < 3 • 3> »0 I P A *ep T/H/M_ 

A<3.A)«-QCB/*.^) 

0<3.*)«-OCB/A.O 


ni3|4  ■•uv.o/  • • w 

A(3*s>-  -QVEV/A.o/HACAP/Sac/A.o 

^Us^r->XSiYT**0Cfltjf(-fe^*i3kL7gEPCizai7(  ia2*Z6t^AA7zay^ 

i^zoipabT^zo^paazkppcazq*  -zc'jpppzoa-zovevAe  zvelp*apb*zcurpi 

2-01 p*c URP-QCA  *CAC-0CB*C3C*C AP* V EUC/EPC  ) 

P la  OIPB*EPT*CBPP ■ 

~ P 2»  EPT  • CB  P PCBP 

«2*?Io*TOP^*ePT*OIPBPCP»/SPC 

-SlTO3^BCTgi  TAB  ) ITvIcuHpo » CU BP? /» CMU7 bW! 

l"S5&5555TwS.;;r.?si:«ToTi;;j;5Sc‘",'<:“ 

R XP»OELT|  AR7/EPC  ^ [raAiC-ll  

'zSl*  ZD IPS* ZEP T*ZCBP° 

ZP2»ZBPT*  ZCeP*ZC8P 

*a3.ZgPT»2CAP»ZCBP ..rj..L,j,Liir 

— ZP*-»  OTOVT  0 B ) * ZEP  T SZo  i p a*  4CB®  / z . 

IS:!;j;^KS5«S"?!;SSJSyi!5SS„c»/it»c 


A(4* I>»-(R9*«VEY*AP3-RIP»/A.0/M 
0(A • I )*♦( as ♦PVEY* AP3-R IPI/*«0/H 
A(4.3)»-BCA/A,0 

“BT*.'3yi^«CA/A.7 


M***l»-0«S*0tPB*E®T/H/H  ♦P3/A.0/H-RC8/A.0  lS*2yoiLT 
0<A.*J»-0.5*0IPB*EPT/H/H  -R3/A.0/H-RCB/4.0  ♦O.S/OELT 

~ SjgWCTPyTTTT*TqWT*gPy*P  |g3V»PPt/EPC 

p-Rv^^o%^s^:to°  

a^PXp2cuPP-BCAic“c-BCB*C9C*CBPPVEUC/EPC» 


CURB  ■O.SAlC(ltJJ*CM.J->»» 
PM1B  »0 ,3.<_C-Cl2_xail£f  ?■»  J-l 1 ?- 

CAB^*0 .5*  c Cl  3 . J I *C  J 3 . J-l  I » 
CBS  *0.S*ICC4.JI*CC4. J-l II 
VEL6  *0.S*ICIS.JI*C<3.J-I l{ 
_PHXlB*0.3*lCl6.J)*Ct.&.J-l U . 
CUPPS* (C( l . JI-CC l. J-l I l/H 
P«t IPB*IC(6. JI-CI6. J-t ))/H 
QMO 1*1 OP  I J I *OMC J— I l > /2.0 


QP*0  I QKII  9 tun  x u I 9 bp  s • * 

pMNONB.  tQMHONt— XLtOMNJN  LJj.U 
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^BIS  PAQB  IS  BB5T  QUALITY  PBACHCABLl 

raoi  oory  pumushd  xo  d do  


QPBaQMB 1 -0  MNONB 
QRa(Q<  J >#0<  J-l  n/2.0 

ep8«c cp ( j )#cp< j-t ) >/».a 

_cau.  eouiu(CAatCBa«c3ca*s.0'jXA8*CQCB»> 

CALL.  BK  AY  < OELT  .EPS  • 3**d  * 33  • AM8 . AMCUP  • AMC  A.  BM8  .BSB  .BMO 1 . BSD  I • 
t RMCA.RSCA) 

AMCa*«AMCA*eOCBB 

_RMCB9-BMCA*CQCa® 

b sca«-B  sc  a * eac  a® 

0BLTA#— 0€L#RS8 
0ELT8a0EL*RSS#SQC8® 

_ABCT(J)*AMe 

IBt TIMe.GT.TMl > SO  TO  35 

ACBN«OC  AP  ( A4.BM  A*(  PH  It  5-®H  1 9 > AB/B/T  ) 

CCPMaOCXP  < - ACPHC*  < PH  U B-PH I A J #B/B/T  » 

_gsBsAEPitacgai __ 

CXPO*<  AL®HAAABPNAA|.®HCACSPN)  ab/b/t 
G*aB  10  A ( CAB/CO A) ##CAM#<  C9S/C0S ) **ZETA 
Alt. 1 >a-l . 0/M 

“SHtilHie #G*  #ex®0/2.0 
aca.aia  AMBAC*  aexpo/2.0 

AC  2 ••) ■— AM8#G*  ACXPQ/2.0 

BI2.A1— AHHAC*  ACXPQ/2.0 

Cta)  AAMBAC*  ACCXfl-CXPOAPHI  IBACXPOAPHIB) 

BO  TO  *3 
CONTINUC 

BOMBaB«B»PCHO/BMO 

BQMCUBaAMOl 

BOMCAaBMCA 

BOHCBaBMCS 

HlfM.U.fMMI  1 fi«.f  At  Ml.  ■ - 

IPfCSB.LB.CALML)  C93»CALML 
ASAKal .0-S2ACUBBB/SH/B/CAS/BMS 
BSBKat • 0“ S3  ACUBPB /SN/B /C53/ 90H8 

__ASAX  MZiS2/ CM/P/CAB  /R*3 

BSBKM2  *S3/CN/B /C6B/93H8 
ASAKMl aASAKP2#CUBP9 
BS8KM|aBSaKM2*CUBPB 


IPf BSBK.Le.CALML)  339KaCALML 
I BfASAiC.Ce  .CAL  MM)  AS  AKaCAL  AH 

IB< BSBK.CC .CALMH ) 9S3KaCALHH 

CAS ( J )■ ASAK 
CBATtOf  J ) a ASAK /BS8K 
CBAT 10 ( t ) aCRA T IOC  2 ) 

CAAVafASAK *_CA 8*  C A 8 )/2  .0 

CBAVaf  BS8K*C8B+CB8  ) /2 .0 

CALL  SACTI VCCAAV .C8AV . ACTA . ACTS. YA. YB . Y AA  * V AS. Y8 A . Y88 i 
TB2aT2B*CA8/C0A 

ia?iay:a7r3?tBi»aYA 

ZZCTBa-C* .0-2 .0*TR3 I *V9 
IBf JP.CQ.2)  GO  TO  «S 

HH&tSWfH&Hr 

IBfJP.EO.I)  GO  TO  *S 
ZZCTAaO 

- *s  — ISWTOa 

ZCTAAaCAM  —ALPHA#  Z7ET  A 
ZBTBAaZETA— ALPHA#  ZZ5TB 

feTAC»CAM  #ALPHC#ZZSTA 

eTBCaZCT i# ALPhC*  Z ZST  8 


AeTBCaZ6Ti#ALPMC*ZZST8 
CAAM t aA SAK#*( ZC TAA-1.0) 
CBAMlaBS8K##(ZeT8A-t.O) 

COHA AaASAK# CAAM t 
C0N8AaBSBK*CBAMl 

C0NAC#ASAK«CACM1 

~CCN8CaBSBK#CBCMl 

C J^ONAaOexPC  ALPHA #B#  C PHI l 9— ®H  tB ) /B/TI 
exPONC  aOCXPf -ALPHC  #p#  C PH I 1 8-PH I B ) /B /T ) 

COHM l*( C8B/C0B )*•< ZSTA-I .0) 

CZNaCZNMl #CAB/COA 
CO  Ha  COHM t #CBB/COB 
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PAGES  IS  BEST  QUALITY  PRACTICABL* 




PREXPABAM84FI0 4CZN4C3HNI  •-  

PaexPSaAMB«FtO*CGH«CZNMl  •"* 

PREXPO* AMB*PCO*CZN4COM 
PANG0*PREXPG*CAAM1 «C?AMl«*XPQNA 

P~CA'TD*PREXPO*CACMt  •C3CVl>?*Pa*R; 

8PH!bPREXP0*( ALPHA*C3NA4*C3NQ4*EXPONA 
l ♦ALPMC4CONAC*CON8C*EX®ON  C ) *P/R/T 
BAN  t *( PAN004ZE  TAA-aCATOAZST AC }«BSBK 

~BCTI«<PANOO*ZETBA-PCATO»zeTHC>AASAK 

Cll»P«EXPO*EXPa*AMCUP/ AM* HAN l «ASAKM2«( -1 .0 ♦CURPB4RM0 I ) 

1 ♦flCTl *8SBKM2*( —I « 0 4CURP iaROHCUR) 

^ 0A2aPHEX®0»gXP0»AMCA/4X8»ORgxaB4EXP0«GAM/C0A 

T+bani  •aSakpi  • ( i .o/Ci¥+-RNC4>  VhctI  •SSBKmTTSom^a 
683*PBEXPO*EXPO*AMC9/AMR*PReXPA4EXPO*Z*TA/COS 
I 48ANl«ASAK*l«RMCa  -MCT  t*9SB<N  1 a(  1 .0/CBB>R0HCB  I 
A(2,l)*-1 .O/H+GIl/H 

~B(a; I )«=A<  27l)  

A(2.2)«8PHt/2.0 
8(2* 21 *BPMI /2. 0 
3t*-GA2/2 


'2 

AC2.4)a-G83/2.0 

8(  2 »A  ) *—683/2  ._0 

A(2,6)«-8PMt/2.0 
8(2.6 ta-SPHI/2.0 

G(  2 ) »PR  EXPQaEXPO-BPH t *PHI l B+8PH laPHIB— G1 1 4CURP8— 6A2 4C48— GB34C8B 

continue  

1 1 iRsaa  < CSTAft-fTST  art 
T2*Tt*EQUlAB 
TR2»T2R  *C A8/C0A 
TR3  « 1 • 0 — T R 


CALL  SIGMA(EPGI.EP3.QM8.a8.SGZN.SGZN0.SIGMt ) 
IPCTIME.GT.TMI ) GO  TO  47 

ACS.lt*  -<ASTA  *8  STAR  I /H 

“8(5. I J*-A(S« H 
A(3.5)*-l  • Q/H 
8(S.S)*1.0/H 


8C6.lt*  -0.8 

ACS. 6)*  - ( I .0— EPS ) 4S I GM I/H 

8 ( A . 6 ) a -AC6.6) 

"T.  < 6 ) a'  -rrrtOT 
60  TO  49 
CONTINUE 

ACS. 11*  -CASTA  ♦BSTARI/H  » T2*RSOt/M 


UflTr7UW1 


A(5.3)a— O.SaTl — 0. S*T2*R  SC A 
8(8.3)*  ACS. 3) 

AC 8 . 41 *0.3 • TI * EOCBP— 0 «34T2»RSCB 

t(S.4)a  A(S.4) 

A CS .8  t«— l .0/M 
8(9.8>a| .0/M 

«T2-T l • CAB AT 1 •1TS8®*C83-T2* ( R SC A*C * 8 ARSC8*CBB 


BC6.lt  a-0  . S+SGZN4  ( -OSL  TI  ) *®H  1 1 P8*  ( 1 .0-SP8  )/H-S  IGM  I *OH  UP8405L  T I/H 
AC6 , 3)a0.8*SGZNO*C  —OSL  TA ) *Prt I 1P84C I .0-EP8 >-0.S*S I GM 1 4PM I IP840ELTA 
BieOta  AC6.3I  

BC6*4)a  AC6.4) 

AC6.6)a-SIGMI*( 1 .0-E®8)/H 

8(6.6)*  3 C GMI *( t • 0-E®3 ) /M  


1 


PHI 15a0.SP(C(6.J)*C(6. J-ll  1 
CURPMICI  t.  Jl-Ct  t . J-l  > >/N 
PHllPBa(CC6.J)-Ct6. J-l >>/H 
~OS»(Q<  J )*Q( J-t I 1/2.0 

epe«<ep< j »*ep(  j-i  n/2.o 
o**ei«(OM<  j)*op<  j-i  i >/2.o 

_QMNCNBa(QMN0N{ J 1 TOMNON  {. J-l  ) 1/2.0 

QMBaQMBl -GMNONB 

CALC  EQUlLt CAB. 088.8008. EQJIAB.EOCflPJ 

call  RXAVCOELT.EPB.QMa.aB.AMd.AMCUR.AMCA.RMB.RSB.RMOl .RSOI. 

J RMCA^R  SCAJ. 

AMCB«-AMCAPE0C9P 
RMC0a-RPC  ApEOCBP 
R5CB»— R SC APEOCBP 

_DOLT A«-OCLPRSB  

OEL  T RaOEL  «R SBPfOCBP 
AMET(  j)«am*» 

IPiTiaa.^r.rai > oc  rn  36 

LEBHaOexPt  Ali>t1  A«(  EHlId-rgHU  1 PfVB/T  1 . 

CEPNaOE XP ( -ACPHCa < PH  1 1 3 -oh I S > PP/H/T > 

EXRaASPN— C8PN 

EXPt)«(  ALPHAPAEPNPALPHCpCEPN  ) PP/P/T 

C4-P  10  P ( C AS /CQA  > PPG  A*iL£MZC  08 1PPZETA 

A<2. I >•— l .0/H 

8(2. 1 1 a t .0/H 

A(2.2)a  AMS *G*  PEXPQ/2.0 

8(2. 2)v  AMB *04 E6XED/ 2j0 

At 2 »6)a— AN8*G4  PEXPO/2.0 
8(2 • A) a— AMB*<i*  PEXPO/2.0 

G(2)aAMSPG4  P<EXR-*XPOPPH118PEXPOPPHI8» 

SO_ TO-  60 

36  CONTINUE 

R CHB«RM9P« CHO/ RMO 
ROHCURaRROI 

BOUCApRMCA 

ROHCSaRMCB 

I PC  CA8.LE.CALML)  CASsCALML 
IP( C88.LE .CALMC ) CMBaCALML 

ASAKal .0-S2PCURP9/EN/P /CAS/RMS 

BSBXal .0-S3*CUPPa/EN/P/C33/R0MB 
I PC  A SAX «LE. CAL ML)  ASAX*CALML 
tP(  BSSX  «L  E .CAL ML  ) 8S8XaCALML 

IFCBSBX. gg.CALMM)  ASAKagALMH 

IPC ASAX .GE.CALMH)  ASAXaCALMH 
CASC J)*ASAK 
CRATTOt J) »A  SAK/8S8K 

ASAXM2aS2/EN/PyCAB/RMJ 

BSBKM2aS3/EN/P/C8B/R'3HS 
A SAKM 1 a A SAKM2PCUPP8 
BSBK  Mia  8SBKM2  pCURP  6 

CAAVa( A«AXPCAB»CABl/3.a 

CBAVaC  BSBK  PCBS+CB8 1/2.0 

CALL  SACT  I V ( CA AV • C3 AV * ACT  A. ACTS « YA . YB . Y AA  « Y AB . Y8 A • VBB I 
TR2aT2R PCAB/COA 

IR3a l .OtT IR-TR  2 

ZZETAal  .SM1.0+TP2  >«TA 
ZZETBa- (A. 0-2. OPTR  3 ) PY9 
IP( JP.E0.2 ) GO  TO  59 

TZPTAml.ap (1.0ATR21 

' Z2ET8*  — ( A • C-2. Op  TP  3 1 
! P( JP.EO. I ) GO  TO  55 

ZZETAaO  

ZZETB-O 
85  CONTINUE 

ZETAAaGAM  -ALPHAPZZETA 

ZETBAaZETA— alpmapzzetm • 

ZCTACaGAM  PALPHCPZZETA 
ZETBCaZETAPALPHCPZZETB 
C AAMI a A SAX PP ( 2ETAA— l .0  ) 

«aaSKKHH*Hh« 

CSC Ml «BS8XPP ( ZSTBC— l .0  ) 

CON  A AaA  SAXPCA  A Ml 

CONBAaBSBXPCBAMl  

CONACaASAXPCACMt 
CONPCaesaxPCBCMi 

EX—ONAaOEXP<  ALPHA PPP( PHI I 9— PH  IB l/R/T ) 

. .EXPONC  *OE  XP( -ALPHCPPP ( PHI  1 J-PHlB l/R/T ) 


I / • . 


151 


Oil  Mg  18  BBT  QUALITY  PlACTlCABLi 
IM  OQEX  IWBttSHBD  XO  DDQ - 


EXPOaCONAAACONBAAEXPONA-CONACACONBCAEXPONC  ’ 

CZNMIaf  CAB/COA)«*( GAM-l .0 ) 

COHM1*: C8Q/C0S )*• ( 2ET4-1 .0) 

CZNaCZNMl  AC AB/COA 

coH«coHMi«caa/coe 

POPXPAaAMBAPIOACZNACOHMl 
PRC XPfla AMS aPtOaCCHaCZ^Mt 

PRexPOaAP8aPIOaCZN#Cia 

PANOOaPRE XPO«C AAM I aC  JAM 1 aSXPONA 
PCATOaPREXPOaCACM  t aC‘«CM l aEXPONC 
BPMl aPREXPOAl ALPH AMC0N4 AACONHAMEXPONA 

J ♦ALPHCaC0NACaC0NBCa*xP0NC>ap/H/T  

8AN1 a(PANQOaZETAA-PCATOaZETAC IABS9K'  ~ 
■CT|aCPANQ0aZFT8A-«CAT0aZ»TaC»aASAK 

fi!laPREXPQaEXP0aAMCUR/AMja8ANiaASAKM2a(-| .OaCUPPBARMOI ) 

I 48CTlaBSBKM2a( -1  .0aC<JQP3a40hCIiR) 

"GAZaPR'ExaO  aSXPO'aa  MCA/VH8'aP'R*x66aEX?n5»GAM/,io  A 
1 4BANI  a AS  AKX 1 a ( t • 0 /C 4 3 ARM  CA  ) A8CT l A8  SflKN  t aRCHCA 
6B3«PRSXPQASXPOAAMCB/AM0APR«xPAAEXPOAZST A/COB 
1 ABAN1 AASAKM14RMCB  ♦8CTl*BSaKMXa( 1 . O/CSBAROHCSI 

“Jk(2.l)a-i  .0/MaGI1/R 

8(2.1 )a-A(2,t) 

A(2*  2) aBPHI /2.0 


ftp  CCNT1NU6 
TiaRSBad 


B(2.J>a-GA2/2.0 
A(2 .4 ) a— G83/2. 0 

8 ( 2 . ♦ ) a— G3  3/2 • 0 

A < 2 . 6 ) a -0P h I /2V0“ 

B<  2 . 6) a— 8 PM  I /2 .0 

G(2  > aPREXPQAEXPO— 8PH I aPHI IBabph I a pm  1 8— G 1 1 4CURP8— GA2aCA8— G83ACBB 


T l a R SB a < C TTaR-0STA9) 

T 2a T 1 4EQU t A8 
TR2*T2R AC AB/COA 
TR3*1 • 0— TR2— TR 1 

ZCURPBa  ( Z ( lVj  > -Z<  T7J=’TTJ'7fl 

CALL  S I CM A ( EPG I . EPB. QMS. QB. SGZN . SGZNO . S IGMT  ) 
IP(TH*E.GT.TMI  1 GO  TO  S3 

8 ( 3.3 1*1.9 

i"> 

iC 
>C 


G(6 I*  -PITOT 
GO  TO  S2 

CONTINUE ^ 

A(5  < I )*  - ( ASYA  aeSYTOl/H  a t2*R SOT/M 
B(S*1  1 ■— A (3.1) 

A (3 . 3) ■— 0. S*T1 —0. 34T24RSCA 

s. 

a' 

8(3. *)•  A ( 3. A) 

A(S*S)a— 1 «0/M 

G ( 3 > «T2-T  1 aCABVTl  aEo^P^CMiT Za  ( R'SC %•£ a'Ba  »SC  Ba CiBY- a*fS5Ta CCi Rpft _ 
A(A . I ) a— 4 . 3-SGZNa ( -02L T I )*SM1 !P6a( l.O-EPB I/hasIGm I aPHI 1P0AOELTI/M 
8(6.1) ■— 0 . SaSGZNA ( -OEL  TI »a»MllPSa(l  «0-EPR  >/ H-MGM  I APHI  IPB  A05LT  l/N 
A(6.3)ap.3a 


ULi X U'n-Lt  1 L'Ulli-. I * tliUL’J 


A ( 6. A)aO. 3*SGZN0a ( -OEL  T3I **HI IPB* ( 1 .0-CP8 )-0.3aS IGN I aPHI IP0AOELTB 
8(6. Ala  A(6.«> 

A(6.6)a-SIGPla(l .0-EPB1/M  

~8 (6 . 6) ■ StGNlat l.O-aai)  A( 

G(6 ) a -PI TOTaPHl lPja(SGZMA< -OELT I ' ACU*®8  aSGZNOal -OBLTA) aCAS 
1 aSGZNOai -0ELT8 )aC  48) « ( l .O-EPOI-S ICRtaPMl |P8a(0SLTiaCURP8aoSLTAa 

^calc^sotptc  ( 5?  3 TTCTrrTTTTPTmx  r»mnpjra~inPZffTT»  r<nrAT7TPtf>n 

CALL  SO IP( Z(  3. J) .Z(4. Jl .ZOlPA.ZOIPB.ZDIPAA.ZDIPAS.ZOtPBA. Z0IPB9) 
CALL  EOUlL( C(3.J).C(A.JI.*3C8 .EOU t AB .EOCBP) 

OMPaOMI j| -OMNON( J ) 


Its  B*6* 


****** 


“ CAL hi  RKAr(0ELT.EP(  3)  .OMMJQI  J)  . AM  , AMCUK  . AMCA  .RM,  RSI  J J . RMCUB.'RSCUR 

t .RMCA.RSCA) 

OELT  A«-QEL*RS( J ) 

DEL  TB=o£l  *R  S<_J  > *6008? 

CAP* (3 .U*C(3. J)-* .0*C( 3. J-t J >CC 3. J-2) )/2. fi/H  “~ 

C8P*( 3.0*C( 4, j )-4.0«C (A , J-t ) +C(4. J -2) )/2.0/m 
ZCAP*( 3.0*Z( 3. J)-4.0*Z( 1. J-l )+Z(3.J-2> ) /2.0/H 
ZC0P*<3  .0*Z(4.  J ) -4. 0*2 (4.  j-n  ( 4.  J-2) > /2.0/H 

“ CURP*(  3.0«C(  l . J)-4  .J*C(  IV  J-l ) *C(  1 .J -2 ) ) /2.0/H 
EPT*EP( J)**TOR 
ZEPT»ZEP<  J )**TOR 
04*( 1.0>TORJ*EPT*OIPA*C4P 

^j4-rr;o  ♦tor  r«e^t*oTP^*t^s 

IFINBC.eQ.O)  GO  TO  82 
C3NJ«C(  3.NJ) 

C4NJ*C(4.NJ) 

~ IP(L.EO.l)  VELMEM.vSCT 

IF(L.E0.2>  VELMEM»-V*I_M 
V ELO IP* V£LMEM-C( 5 . J ) 

VELZIP*VELMSM-Z(3, J1 

^TRN2  BIT  R 2 M -TR2 

TRN30»TR3M-TR3 
IP(NBC.EQ.l)  GO  TO  83 
IPIN8C.EQ.2)  TRN2D*  0.0 
IP ( NBC. GE. 3 ) TRN3 0*0  .0 
I PI NBC. EQ .4 ) TRN20*0.0 

_GQ  TO.  83 • 

CONTINUE 
C3NJ*0.0 
C4N J*0 .0 


VELZ  IP*0.0 
TRN2 0*0 .0 
TRN 30*0.0 

83 CONTINUE. 

VI3.1I*  04* OELT  t/2 .0/H 
A (3.  I )*-04*0SL TI *4. 0/2. 0/H 
8(3.1)*  04 *OELT I *3 .0/2. 0/H 


8(3.3>»-C3NJ 

Y (3 .3 )*  E* ( J )*E®T  *0  IP  4/ 2.0/H 
A ( 3 • 3) *— EP< J ) *EPT*0 I P 4*4 . 0/2 . 0/H 

a<3  .3  l.»__EP  l J ) • EPT *0  IF 4 *.3 .0/2.  0/H+2..0*  VOLUPE/DELItGLAA  MiYSLOXf 

I *EP( J )*EPT*CA»*0IF4*0  IFAA404»0ELTA 

6(31 *2. 0*VCLUME*  Z ( 3. J ) /OELT+2 ,0*TRN20*F ITCT /2 .0/ F+2 .0*OLAAM*C0 A 
1 -ZOIPA*ZEP( J )*ZEPT*ZCAP  -(DLAAM.-VELZIF  )*Z(3.J) 


_ C3NJ*C(  3.  J) 

VIA. 11*  R4*0«LTI/2.0/H 
A I 4 • 1 ) »— 4. 0*R4*OELTI /2.0/H 

_B(4.l)  ■ 3.0*R4*OELTI/2.Q/H 

8(4.3 )*EP<  J ) *EPT*C3P*0 IF8*0IPBA+R4*0ELTA 
8(4*3>*-C4KJ 

r(4.4>*  EP( J )*EPT*0tP8/2.0/H 


8(4.41*  EP( J)*»PT»0 1*8*3. 0/2. 0/H*2 . 0*V0LUME/0ELT*DL8BN+VEL0IF 
1 *EP( J) *EPT*C3»*0I F340IFB8  ♦R4*0ELTB 

G(41 *2.0* VOLUME *Z(4 .J l/OELT  *2. 0*OLBBM* COB-ZEP ( J 1 *ZEPT*ZO IF3*ZCBP 

1 -(OLBBM+VELZIP >*Z(4.J>  v(P( J > *EPT*CBP»Q1FB» (01 FB4*C ( 3 » J ) » 

2~  01 FS8* C(  4 . J ) ) *^4  * ( OEL  T I *CURP  ♦DELTA*C(r3.j’)  ♦OEL  T8*C ( 4 . J > ) 

3 -C*NJ*C(5.J>  ♦2.0*TRN30*PIT0T/P 

CONTINUE 
ALI 3ANQCJJ 


rm 


00  *76  J*2. NJM 1 

476  TRNCuR ( J)*(C(l .J4l )-C(l .J-l 1 )/( A0*2.0*H) 

TRNCUR( I) *(-C( 1 .3)**<0*C(1 .21-3.0  *C ( 1 .1 ) ) /( A0*2.0*H) 

TRNCUH(NJl«(C( 1.NJM21— 4.0*C( 1 .NJM1 1 43.0*C ( 1 .NJ) 1 /( A0*2.0*H) 

I P ( JCOUNT . GT.2 0 ) GO  TO  39 
DO  *79  t* I «N 

OQ.  479  J*2..NJMJ 

|P(OAOS(C(  I .J) 1.LT. 1.00-10)  GO  TO  *79 

1 P(  0*tlS(  C(  I . J 1 -C I ( I.J)).GT.t.  00— 3*0  ABS  (C(  t.J))1  GO  TO  8 
479  CONTINUE 

GO  TO  39 

478  CONTINUE 

00  71  J*t.NJ 

CALL  EOUt L ( C(3 . J) .C (* • J 1 . EOCB .SOUL ( J) .EOCBP ) 

_71 CONTINUE 

00  72  J*1 .NJ 
OMN*OM(  J|-OMNON(  J) 

CALL  EOUlL(C(3.J).C(4.  J >.  E.1C8  .EQMt  AS. E0C0P) 


aaaBrssrr" 


'CALL  PKAY ( OELT |EP( J I I OMM. Q ( J ) . AM • AM  CUR . AMCA . RM, RS ( J >. RMCUR.RS CUP . 

1 RMC  A. RSC  A ) 

A MET ( J ) *AM 

CONTINUE  

OOEL  I- A ST  A*OEL  T/2 . 0 
00EL*CSTAR*0ELT/2 .0 
00  70  J*2.NJMI 

curp  *<c<  i . j»i  )— c n,j-n)/;.o/H 

— *Cu«<P*<  Z ( lTj*l)~Z  ( l.  J-t  l 1 /2.0/H 
CALL  EQUIL ( Z ( 3 • J ) ,Z(4.J).ZE0C*S. ZEOU A8 . ZEOCBP ) 

OM(  j)<Z0M(  J l-OOEL  I * ( C'JRPTZCURP  ) 

Q< JI*ZQ( J)*OOEL*(RS( J)*«0UL(  J ) *ZQS  < Jl*ZEQUAB) 

~ EP(  J )*ZEP(  J )*00EL  t *(  CURP*ZCUQp  P-dOEL*<  R S<  J)  *EOUL ( J ) *ZRS< J)*ZEQUAB  ) 
CONTINUE 

CURF*(-C< 1 .3l*«.0*Cd .2)-3.0*C( I . 1 ) 1/2.0/H 

2CUPP»<—  Z(  1..  3 ) **.0»Z(  1.2)-3.q*Z<  1.  1))/2.Q/H 

CALL  EOUILI Z(3.1).Z(A,l). ZEoTsTTEQUAB, ZEOCBP) 

0M(  I )*ZQM(  I ) -OOEL  I *<  CURF*ZCURB) 

0(1  IsZQd  ) *Q0EL  *<  RS  ( 1 >*E3JL(  1 )*ZRS<  1 )*ZEQUA8  ) 

EP( 1 )*ZEP( l I ♦QOEL I * ( CURF*  ZCURF  > — OOEL* t PS  < 1 )*EQULd  ) *ZRS ( l )*ZEOU AB ) 
CURB* ( C ( 1 ;NJM2>-A.0«C( 1 .NJMi )*3.0*C(1 .NJ) )/2.0/H 
ZCUR8* ( Z ( 1 .NJM2J-4.0 *Z< l • N JM 1 )*3.0*Z<  l . N J )) /2.0/M 
CALL  EOUtL( Z( 3 «NJ ) ,Z(4 .NJ) . ZE CCS . ZEOU AB . ZEOCBP) 

__0M(N  J )*ZOM< NJ ) -OPEL I»< CUR 3*ZCURS)  

&{NJ)*ZaCNJ)*a6EL*<RS(NJ>*=0UL(NJ>*ZRS<  NJ  )*ZEOUAB> 
gp< N J ) *ZEP ( NJ ) *QOEL I* ( CURB* ZCURB) -OOEL* ( R S( NJ ) *E  OUL (N J ) ♦ ZRS ( N J ) * 
t ZEOUAS) 

_00  73  J*1  .NJ _ 

IP(OM< J l.LE.O.O)  QmTj)*0.0 
IP(0(J).LE«0.0 ) 0(J)*0.0 
EP<  Jl*l.O-OM( J>-0< JI-EPGI 

-vm go  T ra 

IF( JPLG.EQ. X ) GO  TO  33 
MEP*0 

00  92  J»l.NJ 

— IFroC'J  )*.G£  .ESPCGlMEPaWSP^T 
[F(MEP.EQ.l)  N JM IN* J 
CONTINUE 

~P- 


00  97  J*!.NJ 

QMMAXaOMI J) *Q( J)*VM/VS  _ 

6mplg*6V:  j j -sp ( j ) *vm/( vs-vm  > 

IF( OMPLG.LE.O.O)  0MPLG*0.0 
OMCCV«OJS*(  QMMAX-QMOLG) 

QMNCN< J )*0MPLG*0MC0V 

OMM*OM<  J) -0*N0N( J ) 

CALL  EOUI L ( C(  3 . J ) .C(4.JI  .EOCB.EOUH  Jl  .EOCBP) 

CALL  RKAY  (OELT  .EP<  J)  . OMM. 0(J  ) . AM  ■ am  CUR  » AMCA  . RM.  R S(  J ) < RMCUR  .RSCUR  . 

1 RMC  A . R SC A I 

JTmET  ( J ) »A  M 

97  CONTINUE 

00  ABl  J*2 • N JM 1 

AREA(J)»A0*((I.0-EP(JI)/(t.0-EPt))**<2. 0/3.0)  

■“*  TRNCUR ( J)*(C(1 «J*l )-C( 1 .J-I ) ) /< AREA ( J ) *2 .0*M ) 

ABl  CONTINUE 

AflEAd  ) ■ A 0*  ( ( 1 .0-£P(l  ) )/(  1 .O-EPI  ) )•*(  2.0/3.  01 

AREA(NJ)  *A0*(  ( l .O-EP(NJ)  l/Jl  . O-FPM.  >••_<  3..0/3.0.) 

YRNCUR(  1 I ■ ( -C(  1.3)  *4.  )«C<  l .2)-3.0*C(l  . 1 ) I /(  AT»EATnV3T.7*M) 

TPNCURI NJ) *(C(  1 .NJM* ) -A .0*C( 1 .N JM 1 ) *3 .0 *C  d .N J I ) /( AREA( NJ )«2. 0*M) 
IF I JCOUNT • LT.  I 3)  GO  TO  S4 

WRITE! A. 202 1 . 


5w«w  — 


94 


NOT  CONVERGE* /// ) 


20X • ' 

CO  TO  59 
00  59  1*1. N 
oc  sa  J *2  . N JMl 

c^g;;fig^rjinrr7rrj;-j  * u h » 


•»*•«<«  i.  j>>  i oo  ro  9 

00  78  J»1 ,nj 


78 


.Arm  . 


CALt^SArrr  —9.  rA/,ACTB0*** ) 


wm  EM*F|  TCT/CONmem  

— gjjpliliiil gr— - 

74  SS?y?Nl J1»Q<J>/IVS*{09( J)/VM  


00  75  J-I.Nj  — 

ISSSfl'. ?>■  ’./.'iSTiTTi Ji'»l«T,jia»«.i  ,. 


8*  gSTKie™'7-7-'^^ 

SUM! *0.0 
SUM  2*0.  0. 


SUM 3-0.0 
SUM*— 0. 0 
SUM  5-0. 0 

_SUMAaO.fi 

00  88  j-t , n j 
SUM1-SUM1*CP(  j) 
SUM2-SU-2 ♦ QM  { j j 
Ji 


SUM*- sum*  *c  ( 3,'Jj 

f^s^uMs^cr*,  j| 

-88 cmx^u?6^0NVRN<  J' 


0 7N J -OF UOatTnTT 

SUMt-SUMl /OFNJ 

SUM2-SUM2/0FNJ 

-2aK:S'!S'22JH- 


'vi-nj. 

SUM*«SUM*/OFNJ 
SUM5-SUMS/0FMJ 
SUM6-SUM6/0FNJ 

-87 FORMAT  ( / • EP- • .Oil  » .*  SUMS  • SUM» 

,iF(Tisi:oT°Tu?^,3^isNh,,°,i-7'»7,2x,'0*i,ol6*'7,2J‘’,C4-^’7- 

^POT*gS<C<2,N-,,>*Cr»;,3rit">  00  TO  990 


81 


82 


00  61  J-I.NJ 

-♦a??a.-vfci»  

•MINT  62,  TM 

•i!!ig|^p«gQ«oi-*°oli3L43ION  4T  Tlwg,,»<rao»io//> 

CO  TO  Qo—  — L"t  — ■■  . 


CO  TO  998 
•3  CONTINUE 

- Df£TI?oeLT*‘TUIH*T,T^P,,  00  T0  ««0 


153 


flftl  MOB  u BBST  qUALin  FB4iG91SiM 


•V’  .’UK 

c ij-srt 


juuu 


NTMtaNT-1 

IF1NTI.LT.NTM1 | DeLTaT0EL/2.0**!NT-NTl 1 
IF1NT1  .EQ.NTMl  ) DELTaTOEL“T t ME 

NTI»NTI_*1 

t^(TIMB»OT.(TL I M— TDEL) ) GO  TO  999 
TIt*E»  TIME  ♦CELT 
GO  TO  2 

-999  - CONTINUE  

IFINCYCLM.EC.O)  GO  TO  990 
00  129  Jal.NJ 
Ctl. J)>0.0 

cxa^ii  «5UM9 

C I *» J )• SUMS 
129  CONTINUE 

MIT  0T«0 . 0 

IPITIMS.GT.ITLIM+TSTOO) ) GO  TO  988 
TIMEaTI ME+OELT 
GO  TO  2 

.909 CONTINUE  

IFINCYC.GE.NCYCLM)  G^TO  990 

LaC*l 

IFII.E0.3)  L*t 

— * 

DELTaTOEL/2  .04*! NT— NT I ) 

IFINCYCLM.EO.O)  GO  TO  990 

NCYCbNCYCM  __  

IMIL.EQ.t)  FITO'TiCtTOf 
I F! L.EQ .2 ) FITOTa-CITOT 
IF1L.E0.11  PRINT  132. TIME 

IfA  L-.EQ  .2  »_«*S.IJLT_J  2 \jLLU*S— 

131  FORM  AT  ( // 1 H »20X.  • THIS  IS  »±0R  CATHODIC  STARTING  TI  M£a*  ,F  I 0.3/ /) 

132  FORMAT! //IN  ,20X.«THIS  IS  FOR  ANOOIC  STARTING  TIME**  .F10.5//1 
IF(L  »EQ  .2  I GO  TO  6 

GO  TO  A 

990  CONTINUE 

IF1LPUN.EQ.0)  GO  TO  998 
00  99  J»1 .NJ.LPUN 

*"  S js<  t,.i  wj_cui/»jj;ut^i.±apt4)  ; 

WRlfEC7.l 03>TC6. J ).*P< J > .OM<j) ,Q<  j > ,RS<  J > 

99  CONTINUE 

WRITE!  7.103  HOWNONlJ). J» I. NJ.LPUN) 

WRITE  17.1031! EOUL (J).J*1. N J ,LPUN  ) 

WRITE!  7. 204)  TIME.OELTiTLtM.NTI' 

WRI TE! 7. 103)  SUM4. SUMS 
998  STOP 

ENO 


SUBROUTINE  8AN0!J) 

— IMPLICIT  «E*L»9 ( 

O I ME  NS  1 ON  A! 6. 6) .8! 6. 6). CCS. 93) .016 .1 3 ) .G!6 1 .X 16 .6 1 .Y 16 . 6 > < 
I E! 6.7, 93 ) 


gOMMON/aORKA/A.B.O.S.X.Y 

COMMON/WORK8/N . N J 

FORMAT! ISHOOSTERMaO  AT  Ja.141 

• IF1J-2I  1.6.8 

NP1 aN»l * 

00  2 I *1 «N 
Oil. 2 *N At )aG! I ) 

PO.-.E.  L-J  . N 

LPN  ■ LAN 

O! I »LPN ) a X! I ,L) 

CALL  MAT I NV ! N, 2*N  a 1 . 08 TERM  I 
IF! DETERM  14,3.4 
WRITE  16,101)  J 
00  S K«1,N 

E1K.NP1 .1 )aO!K.2*N*l  1 
00  S L«l.N 

Eiif  .c.  ni-TfnrTcr 

LPNaLAN 

X|K.L)a-0! K.LPN) 

RETURN 
00  7 tal.N 
00  7 Ka 1 . N 

00  7 L»I.N  


QlJJJjlW  * 
TODDC 


00  17  K*1  .M 

0(!.K)*0( l.K)-E*0< JCOL.X) 

CONTINUE 

.RETURN 

ENO 

SU8R0UT INS  SCON01 XA. X3 .CNO.CNOA .CNOB) 
IMPLICIT  RCAL*8< A-M.0-21 

_cpmmqn/_wqs  m/aggaso 

CNOAaQ.O 

iPixa.oe. 0.009)  go  to  3 

IP! XS.GE. 0.0037)  GO  TO  2 

IP(XB.Gg.O.OOlS)  GO  TO  1 

B 1*199.3362 

CNO*0.036S5+S1 *<X8-I .3 1 3* 10. 0 *• < -4 ) ) 
CN0B*81/CN0 

CO  TO  4 

81*1 34.09  73 

ONO*0. 300*  ai*(X8-0.00l3) 

CNOB  *8 1 /C  NO 

_GC.  TO  .4 

8 1*3. 141 39* (XB-O.  004*3  J/0.01 03/2.0 

CNO*0 .6372*0003 1 8 1 1 

CN08*— 3 1 .293*3 . 1 A 139  AOS  INI  l31  l/CNO 

-GO  IQ  .A 

CNO* 0.6 1—37 .BOB* ( X9-0  *000 ) 

CMOS*— 37.803/CNO 
CONTtNUC 

_CNO»CNQ*EEECNQ 

RETURN 

CHO 

SUB R OUT  INC  EQU I L ( C A . CS . S08 . CO «8 . EOP ) 

JL.mpjl1CJ.T_9 tonal  a-h.q-z) 

COMMON/WORKH/TOT.L 
IP'CB.GE. 0.009)  GO  TO  2 
IFICB.GE. 0.002)  GO  TO  I 
E 08*0.0 
E0A8*CA— COB 
COP-0. 0 

_eOB^2?0-t 

IP< EQ8.LT. 0.0)  eoa*o.o 
• E0A8-CA— €08 

Sfr  

EOB*l.a*t  0.0** (-4 )*0.0  30*CB*5.0*C8*CB 
EQAB-CA— EQ8 

E0P-0.03MQ.04CB 

■CONTINUE 

tPlL.EO.l .ANO. EQA8.LT. 0.0)  EOAB-O.O 
IP( L.EQ .2 . ANO.SQAB  *GT .0.0 ) CQAB*0.0 


SUBROUTINE  5ACTIVIXA.X3.ACTA. ACTS. VA.YB.VAA.VA8. VBA. VB81 

IMPLICIT  REAL*  3<  A-H.Q-Z ) 

“COMMON/ WORK E/NACT 
BA* 1.1 380 
BET  12*- 0.0 1A2S 

-ffiaiAW* 

«MA*21t.S3 
W -8*36.10 

WMW*18.  016 

VOLW*18 .07 
810CN*46.70 
8 20EN«- 623.30 

A 1 DEN*  l 28 .473 

“A20CN*— 40 1 .68 

OCNTV*! •0*810CN4XB*B20CN*X4**2*A10CN*X4*A20SN*XA**2 
CONW*( DENT Y-XA4WM A-XR*«<48 ) /WM  4 

XMOL  A*( XA/CONW/wMW )*t  000 .0  

'XMOlA*( X8/C0N*/WMW> *1)00.0 

IPINACT.NC.O)  GO  TO  1 
SOT  T A * ( 3. 0* XMOL A ) • *0. 50 

3QTT8*XMOLB**0 .30  _ _ _ _ , ,L-M  . 

“COEEA*-2.0*ALPM*SQfTX/Tin'3T^'A'k30T5rA)*47tf*BtTli?*XMCCT/3f73 
C0EE8— ALPH*SQTTB/1  1 .0*30  TT3I  *2.0*R  CT134XM0L8 
SOLM*tO.O*«C-S> 

I E(  30TT  A. L E . SQLM ) GO  TO  2 _ 

GA*8.0*8ET1  2/3.0  - 1.  0* ALP-*/30T TA/ ( l . 0*8A*SQTTA ) *42 
V 40*  1 .0  *G  A *XMQL  A 
GO  TO  3 

V AO*  1.0  ..  _ . . _ . _ 


mat  if  BBSS  QiuLinr  mAcricA£i*tf 

QfgSt  JMBCSHE)  TO  COO 


UU  uo 


I F(  SQTTB.LE  .SOL*  ) SO  ro  * 

<58*2. 0*85  T1  3 - *LPM/2.0/SQTT8/(  1 .0*SQTTB)**2 
Y80*  I «0*G8*XM0LB 

GO  TO  5 

VBOal.O 
60  TO  5 
CONTINUE 

JOTAL»XNOLB*3.  0AXM0L4 • 

50T0T« <>OT*U )*  *J  .Uo 

C0EPA*-2.0*ALPH*SaT0r/Il .0*8A«SOTOT  >*S.0*BET1  2*XMCLA/3.0 

I »2.0»<BSTI2»2.0»3CTn>*XM3LB/3.0 

" C0EF8»-ALPH*SQTOT/n  .0  *SO  TOT ) a 2 .0  *SET  1 3 *XMQL8 
I ♦CBET12*2.0*BET13>*XM0LA 
GAaB.0*8ETl2/3.0  - 3. 0*4LPH/SaT0T/< 1 .0*8A*S0T0T) **2 

68*2 .0*  OCT l 3 - ALPH/2 .0/50TOT/< 1 .Q.SQTOT ) »*2 

~ T AO*  t . 0*6A  #XNOLA 
Y80- 1 .0 *G8*XM0L8 
CONTINUE 

GAMAaOEKP(CO*PA ) 

~«f»B«OExVTC0e^T 

ACTA«XMOLA*GAMA 

ACT8aXM0L8*GAM8 

YAaYAO/CONW/VOL*  

Y8»VB0/C0N*/VCL» 

CONMTRaCONV **m**conw*vol* 

V88«6B* I 000 .O/CONaTP  -Y3* ( 31 0BN—MM8*A20EN*2 .0*XB )/WM«/CON« 
YAAaQ.Q 


YBAaO.O 

RETURN 

ENO 


SUBROUTINE  SPATPC T I . T 2. CO. XA. XB. YA. Y8. YAA .YA8.YB A . VBB I 
IMPLICIT  PEAL*8( A-m.O-ZI 

•mr  »-rr3  w;i6/ysrjT.a 

TR2»T2*XA/C0 

TR3al.0-Tl-TR2 

YAal  .S*RTP«<-1 .0-TR2I 

~9  AAa- r.S*RTPaT2/CJ/YA 
YBaRTPa  <*.0-2.0»TP3) 

YBA*2.0*RTP*T2/C0/Y3 

Y A8*0 «0 

TBcraaro 

RETURN 

ENO 


SUBROUTINE  S0IP1CA.C3.YA.YB.YAA.VAB.YBA.YB8 ) 

"twpci  c rf~ffrAir»  * t a - h . r-z i — — : 

XAaC API  000.0 
XBaCB* 1000.0 
YAaO. 606*10  .0**1  -S  > 

~%MK0MaS6.  I 0 

HMH20* I 8. 016 

VISCaO. 8898 1 2*0.073331 *X9*9.0 16»01*XB**2 

1P1XR.LE.  3. 866 ) yISC«0.3»37  0»0.090AS49*XE*0.01I96-32*; 

l¥\  xBTGT  . f.  333  I VT  5 C * >,♦ JaaSO- l . 5T*B  2 * xS + o.  if  967  /*xS*1 

OENSat .0136*0. 03S78*X3 

IP(  X8.LE.S.0I  0CNS-0.1971 *0.0 3939* XB 

C«»  C OENS-CB  **MXOH» /WMH23 

~CT»C  *♦  2.0*  CB  " • 

YBa7.3S8*I 0.0* *( -8 > *293. 16*CT/C*/VI SC 
YAAaO.O 

YA8-0.0 

~Y8*-0.<5 

VB8-0.0 

RETURN 

ENO 

.SUBROUTINE  SOI  PP1X  A.  X3  . TA  44.  T AAO.Y  ABB  . YSAA«.YSA8  . YBB91.. 
IMPLICIT  REAL*3< A-M.O-ZI 
VAAAaO.O 
YAABaO.O 

jrA88.aQ.il 


1*1#  MOB  is  BBS!  QUALITY  FVtCTI 


uu 


■ 


YBA A*0.  0 
YSABaQ.O 
Y 8B 8*0.0 

.RETURN 

ENO 


IMPLICIT  REAL**! A-H.O-Z) 

COMMON/ MORKE/SGMO.SGSO. SSSO.SMLIM 
COMMON/wORKH/TOT.L 

_r.T»  TOT ; 

TTMl-TT-l.O 

EMTMlaO.O 

eSTMfaO.O 


IPIEM.GT.EMLIM)  EM  TM 1 a EM* ATTN  I 
I PI  ES.GT.O.O)  ESTMl«5S**TT'Il 
IPIEG.GT.O.O)  EGTMl a«G**TTM I 

_S  1 6a  SGMOAEM  TM1  *CM>SG5fli£Sni  *.estafi.6(tfE6*e6If«J 

SGOMaSGMO*TT*EMTMI 

SGOSaSGSO«TT*6STMI 

RETURN 


SU6R0UT IN E RKAY I CELT. E0S1L. EPZN • EP ZNO . AM . AMOCUR . AMO CA  »RM»  RS. 

1 RMOCUR.RSOCUA.RMOCA.RSOCA) 

I.MPL  ICIT  R SAL* 91  A^H , 3-Z ) _ — 

COMMON/w0RKC/AO.TIM«.RMO.RSO.4xt.E?I*EPG( .E^NUCsTEPNUCM ,E OUI  ‘ 
X EOUMAX.EPLIM.NQZN.MM. J 

COMMON/MORKM/TOT.L 


S«2. 0/3.0 
EPCaEPG I 
RSOaRMO 

_RR TALaRXI 

ATOTa*  0*  ( I l.O-EPSOD/I  i.O-SPt  ) > **P 
IP(E®ZN.LE.EPLIM)  S®ZN««PLIM 
IPIEPZNO.LE.EPLIM)  EPZNOaEPLIM 


EPSaEPZNO 

IP(L.EQ.t)  EPSaEPZNO«-*PG 
IPIL.EO.2)  EPMaEPZN+EPG 

_ EPMfiaEEM*  «P 

EPSO«EPS**S 

EPMS«EPMP*EPSQ 

EPMPRNaEPMP/EPMS 


NEPMaO 

NEPS«0 

IPIL.EO.2)  GO  TO  4 

IPIEPSPRN.GE.EONUCS ) SO  TO  7 
IPIEQUIAO.LE.l .00-10)  EOUI ABa 1.00-1 0 
IPIEOUIAB.GE.EOUMAX)  EOUI ABaEOUMAX 
EPSPRNaEPNUCS* (EOUI A3/EQUMAX )a*Ng®S 
“GO  TO  7 
CONT INUE 

IP! NOZN.EO.O ) GO  TO  6 

I PI EOUPRN.LE.  1.00-10)  EOUPRNal .00-10 
EPMPRNaEPNUCM*EOUPRN**NEPM 

GO  _T0  7 

CONTINUE 
IPIJ-I)  It. IS. 12 


EPMMTXal  ,'j-EPSO 

J£i  ERME&tt  xLE.»  E£»« 

GO  TO  7 


IPIEPMPRN.6E. EPNUCM)  SO  TO  7 
EOUPRNa 1 • O-OAB  SI E QU I AS ) /! OUM AX 
_IPieOUPRN.LE.  1.00-13)  EOUPflNal^i 
EPMPRNaEPNUCM*EOUPRN**NEPH 
CONTINUE 
AMaATOT*EPMPRN 

AS*ATOT*EPSPRN 

“R'NaPMO*  ( A TO  f-AM ) /6U6G(  A TO  T/ AM  i 
RSIaRSO •( A TOT— AS ) /OLOGI ATOT/AS) 

R S» 1 .0/1 l .O/AS/RxTALM.O/RSl I 
.OTaOELT  . 


y*  war  is  Bast  quality  practicable'' 
18QI  QflDC  AIWLSHAD  TO  ClQ  


tPIMM.LE.il  OTaO.O 
ASTAR»9.l  S/2.0/96487.0 

oel« ia.si «o r/2.0 

_OELT  taASTARAOT/2,g 

del taw-oel*rs" 

P0R0*-2.0/3.0/<  l.O-E’I  I 
ATOTIaATOT APORQAOELT  t 
ATOTAaA TO T APQRC*DE L T A 
"Apocupa  (POpa-p*ePSBPN/;pV)'#AM*oELTr 
AMOCAat  POROaSaEPSFRN/EPS  » a AMAQELTA 
ASXa(PQBO*P»ePMBBN/£OM) AASAOELT I 

ASA«JPORa-SaePPFPN/*2.52*AS*Q.Et,TA 

Tl  MC I a | o'/i  A TOT -A  M l -l . 0 / AT  OT/  OLoGT  Af  OT/ AMI 
RMC2»I.0/AM/0L06(  ATOT/AMI-l.O/I  ATOT-AM) 
RMOCURaRMCl  *ATOT  I *RMC2 AAMOCUR 
RMOCAaRMCl  AATOTA.-PMC2AAMOCA 

"ASCOal  «0/t  ATOT-AS  >-l  « O/A TO T/OLOGl  AY(3f  7a?P 
ASCUal .0/ < ATOT-AS >-l .O/AS  /DLOG ( ATOT/ AS ) 
RSC l*ASCO A< l .0-RS/AS/9XT4LI 
♦ASC0*REAAT3T/RXTA 


l • ATOT I ARSC2AAS  I 
RSOCA  mRSCIAATOTAa  RSC2«ASA 
RETURN 

ENO  _ • 

“SUBROUTINE  !f«  £ lL”<  tffcTVj'ix ) 

IMPLICIT  REALMS ( A— M»  0-2 ) 

DIMENSION  C ( 6 • 93 ) • 1M(  13  > • 0 ( 93  ) « EP( 93 ) 
COMMON/WORKB/N.NJ 

Common/ Work o/C 

COMMON/mORKG/OM.O.EP.H .VOLUME .EPGICC. JPLG 

COMMON/ vORkh/TOT , L 

NJMIaNJ-l 

SUMQMaO.O 

SUMOSaO.O 

00  I jaNJMN  .NJM1  

"SUMQMaSUMQM^O  .S  *M*T7'«TJT'f  3**  ( J + iTl 
SUNOS*  SUM OS* 0. S AHA  I 0 1 J I TO  I J>  t ) I 
CONTINUE 

NJSML*NJMN+$UMQS/H/3( n jmn  > *0 


»ni:Pt.n^r>w»ri.iif:M’ri:yii].iMRiwi:wLT,ip.,,TT?iMP»rTii,TnTi^ 


NJSHMaN  icmi 

IP<OM<njSMLI.GT.O.OI  so  to  2 
NJSWMa2.QASUMQM/H/QM( NJMN)  AN JMN  »0.* 
"CONTINUE  * 

00  3 JaNJMN .NJSAL 
OIJ)aO(NJMN) 

MM 


2PK 


ZTTfYl  *1 J 1 CPlTTi 


t -NJMN} 
CONTINUE 
IFIJ.LT.NJI 


upintr-n* 


CO.JIaCI  3.NJI 
CIA. JI*C(4.NJI 

CIl .J)*C(l .Nj)  

"CO.JIaCI  5.NJI 

CALL  SCONOI C (3.JI.CIA.J). CONONJ .CONOA .C 0NO8 I 

C(2.J)*C(2.NJ)-C(1 * J I AOPLO  AT I J— NJ I ah/CONONJ/E  P<  J I AA I 1.0*T0T> 
CI6.  J)aC(6.NJ| 


CONTINUE 

VCL  UMEa VOLUME-MAOFLOA T I N J 3 AL -N Jl 
NJaNJSML 

"RETURN 


wn  f^at  is  BUST  ^UAiim  PiACIIiCAflU 

jaoa  (wry  m—ttsiflp  xodoo  — - 
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